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During  this  period,  the  central  technical  milestone  of  the 
program  was  achieved:  the  free  electron  laser  was  made  to 
oscillate  for  the  first  time  on  a  storage  ring  device. 


Laser  operation  on  ACO  proved  to  be  moderately  reliable  and 
easy  to  obtain.  We  were  therefore  able  to  extract  an  enormous 
volume  of  information  on  the  operation  of  the  system.  In  contrast 
to  the  other  operating  FELs  ,  the  steady  availability  and  good 
reproducibility  of  operation  has  permitted  the  investigation  of 
questions  important  to  the  future  of  the  field  but  which  are 
difficult  to  observe  in  present  machines. 


The  laser  operation  was  discovered  to  exhibit  a  self 
pulsing  behavior  under  conditions  of  free  running  operation.  We 
have  measured  the  dependence  of  the  pulsing  phenomena  on  the 
current,  and  the  cavity  length  detuning.  We  discovered  that  the 
pulses  could  be  modulated  by  driving  the  electron  beam  in  either 
transverse  position  or  RF  frequency.  Theoretically,  the 
dependence  of  the  peak  power  on  the  driving  frequency  provided  an 
unambiguous  identification  of  the  mechanism  of  the  pulsing 
pne.nomenon.  We  report  these  results  in  the  papers: 

First  operation  of  a  storage  ring  free  electron 
laser,  SPIZ  Vol  453  (1934) 

Results  of  the  free  electron  laser  oscillation 
experiments  on  the  ACO  storage  ring,  J.  Physique  45, 
(1984) 

First  Operation  of  a  storage-ring  f ree-electron 
laser,  Phys .  Rev.  Lett  51,  (  19  8  3' 

A  further  discovery  during  the  contract  was  the  existence 
of  transverse  mode  coupling  due  to  the  finite  cross  section  of 
the  electron  beam.  We  devised  a  technique  to  measure  this  effect 
which  relies  on  the  interference  of  the  excited  modes  with  the 
fundamental.  We  were  able  to  show  that  the  higher  order  modes 
excited  in  this  way  could  be  predicted  by  the  theory  in  a 
particularly  simple  case,  and  that  the  mode  evolution  had  the 
appropriate  behavior  as  a  function  of  longitudinal  position. 

These  results,  both  experimental  and  theoretical,  are  reported 
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Measurement  of  stimulated  transverse  mode  mixing  in 
a  free  electron  laser,  IEEE  J.  Quant.  Elect.  (1985) 

Transverse  Mode  dynamics  in  a  free  electron  laser, 
Appl .  Phys  S33,  (1984) 

A  series  of  other  issues  were  addressed,  including  the 
effects  of  divergence  on  the  Madey  theorem.  We  were  able  to  show 
that  indeed  the  divergence  of  the  optical  mode  shifts  the  gain 
and  spontaneous  spectra  in  a  systematic  way.  These  results  are 
described  in: 

'Measurement  of  the  violation  of  the  Madey  Theorem 
induced  by  a  diverging  wave,  IEEE  J.  Quant.  Elect 
(1935) 

Observation  of  the  diffraction  induced  violation  of 
the  Madey  theorem;  in  the  Proceedings  of  the  1984 
Castelgandalfo  Conference,  Nucl .  Instr.  &  Meth. 

A237,  (1985! 

The  analysis  of  the  optical  klystron  experiments,  including  the 
application  to  energy  spread  measurement  in  a  storage  ring,  are 

in  ■ 

Optical  klystron  experiments  for  the  ACO  storage 
ring  free  electron  laser,'  Appl.  Phys.  B34  (1984) 

The  bunch  lengthening  analysis  presented  in 

Characterization  of  free  electron  laser  bunch 
lengthening  on  the  ACO  storage  ring,'  Appl.  Phys.  3 
; 1984  )  , 

completes  the  comparison  between  the  theory  and  the  experiment, 
and  describes  a  unique  experimental  approach  for  measuring  the 
momentum  compaction  factor  in  a  storage  ring.  The  other 
publications  for  the  period  include: 

Inhomogeneous  Broadening  calculation  with  a  single 
integral;  HEPL  Technical  Note  TN-34-2 

Realization  of  a  variable  aperture  diaphragm  working 
in  ultra  high  vacuum.  Applied  Optics  23,  (1984) 

Copies  of  these  papers  are  attached. 


♦V  V 

Realization  of  a  variable  aperture  diaphragm  working 

in  ultra  high  vacuum 

M.  Velghe,  D.  A.  G.  Deacon,  and  J.  M.  Ortega 

An  adjustable  iris  has  lieen  develuped  to  operate  in  the  10'  '"Torr  vacuum  environment  of  a  storage  ring. 
It  was  baked  out  at  ,'100°('  and  has  hoen  used  more  than  ItH)  times.  Measurements  have  been  carried  out 
using  the  spontaneous  emission  of  an  undulator  set  on  the  storage  ring.  They  show  that  the  iris  could  be 
centered  and  adjusted  accurately  from  outside  the  vacuum  chamber  and  used  as  a  mode  selector  for  free 
electron  lasers. 

I.  Introduction  It.  Description  of  the  Device 


Laser  transverse  mode  selection  by  an  iris  placed 
inside  the  optical  cavity  has  been  extensively  studied 
in  the  past  (Ref.  1  and  references  therein).  Since  the 
optical  gain  of  a  free  electron  laser  ( FKL)  may  develop 
in  several  transverse  modes,- the  Orsay  FKL  group  de¬ 
cided  to  install  an  iris  inside  the  laser  optical  cavity. 
Although  laser  oscillation  was  finally  obtained1’  on  the 
TKMun  mode  without  the  help  of  that  iris  due  to  the 
extremely  low  value  of  the  gain,'  future  high  gain  sys¬ 
tems  are  expected  to  operate  multimode.  It  is  useful 
to  look  at  the  technical  solutions  which  have  led  to  an 
ultra  high  vacuum  variable  aperture  diaphragm  and  to 
,  its  practical  operation.  This  is  the  purpose  of  this 
|  0  raper. 

Due  to  the  weakness  of  the  optical  gain  in  the  Orsay 
experiment  on  the  ACO  storage  ring  it  was  impossible 
to  insert  any  optical  element,  such  as  a  window,  between 
the  laser  medium  (high  energy  electrons  passing 
through  the  undulator)  and  the  cavity  mirrors.  For  that 
reason  the  iris,  as  well  ns  the  cavity  mirrors,  has  to  he 
placed  in  a  vacuum  chamber.  The  pressure  inside  this 
chamber  is  then  necessarily  as  low  as  the  pressure  inside 
the  storage  ring,  that  is,  between  10~‘Jand  l()-ll,Torr, 
since  the  two  chambers  communicate  (Fig.  1). 

The  technical  problems  posed  by  the  requirements 
for  three  fine  mechanical  motions  (two  translations  and 
one  rotational  for  the  iris  aperture)  are  severe  in  this 
environment.  In  addit  ion  the  whole  chamber  has  to  he 
baked  to  at  least  200°C  (for  degassing)  to  reach  the 
desired  vacuum.  No  such  device  is  commercially 
available  (existing  iris  diaphragms  are  designed  for 
pressures  above  10~7  Torr). 

When  this  work  was  Hone  M.  Velghe  was  with  Universite  He 
Paris  Slid,  l.ahoratoire  He  Photophysique  Moleculaire,  91-105  Orsay, 
Km  bit;  he  is  now  with  Stanford  University,  High  Knergy  Physics 
I.ahoratory,  Stanford,  California  91:105.  D.  A.  (*.  Deacon  is  with 
Deaton  Hoearch,  75-1  Diimartlinc  Way,  Sunnyvale.  California  9-1087, 
and  I.  M.  Ortega  is  wiih  Krole  Supcrieure  He  Physique  el  Chimie,  10 
rue  Vauquetin,  75291  Paris  CEDKX  05,  Frame. 
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We  have  designed  and  tested  a  device  which  satisfied 
the  Ihree  following  requirements  (Fig.  2): 

( 1 )  Given  the  limited  space  available  in  the  chamber 
the  iris  has  to  he  movable  in  a  direction  perpendicular 
to  the  light  beam  axis  so  that  it  can  be  replaced  by  a 
mirror.  'Phis  mirror,  oriented  at  45°  with  respect  to  the 
beam  axis,  is  used  to  extract  the  light  when  the  optical 
cavity  is  not  in  use.  Therefore,  the  iris  and  the  mirror 
are  mounted  on  the  same  axis,  and  switching  from  one 
to  another  is  made  by  a  translation  stage  operated  from 
outside  the  chamber. 

(2)  The  iris  centering  has  to  he  adjusted  mechanically 
over  a  distance  typically  of  ±5  mm  in  any  direction  in 
the  plane  perpendicular  to  the  light  beam  axis.  Vertical 
adjustment  is  ensured  by  the  micrometer  stage  men¬ 
tioned  above.  The  horizontal  displacement  is  realized 
by  a  second  stage  mounted  perpendicularly  to  the  first 
(Fig.  2).  The  two  displacements  are  allowed  by  the 
same  ultra  high  vacuum  bellows  (Fig.  2,  bellows  1). 

CD  'Fhe  aperture  of  the  iris  is  controlled  from  outside 
the  vacuum  chamber.  The  rotational  motion  needed 
to  change  the  iris  dimension  is  accomplished  by  means 
of  a  1 5-ram  vertical  linear  translation,  given  by  a  third 
translation  stage,  allowing  any  aperture  ranging  from 
1.6  to  38  mm  in  diameter  (Fig.  2,  bellows  2). 

The  mechanical  assembly  of  these  various  parts  is 
shown  in  Fig.  3.  The  whole  chamber  may  be  baked  at 
300°C.  The  ultra  high  vacuum  requirements  mean  that 
grease  cannot  he  used  to  lubricate  mechanical  parts. 
We  shall  now  discuss  the  solution  we  used  to  move  the 
delicate  pieces  of  a  diaphragm  under  ultra  high 
vacuum. 

III.  Lubrication  of  the  Diaphragm 

Commercial  irises  are  usually  made  of  aluminum  or 
brass  and  arc  painted  or  coated  with  Teflon  for  lubri¬ 
cation.  None  of  these  materials  is  usuable  in  ultra  high 
vacuum  because  of  their  outgassing  and  structural 
changes  which  occur  during  bake  out.  We  have  fabri¬ 
cated  a  greaseless  stainless  steel  iris  diaphragm. 

Initial  tests  at  4  X  10-10  Torr  terminated  abruptly 
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Fig.  1.  General  setup  of  the  undulator  experiment  mi  the  MX) 
storage  ring  FEL. 


Fig.  2.  Design  of  the  iris  rear  mirror  device  within  its  vacuum 
chamber. 


when  a  spring  clip  which  held  together  the  two  movable 
posts  fell  off.  The  bake  out  temperature  of  300°C 
probably  annealed  the  piece  which  lost  its  elasticity. 
Upon  extracting  and  inspecting  the  iris  we  also  observed 
that  the  two  moving  parts  that  rotate  coaxially  showed 
some  mutual  scratching.  Since  lubrication  was  not 
possible  we  had  to  specify  a  certain  clearance  between 
the  two  pieces.  One  of  the  two  parts  undergoes  a  strain 
perpendicular  to  its  rotation  axis  and,  therefore,  rubs 
strongly  against  the  other  piece. 

Device  modification  was  necessary  to  avoid  these 
problems.  We  used  a  dry  lubrication  with  very  low 
outgassing.  In  this  process  a  jet  of  a  molybdenum  di¬ 
sulfide  plasma  is  projected  onto  both  moving  parts. 
The  plasma  is  adsorbed  by  the  surfaces  to  a  depth  of  ~1 
pm.  In  addition  the  spring  clip  mentioned  previously 
was  replaced  by  a  threaded  ring  which  holds  the  two 
moving  parts  more  firmly  and  which  can  be  baked  out 
at  300°  C.  Finally  the  stainless  lamellae  of  the  iris, 
which  bad  developed  little  burrs  due  to  embossing,  were 
carefully  polished  with  a  soft  stone. 

Alter  these  improvements  the  iris  diaphragm  has 
worked  properly  for  more  than  100  cycles  (opening  - 
closing)  in  a  vacuum  of~5X  10“ 1,1  Torr  after  a  bake  out 

at  ;ioo°c. 


IV.  Experiment 

To  test  the  iris  operation  we  measured  its  effect  on 
the  light  stored  in  the  optical  cavity.  The  light  was 
emitted  by  the  electrons  of  the  AGO  storage  ring  passing 
through  an  undulator’6  (see  Fig.  1).  We  centered  the 
iris  in  the  colored  rings  emitted  by  the  undulator  (Fig. 
4).  The  fine  adjustment  was  made  by  maximizing  the 
intensity  transmitted  by  the  output  mirror  of  the  cavity 
of  the  resonant  wavelength  (A  =  6300  A  in  our  case). 
We  then  measured  the  intensity  of  the  light  at  this 
wavelength  as  a  function  of  the  iris  aperture  (Fig.  5). 
The  measured  intensity  drops  sharply  when  closing  the 
diagram  below  a  diameter  of  ~4  mm  due  to  the  very 
high  quality  of  the  optical  cavity  whose  losses7  are  very 
low  (1.2  X  10':'  in  this  experiment).  Since  the  light  is 
stored  for  ~ll):’  round  trips  in  the  cavity  a  small  loss  at 
each  pass  on  the  diaphragm  produces  a  large  variation 
in  the  resulting  intensity. 

The  spontaneous  emission  in  the  undulator  of  the 
electron  beam  at  235  MeV  diverges  with  a  half-angle 

„  /  2.8A  y/» 

0^ - 1  ^  0.66  mrad,  (1) 

where  N  =  17  is  the  number  of  periods  and  Aq  =  7.8  cm 
is  the  magnet  period.  This  angle  is  only  a  factor  of  2.7 
larger  than  the  divergence  angle  \/ttWq  of  the  TEM(K) 
mode.  In  these  conditions  the  spontaneous  emission 
is  expected  to  be  stored  on  the  first  three  of  four  cylin¬ 
drical  modes  of  the  cavity.8 

The  intensity  transmitted  by  the  system  can  be  easily 
calculated  for  various  cylindrical  transverse  modes  if 
one  makes  the  assumption  that  the  form  of  the  modes 
is  unchanged  when  closing  the  diaphragm.  This  as¬ 
sumption  is  probably  not  too  bad  in  our  case  since  the 
apertures  of  the  iris  for  which  the  stored  intensity 
vanishes  is  much  larger  than  the  beam  waist  at  the  iris 
location  (a  factor  of  I  5  in  our  case).  The  loss  at  each 
passage  through  the  iris  is  then  found  by  integrating  the 
mode  intensity  in  the  space  surrounding  the  aperture. 
Since  the  iris  is  only  separated  from  the  cavity  mirror 
by  ~25  cm  this  power  reduction  occurs  only  once  each 
round  trip.  If  we  include  the  power  lost  at  the  mirrors 
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Fig.  4.  Spontaneous  emission  produced  hy  235-MeV  electrons 
traveling  through  the  undulator. 


iiL'tut:  i  ini.  oiuit  patterns.'-1  ThcriTure.  this  approximate 
I  lu-iuA  at.  - 1  i  i  lies  the  data  sal  isl.icl  nr  it  y  and  allows  th.tt 
tile  iri..  has  been  i  entered  ami  operated  properlv. 

V.  Conclusion 

We  have  const rneleil  a  variahle  aperture  iris  capable 
of  operating  in  the  very  difficult  conditions  of  an  ultra 
high  vacuum  environment.  This  device  has  been  useful 
in  the  analysis  of  the  mode  structure  of  the  light  emitted 
by  our  undulator.  Also  it  appears  that  it  will  be  able  to 
serve  as  a  mode  seleetor  on  the  future  storage  ring  free 
electron  lasers  particularly  in  the  VI I V  spectral  range 
where  a  high  vacuum  will  lie  probably  needed,  for  var¬ 
ious  reasons1-7  in  the  whole  optical  resonator. 
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Fig.  5.  Spontaneous  emission  intensity  vs  aperture  of  the  iris. 


and  sum  over  all  the  round  trips,  we  find  that  the  in¬ 
tensity  in  the  lowest-order  mode  is  reduced  from  its 
value  for  the  iris  fully  open  by  the  factor 

/nolul  _ _ r _  _  ^2) 

/oo(*)  I'  +  exp(-2a2/ir2) 

where  a  is  the  iris  aperture,  w  =  1.03  mm  is  the  size  of 
the  mode  at  the  iris,  and  I'  is  the  total  cavity  losses  per 
round  trip.  For  the  TKMio  cylindrical  mode,  the  ex¬ 
pression  is 


‘  u>l )  /  a 

I’  +  It  f  —j  c.xp <— l?or-/it>-) 

Fitting  these  curves  to  the  data  of  Fig.  5,  we  can  see 
that  55%  of  the  output  power  is  contained  in  the  low¬ 
est-order  mode,  and  -13%  of  the  remainder  is  contained 
in  the  THM,,,  within  the  accuracy  of  the  fit. 

This  result  is  in  good  agreement  with  a  previous 
analysis  using  a  mathematical  analysis  of  the  sponta- 
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Abstract 


The  complex  problem  of  calculating  the  inhomogenous 
broadening  in  an  FEL  with  arbitrary  energy  spread,  and  hori¬ 
zontal  and  vertical  angular  spreads  is  reduced  to  a  single  well 
conditioned  integral  which  can  easily  be  performed  on  a  micro¬ 
computer  . 


The  calculation  of  the  gain  in  an  FEL  driven  by  a  "hot" 
beam  in  a  recurring  problem  for  FEL  system  designers.  A  know¬ 
ledge  of  the  gain  reduction  produced  by  a  given  energy  spread 
and  emittance  is  essential  in  designing  short  wavelength  lasers 
or  short  interaction  length  systems  in  which  the  available  gain 
approaches  the  net  cavity  losses. 

The  problem  is  conceptually  simple,  as  it  involves  nothing 
more  than  a  convolution  of  the  single  electron  gain  spectrum  with 
the  effective  energy  distribution  of  the  electron  beam.  No  closed 
form  expression  exists  for  the  convolutions,  and  results  must  be 
obtained  numerically.  However,  the  general  case  involves  separate 
convolution  integrals  for  the  energy  spread,  the  horizontal  angular 
spread,  and  the  vertical  angular  spread.  A  numerical  solution  to  this  prob-- 
lem  involving  three  integration  steps  can  become  impossibly  time  con¬ 
suming  on  the  computer,  and  complex  to  set  up  due  to  the  divergences 
in  the  effective  energy  distribution  of  the  angular  spread,  and  due 
to  the  divergence  of  the  denominator  in  the  gain  expression  at  v=0. 

In  this  paper,  I  present  a  solution  to  these  problems  using  the 
properties  of  the  Fourier  transform,  which  reduces  the  triple  inte¬ 
gration  with  divergences  to  a  single  integration  of  an  analytic  functio 
The  resulting  expression  takes  only  a  few  seconds  to  evaluate  on  a 
microcomputer,  and  can  therefore  be  incorporated  as  a  subroutine  into 
a  larger  progream  which  evaluates  the  performance  of  a  proposed  design. 
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I  assume  the  "hot 


electron  beam  can  be  approximated  by  a 


Gaussian  distribution  in  energy,  angle,  and  position.  The  resonance 
parameter,  \> ,  of  an  arbitrary  electon  is: 

»  s  27 rfJ  l^l  -  £<•»•■*>*) 

for  a  planar  undulator  where  N  is  the  number  of  periods  in  the 

undulator,  X  and  X  are  the  undulator  and  light  wavelengths, 

and  y nc“  and  K  are  the  electron  energy  and  trajectory  para- 
A  2 

meter  K  =  2  a  eD/mc  X  .  If  we  take  as  our  reference  the  resonance 

o 

parameter  x  of  an  axial  electron  at  the  central  energy  yq  / 


X  2  T-TiN 


(2) 


we  find,  for  v  <<  N, 

v  -  x  +  H-TTtJ 

•  f  +  r 

where  5y  =  y  -y  .  The  gain  curve  of  the  hot  beam  must  be 
expressed  in  terms  of  x. 


The  numerical  distribution  P(5y,8x,9^)  of  electrons  is 
assumed  to  be  separable. 


P^.y.v) 


The  average  gain  is  then 


“’HI 

p  f  \ 

--  ct^P^)  Jte^coj  (5) 


which  can  be  reduced  to  the  convolution  integrals 


<6->  =  1(2)  *  A*(*)  •*  -V*)  *  G-(v(*>) 


where 


EM  *  ==;  e 

Jzttct 


- 


GY 

<r  =  V-7T/V 


Ac*)  = 


-  x/cC 


hr  ax. 


a  ‘t-TTtjfC 

i+  f 


and  0 (x)  is  the  step  function  0 (x  >  0)  =  1;  0 (x  <  0)  =  0. 


Let  us  neglect  the  constant  in  G(v),  and  use  the  analytic  expression 


for  the  small  signal  gain  spectrum,  including  the  phase  shift 
of  the  slowly  varying  wave 


G-M  - 


c<n  V  -  — 


y>  ^  \7  y  -  %  5,'m 


-  I 


(9) 


It  is  clear  from  (8)  and  (9)  that  care  must  be  taken  in  performing 
the  three  integrals  in  (6)  to  correctly  handle  the  divergent  deno¬ 
minators  . 

Equation  (6)  can  most  simply  be  handled  by  taking  the  Fourier 
transform,  and  using  the  convolution-product  theorem: 


<  6-(x)> 


U  e 


£(<*) Ay  660 


(10) 


-s 


where  E(ct),  (a),  etc.',  are  the  Fourier  transform  of  E(x),  A(x), 
etc. 


tU)  =  j<fx  tU)  e'** 


-<*l<ryz 
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(ID 


A  (*)  --  0 


-t  i  cut 


■% 


-  (l  +  A1*1')  c*p{'  i  (12) 


The  transform  of  G(v)  can  most  simply  be  obtained  if  we  use 
the  expression: 


5 
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Ecuation  (15)  permits  the  calculation  of  both  the  real  and 
imaginary  parts  of  the  inhomogeneous  broadened  gain  spectrum  as 
a  function  of  x,  which  was  defined  in  (2)  as  the  resonance  para¬ 
meters  of  an  ideal  filamentary  beam.  It  is  valid  in  the  small 
signal  regime  for  an  undulator  with  zero  taper.  In  the  large  signal 
regime,  the  inhomogeneous  broadening  is  of  less  interest  since  the 
gain  is  above  threshold,  but  one  can  obtain  the  gain  reduction  by 
convolving  the  numerically  obtained  gain  spectrum  with  the  beam 
function  B (x) . 
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(16) 


The  other  effects  which  reduce  the  gain  are  independent  of 
the  inhomogeneous  broadening  (15).  The  longitudinal  effects,  of 
course,  do  not  couple  with  the  inhomogeneous  broadening  term,  and 
have  been  calculated  in  [1]  and  [2].  The  same  is  true  for  the  trans¬ 
verse  effects,  which  can  be  included  as  a  multiplicative  filling 
factor  [3]  for  low  divergence  optical  and  electron  beam.  The  optical 
mode  divergence  will  at  some  value  reduce  and  distort  [  4]  the  gain 
curve;  the  effects  of  large  electron  beam  divergence  are  also  un¬ 
desirable,  although  they  have  not  been  calculated. 


T 
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The  effects  of  the  quadratic  dependence  of  the  unaulator  field 
on  one  of  the  transverse  dimensions  could  have  been  included  in 
(15)  in  an  approximate  way  (the  position  distribution  is  identical 
to  (8)  with 


CL 


(ZTT)Z  h}<T* 

0+  i')  a: 


(17) 


but  this  effect  should  really  be  handled  in  a  transverse  mode 
calculation  such  as  [3].  The  assumption  of  separability  made  in 
[4]  is  valid  for  most  systems,  and  produces  a  worst-case  estimate 
for  coupled  systems. 
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Appendix 


These  tables  list  the  values  of  <  G (x  )  >/G  and  x  for 

o  max  o 

values  of  a,  a  and  a  between  0  and  5.  <  G  >/G  for  each 

x  y  max 

entry  is  listed  above  xq,  which  is  the  resonance  parameter  at 

the  peak  of  the  inhomogeneously  broadened  gain  curve.  The 

normalization  G  =  .0675  is  the  maximum  value  of  the  gain 
max  ^ 

curve  attained  with  no  broadening  effects. 
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ABSTRACT 


This  paper  summarizes  the  basic  low-current  laser-induced  bunch-lengthening 
measurements  that  have  been  made  on  the  ACO  Storage-Ring  Free-Electron  Laser 
(SRFEL).  The  measurements  provide  verification  of  both  the  functional  dependence 
and  absolute  magnitude  of  SRFEL  theoretical  models.  The  method  of  measure¬ 
ment,  which  is  is  explained,  exploits  frequency-domain  techniques  and  is  capable  of 
accuracies  comparable  to  those  of  a  streak  camara.  The  measurements  are  in  good 
agreement  with  existing  SRFEL  theory  and  provide  an  important  base  for  future 
work. 


PACS:  42.55T,  29.20,  41.80D 
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INTRODUCTION 

Operation  of  a  free-electron  laser  as  part  of  a  storage  ring  [SRFEL]  will  cause 
an  increase  in  the  nominal  energy  spread  of  the  electron  bunch  if  the  optical 
ph  ase  of  the  electrons  is  uncorrelated  from  pass  to  pass.  This  increase  in  energy 
spread  results  in  an  enlargement  of  the  bunch  length.  The  increased  bunch  length 
and  energy  spread  cause  a  reduction  in  gain  and  can  lead  to  the  saturation  of 
the  SRFEL  for  low  gain  systems.  As  part  of  the  LURE/Stanford  ACO-SRFEL 
project  this  FEL-induced  bunch  lengthening  was  measured  for  the  first  time  on  the 
superconducting  undulator  mounted  on  ACO  in  1981  [1,2].  In  order  to  measure  this 
FEL-induced  bunch  lengthening,  we  have  developed  a  frequency-domain  technique 
which  measures  induced  bunch-length  changes  to  1  part  in  103  and  absolute  bunch 
lengths  to  a  precision  of  lCo  as  established  by  data  taken  on  stored  electron  beams. 

In  this  paper  we  present  an  outline  of  the  measurement  theory,  experimental 
technique,  and  a  complete  description  of  the  basic,  low-current  results  of  the  laser- 
induced  bunch  lengthening  obtained  on  ACO.  We  present  results  which  provide  a 
graphic  demonstration  of  the  scaling  of  FEL  interaction,  of  the  storage  ring  scaling, 
and  absolute  magnitude  of  the  FEL-induced  bunch  lengthening  as  a  function  of  laser 
intensity.  An  appendix  is  also  included  which  summarizes  the  techniques  used  for 
measuring  the  storage  ring  parameters  necessary  to  compare  the  data  with  theory. 

At  high  currents  on  ACO  electron  beam  interactions  with  the  environment 
cause  anomalous  bunch  lengthening  [3].  This  causes  strong  nonlinear,  anomalous 
coupling  with  the  FEL  altering  both  the  quantitative  and  qualitative  heating  effects. 
The  nature  and  scaling  of  the  FEL  interaction  at  high  currents  on  ACO,  which  is 
indicative  of  fine  structure  on  the  electron  beam,  have  been  reported  elsewhere  [4- 
7].  We  will  restrict  this  paper  to  those  results  which  apply  to  the  verification  of 
existing  SRFEL  theory. 


Measurement  Theory  and  Procedure 

All  measurements  presented  in  this  paper  were  made  w'hile  operating  the  SRFEL 
as  an  amplifier  to  an  external  laser.  This  is  the  same  configuration  used  for  the 
gain  measurments  [8,9].  Figure  1  is  a  schematic  diagram  of  the  ACO-SRFEL  and 
external  laser.  An  argon-ion  laser,  D,  is  aligned  coaxial  to  the  electron  beam  in 
the  undulator  through  a  mode  matching  telescope  and  optical  transport  system,  E. 
Synchrotron  light  from  a  bending  magnet,  F,  is  then  detected  on  a  fast  photodiode. 
The  bending  magnet  chosen  was  to  prevent  any  spurious  light  from  the  laser  or 
undulator  from  reaching  the  photodiode. 

Figure  2  is  a  schematic  of  the  electronic-detection  system  used  for  the  bunch¬ 
lengthening  measurements.  The  synchrotron  light,  focused  by  cylindrical  mirrors, 
is  detected  by  a  fast  photodiode  (B  AH  OC3002LN,  or  ITL3181)  and  sent  via  a  low 
dispersion  cable  (RG  331)  to  an  RF-spectrum  analyzer  (AILTECH  Model  757)  which 
is  interfaced  directly  with  a  computer.  The  interface  between  the  spectrum  analyzer 
and  the  computer  is  used  for  the  aquistion  of  global  spectra.  When  monitoring  the 
bunch  length  changes,  the  analyzer  is  operated  as  a  phase-locked  receiver  tuned  to 
a  harmonic  of  the  orbit  frequency.  The  analyzer  output  can  either  be  input  to  a 
lock-in  amplifier  for  synchronous  detection  at  the  frequency  of  the  laser  chopper,  or 
time  averaged  measuremerfts  can  be  made  using  an  averaging  digital  oscilloscope. 
In  all  cases  the  detected  signal  is  recorded  on  the  data  aquisition  system  together 
with  the  storage  ring  parameters  including  the  RF  voltage,  the  storage  ring  current, 
and  the  output  power  of  the  external  laser. 

After  injection  of  ACO  is  completed  the  undulator  is  closed  and  the  magnetic 
optics  of  the  machine  are  compensated  for  the  focusing  induced  by  the  undulator. 
ACO  is  operated  at  the  coupling  resonance  of  the  betatron  tunes.  This  assures 
nearly  round  beams  for  facilitating  interpretation  of  the  filling  factors.  Corrector 
coils  are  used  to  change  the  orbit  from  the  nominal,  injection  orbit  to  one  which  is 
centered  in  the  undulator  as  established  by  imaging  the  electrons  in  the  undulator. 
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With  the  do  ctron  optics  of  the  machine  established  through  the  undulator,  the  laser 
is  aligned  onto  the  electron  beam.  After  optimization  of  the  magnetic  field  and  laser 
alignment,  bunch  lengthening  data  are  taken. 

Since  the  cumulative  effect  of  the  additional  FEL-induced  energy  spread  is 
the  principal  concern,  the  description  of  the  electrons  in  the  storage  ring  can 
be  restricted  to  the  synchrotron  motion,  including  synchrotron  oscillation  and 
damping,  in  the  longitudinal  dimension.  Several  works  exist  which  give  a  complete 
description  of  storage  ring  theory  [10-12],  In  linearized  form  the  equation  for 
synchrotron  motion  has  the  form  of  a  damped  harmonic-oscillator  driven  by  a  noise 
source  [  1 0] : 


(fl 

dt 2 


2  dc 

H - -Jl  +  fl~  ?  — =  6(q  +  8ef  el. 

T.  dt 


(1) 


In  this  equation  <;  is  the  longitudinal  phase  of  the  electron  measured  with  respect  to 
the  synchronous  or  design  electron  in  units  of  the  orbit  period  T0.  The  incoherent 
energy  damping  time  is  re  ,  and  J?  is  the  angular  synchrotron  frequency  of  the 
electrons.  The  two  terms  on  the  right-hand  side  of  Eqn.  (1)  are  the  driving  terms 
arising  from  fluctuations  as  a  result  of  the  discrete  nature  of  the  synchrotron 
radiaton  6eq  [10],  and  the  interaction  of  the  electrons  with  the  optical  field  in 
the  FEL,  bejel  [13].  Under  idealized  conditions,  the  quantum  fluctuations  give 
rise  to  steady-state  Gaussian  electron  distributions  in  energy-  crt  and  phase  crf  [10], 
The  scaling  factor  between  the  two  distributions  is  the  ratio  of  the  momentum- 
compaction  factor  and  the  angular  synchrotron  frequency 
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The  first-order  interaction  between  the  electrons  and  the  optical  field  bejei  in 
Eqn.  ( 1)  causes  an  increase  of  the  nominal  energy  spread.  In  storage  rings,  since  the 
electrons  are  recirculated,  the  first-order  damping  of  this  increased  energy  spread 
will  occur  from  the  normal  synchrotron  damping  of  the  storage  ring.  (There  is  a 
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possible  FEL-induced  damping  but  it  is  of  higher  order  and  becomes  important  only 
at  optical  intensities  much  greater  than  those  accessible  on  AC'O  [11-16].) 

Generally  the  optical  phase  of  an  electron  in  a SRFEL  will  be  uncorrelated  from 
pass  to  pass  in  the  undulator.  This  will  be  true  unless  deliberate  action  is  taken 
in  the  design  of  the  storage  ring  as  is  the  case  with  the  isochronous-storage-ring 
[17],  Because  of  the  uncorrelated  contributions,  the  FEL  interaction  will  behave  as 
an  additional  stochastic  noise  term  in  Eqn.  (1).  Its  contribution  will  be  scaled  by 
the  energy  damping  time  of  the  ring  [ij.  In  the  small-signal  case  the  Central-Limit 
Theorem  [18]  provides  the  value  of  the  expected  increase  in  the  steady-state  energy 
spread  in  the  presence  of  the  FEL: 

o\—ato  2 +  (3) 

O 

Here.  rr,a  is  the  nominal  relative  energy  spread  of  the  storage  ring  resulting  from 
quantum  fluctuations  [10].  The  optical  phase  of  the  electrons  V  is  averaged  over 
to  obtain  the  second  moment  of  the  FEL  interaction  {(<5c/ef  )2)o-  At  higher  power 
levels,  where  the  FEL  term  in  Eqn.  (1)  is  comparable  to  cr(o,  the  small-signal  case 
is  no  longer  valid.  Under  these  circumstances  the  energy  dependence  of  the  FEL- 
induced  spread  must  be  taken  into  account  as  Renieri  and  others  have  done  for  the 
idealized  SRFEL  operating  as  an  amplifier  and  oscillator  [13,15,19-21]. 

A  major  consequence  from  the  accumulated  energy  spread  is  the  mechanism  of 
gain  saturation  of  the  SRFEL  for  low-gain  systems  [6].  The  gain  may  be  depressed 
and  saturate  in  one  of  two  processes:  (1)  the  increased  energy  spread  of  the  electrons 
fills  the  energy  acceptance  of  the  undulator  reducing  the  number  of  gain  contributing 
particles,  or  (2)  the  increased  energy  spread  will  increase  the  bunch  length,  and 
thereby  decreases  the  charge  density  and  the  gain  which  is  proportional  to  the 
charge  density.  Both  of  these  mechanisms  are  central  to  the  operation  of  SRFEL's, 
and  underscore  the  importance  of  the  bunch-heating  phenomena  [22]. 

The  synchrotron  radiation  from  a  normal  bending  magnent  in  storage  ring  is 
generally  assumed  to  have  no  coherent  or  interference  effects  present  [23j.  The 
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bending  magnet,  may  be  considered  as  an  isolated  element  of  the  system,  and  the 
electron  phases  are  uncorrelated  [21] .  Since  the  emitted  radiation  is  incoherent,  the 
number  of  photons  emitted  through  a  given  solid  angle  of  the  bending  magnent  is 
proportional  to  the  number  of  electrons.  Equivalently,  the  instantaneous  intensity 
of  the  synchrotron  radiation  /„(/)  is  proportional  to  the  charge  density  p(t): 

I,(t)ocp(t).  (-1) 

The  time  dependence  of  the  synchrotron  radiation  is  characterized  by  a  series 
of  pulses  separated  by  a  time  T0/k.  where  T0  is  the  orbit  period  of  the  storage  ring 
(73. 1  nanoseconds  on  AGO)  and  k  is  the  number  of  bunches  stored  in  the  ring.  The 
pulses  will  have  a  temporal  length  aT  which  corresponds  to  the  longitudinal  length 
of  the  electron  pulse.  One  way  of  measuring  the  length  of  the  electron  bunch  is 
to  simply  measure  the  length  of  the  optical  pulse  Is(t)  emitted  by  the  electrons. 
However,  in  practice  unless  a  streak  eamara  is  used,  this  technique  is  incapable 
of  providing  the  precision  required  for  adequate  bunch  length  and  induced  bunch¬ 
lengthening  measurements. 

In  order  to  obtain  the  precision  needed  for  measurement  we  adopted  a  frequency- 
domain  technique  which  measures  the  shape  of  the  Fourier  transform  of  /„(/).  The 
Fourier  transform  of  the  charge  density  and  the  synchrotron  radiation  intensity  are 
also  proportional: 


/«M  «  (5) 

If  a  detector  is  placed  in  the  synchrotron  radiation  from  the  bending  magnet 
the  signal  recorded  /(/)  will  be  proportional  to  the  convolution  of  the  detector 
time  response  to  a  delta  function  and  the  detected  signal:  /(/)  oc  /;(/)  ©/«(/).  In 
order  to  measure  the  pulse  (bunch)  length  the  response  of  the  the  detector  must 
be  deconvolved  from  the  recorded  signal.  This  is  relatively  simple  if  the  signal  and 
detector  response  are  Gaussian,  but  in  practice  detector  response  is  rarely  Gaussian 


[7].  In  the  frequency  domain  the  temporal  convolution  becomes  a  product  and  7a(«j) 
can  be  obtained  quickly  from  Tp)  by  dividing  out  the  detector  response  P(u). 

If  /.,(/)  is  modeled  the  inverse  bunch  length  can  be  estimated  by  fitting  the 
compensated  spectrum  Fpl/Tp)  to  the  apropriate  envelope.  Under  idealized 
conditions  the  longitudinal  electron  distribution  will  be  Gaussian  [10],  so  that 
the  transform  envelope  is  also  Gaussian.  Figure  3  is  the  uncompensated  signal 
log  /'(-')  from  an  actual  electron  beam  stored  in  AGO.  The  spectrum  was  taken 
during  a  typical  synchrotron  radiation  run  with  moderately  high  current.  The 
electron  energy  was  536 MeY  and  the  RF  accelerating  voltage  was  17.5 kV.  After 
compensation  for  the  photodiode  response,  the  fit  to  a  Gaussian  envelope  gives  an 
electron  bunch-length  of  ar  —  665  psec  ±1.  1*7.  The  mean-squared  correlation  of 
the  lit  R-  >  0.002  which  gives  a  confidence  factor  in  the  Gaussian  lit  very  near 
100*7.  At  lower  currents,  such  as  those  used  for  zero-current  bunch-lengthening 
measurements  at  210 MeV.  the  typical  bunch  length  measured  is  ~200psec  with 
an  uncertainty  of  -~2*7'  and  a  mean-squared  correlation  R2  >  0.98.  This  is  in  good 
agreement  with  the  theoretical  previsions  for  AGO  [25]. 

For  FFL-induced  bunch  lengthening  the  fractional  bunch-length  change  is  the 
principal  quantity  of  interest.  The  method  of  measurement  is  simplier  for  induced 
changes  than  that  used  for  measuring  absolute  bunch  lengths.  If  the  induced 
bunch-length  changes  are  relatively  small,  or  maintain  their  functional  form,  scaling 
relations  from  Fourier  tranform  analysis  require  that  the  product  of  the  second 
moments  remain  constant  because  the  total  area  under  the  pulse  remains  constant 
[26],  For  the  storage  ring  this  is  applicable  since  the  total  charge  of  the  electron 
bunch  is  conserved.  This  invariance  remains  valid  even  under  more  general  pulse 
modifications  as  a  special  case  of  Parseval’s  theorem,  and  is  related  to  the  bandwidth 
theorem  and  Schwartz’s  inequality  [26],  If  a  single  harmonic  of  the  comb  spectrum 
of  /'p)  is  examined  the  change  in  the  amplitude  will  be  inversely  proportional  to 
the  bunch-length  change.  On  AGO,  where  the  direct  measurements  of  the  bunch 
length  have  shown  that  the  electron  distribution  is  essentially  Gaussian,  a  change 


from  alt  to  crtf  results  in  a  change  of  F{oj)  of 


,  (r,w\  -_'2a^= 

lnl7^J=  — 2 — ■ 


where  Arr:  =  a j  —  of.  The  change  in  the  signal  amplitude,  induced  by  a  change  in 
the  bunch  length,  scales  as  the  square  of  the  frequency  at  which  the  measurement 
is  made.  It  is  important  to  note  that  measurement  of  the  fractional  bunch-length 
change  is  independent  of  the  response  of  the  detector  system. 

Figure  -4  is  a  plot  of  the  detected  bunch-lengthening  signal  from  ACO  as  a 
function  of  the  detection  harmonic  during  a  single  bunch-lengthening  measurement. 
The  least-squares-fit  parabola  shows  clearly  that  for  harmonics  as  high  as  100vJo  the 
low  current  bunch  on  ACO  remains  Gaussian.  This  justifies  the  use  of  the  Gaussian 
model  for  the  electron  bunch.  The  error  bars  are  calculated  directly  from  the 
fluctuations  on  the  measurement  and  represent  ±lcr.  Since  the  fluctuations  of  the 
electron  bunch  length  are  the  source  of  the  uncertainty  the  signal-to-noise  should 
remain  essentially  constant  as  long  as  the  noise  levels  of  the  detection  system  are  not 
significant.  Measurements  at  higher  harmonics  are  dominated  by  the  decreasing, 
continuous-carrier  harmonic  level.  The  noise  level  of  the  detection  system  begins 
to  dominate  at  extremely  high  harmonics  for  which  the  carrier  amplitude  is  too 
small.  This  effectively  imposes  a  limit  assuring  that  measurements  on  ACO  will  be 
performed  in  the  frequency  regime  where  the  electron  RF  spectrum  is  Gaussian. 
The  data  of  Figure  4  were  taken  using  synchronous  detection. 

Combining  Eqns.  (2),  (3),  and  (6)  the  induced  fractional  bunch-length  change 
Arr;  may  be  directly  related  to  the  theoretical  predictions  for  free-electron-laser 
bunch  lengthening: 


This  permits  direct  verification  of  FEL  bunch-heating  models  without  recourse  to 
the  absolute  bunch  length  which  must  be  compensated  for  the  dectector  response. 
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Tin*  sensitivity  and  precision  achievable  by  the  frequency-domain  measure¬ 
ments  is  superior  to  that  achievable  with  the  available  temporal  methods.  The 
detection  system  response  is  largely  eliminated,  or  is  easily  compensated,  for  even 
non-Gaussian-response  systems.  Rapid  oscillations  of  the  electron  bunch  results 
only  in  sidebands  at  the  carrier  harmonics  which  can  easily  be  eliminated.  The 
frequency-domain  technique  is  applicable  to  situations  where  disimilar  multiple 
bunches  are  present  in  the  storage  ring  [7],  Additionally,  the  bunch  lengthening  can 
be  monitored  eontinously  for  time-dependent  response  [6],  While  streak  cameras 
have  comparable  resolution,  precise  absolute  measurements  require  multiple  correc¬ 
tions  for  electron  transit-time  spread  and  film  density.  No  such  corrections  are 
required  in  our  frequency  domain  measurements  in  which  transit-time  spread  in 
the  photodiode  has  no  effect  on  the  detected  signal. 


Basic  Experimental  Results 

The  measurements  presented  in  this  paper  were  made  with  the  NOEL  permanent- 
magnet  undulator  and  optical  klystron  mounted  on  AGO  [27].  They  were  also  made 
at  total  average  currents  of  less  than  1  inA.  This  current  level  is  dictated  by  the 
thresholds  of  anomalous  bunch  lengthening  on  AC'O  [2,1].  The  FEL  oscillator  could 
never  aeheive  ignition  if  operation  was  ever  attempted  at  this  current  level  so  all  of 
the  bunch-lengthening  data  were  obtained  while  operating  the  SRFEL  as  an  amplier 
to  an  external  laser  [8,0].  All  of  the  measurements  presented  here  were  obtained 
with  synchronous  detection  with  the  exception  of  the  time-resolved  response  of  the 
electron  bunch  (Figure  '>). 

Table  I  outlines  the  experimental  parameters  present  at  the  time  the  time- 
resolved  data,  Figure  5,  were  aquired.  These  parameters  are  typical  for  all  of 
the  basic  results  presented.  Any  differences  will  be  noted  with  the  individual 


measurement. 


Figure  5  is  the  time- averaged  response  of  the  induced  quadratic  bunch-length 
change  as  a  function  of  time.  In  the  small-signal  stochastic-heating  model  of  SRFEL 
the  induced  energy  spread  has  the  characteristic  scaling  given  by  the  ratio  of  the 
energy  damping  to  the  orbit  period  of  the  ring  as  can  be  seen  in  Eqn.  (7).  Note 
that  a  decrease  in  signal  logFfu;)  is  an  increase  in  bunch  length.  The  measured 
characteristic  damping  time  is  re/2  =  30  db  3  msec.  The  error  on  the  indivudual 
measurements  of  the  damping  time  were  actually  much  smaller  than  this  and  the 
majority  of  the  uncertainty  arises  from  shifts  in  the  mean  energy’  of  ACO  from 
run  to  run.  The  damping  time  was  obtained  by  making  a  nonlinear  fit  of  the 
time-averaged  response  from  several  runs  to  the  equation: 

/(0  =  -di  exp(M2)  +  A0.  (8) 

The  fits  were  made  using  a  \2  reduction  search.  The  confidence  ratio  of  the  fit 
is  >  OS'7.  The  theoretical  damping  time  for  ACO  at  2-10 MeV  is  re  =  05  msec. 
Thus  the  measured  characteristic  time  is  in  excellent  agreement  with  theory.  The 
induced  bunch  lengthening  on  Figure  5  corresponds  to  an  increase  in  bunch  length 
of  --•■>'7.  The  data  were  obtained  by  using  a  digital  oscilloscope  and  averaging  over 
0  1  scans.  The  magnitude  of  the  bunch  lengthening  and  the  exponential  response 
are  consistent  with  existing  bunch-lengthening  models  [1,13,15.20], 

Figure  0  plots  the  detected  quadratic-bunch-length  change  as  a  function  of 
the  RF-accelerating  voltage  for  a  constant  external-laser  intensity.  The  constant 
of  proportionality  between  the  scaling  of  the  energy  spread  and  the  bunch  length 
is  given  by  the  ratio  of  the  momentum-compaction  factor  and  the  synchrotron 
frequency  as  seen  in  Eqn.  (2).  Additionally, 

J?ocf_1.  (9) 

-1/2 

The  nominal  bunch  length  is  proportional  to  \/Y  ,  where  V  is  the  peak  RF- 

accelerating  voltage  of  the  cavity  [10],  The  detected  signal  Act?  should  also  obey 
the  l/Y  scaling  as  long  as  the  FEL-induced  effects  on  the  beam  do  not  alter  the 


synchrotron  frequency,  or  the  coupling  between  the  energy  spread  and  the  bunch 
length.  The  solid  line  in  Figure  6  is  the  least-squares  fit  to  Act2  —  bV .  The  errors 
for  each  of  the  data  points  is  smaller  than  the  symbols  representing  them  and  so 
are  not  included.  The  proportionality  of  Eqn.  (9)  should  hold  as  long  as  the  RF 
over-voltage  is  large  compared  to  the  energy  loss  per  turn,  ( U0  =  280 eV  on  AC'O 
at  2IOMeY).  At  RF  voltages  comparable  to  the  energy  loss  per  turn  the  potential 
well,  which  determines  the  synchrotron  frequency,  contains  higher  order  terms. 
These  higher  order  terms  become  comparable  to  the  quadratic  terms  resulting  in  a 
departure  from  the  1/1'  scaling  as  V  — *•  UQ- 

This  scaling  verification  is  important  since  numerical  simulations  indicate  that 
at  high  laser  intensity  a  reduction  of  the  synchrotron  frequency  should  occur  [21]. 
This  will  change  the  coupling  between  the  bunch  length  and  the  energy  spread. 
Because  of  the  possible  coupling  modifications  it  was  necessary  to  verify  that  the 
measurements  maintained  the  proper  functional  dependence  of  scaling  for  quanta- 
tive  comparision.  Unfortunately,  the  argon-ion  laser  and  optical-transport  system 
on  AC'O  are  incapable  of  providing  sufficient  optical  intensities  to  permit  the  ex¬ 
amination  of  possible  synchrotron-frequency  depression.  Because  of  this,  and  as 
Figure  (j  indicates,  all  measurements  remained  in  a  region  of  intensity  where  the 
coupling  between  the  bunch  length  and  energy  spread  remains  constant. 

Figure  7  is  a  plot  of  bunch  lengthening  as  a  function  of  the  transverse  alignment 
of  the  external  laser  with  respect  to  the  electron  beam.  Both  the  laser  beam  and 
the  electron  beam  have  finite  transverse  dimensions  so  that  the  laser  intensity  seen 
by  an  electron  will  not  be  constant  but  will  vary  from  pass  to  pass.  The  finite 
transverse  beam  size  of  the  laser  reduces  the  average  intensity  by  a  factor  which 
originates  from  the  integral  over  the  transverse  distribution  of  the  electrons  and 
transverse-laser  mode.  If  the  beams  are  coaxial  the  filling  factor,  7X ,  is 


for  transverse  electron  beam  dimensions  and  cry,  and  a  laser  beam  waist  w.  If  the 


beams  art'  eolinear.  that  is  parallel  but  having  a  transverse  displacement  between 
their  maxima,  the  filling  factor  has  the  form 


7  —  Jxh- 


(11) 


where 


7i  =  exp 


+ 


Vo 


(12) 


2(  w2/i  +  )  2(;t-2/4  +  er2) 

The  transverse  displacement  being  x0  and  yQ.  The  FEL-induced  bunch  lengthening 
should  have  a  Gaussian  dependence  as  a  function  of  transverse  alignment  of  the 
electron  and  optical  beams  if  both  have  Gaussian  profiles. 

Figure  7  verifies  the  colinear  filling  factor  7;.  The  solid  curve  is  a  weighted 
least-squares-fit  Gaussian  to  the  data.  The  errors  bars  are  the  calculated  statistical 
uncertainty  of  the  measurements  and  represent  ±lcr.  Its  cr  =  ,75mm  which  is  in 
good  agreement  with  the  predicted  value  of  (  0  mm)2  +  (.4  mm)2  =  (.72  mm)2,  where 
w/ 2  was  measured  at  .6mm  for  the  laser  and  crx  of  the  electrons  was  measured  at 
.4  mm  for  the  electron  beam. 

Madey  has  demonstrated,  as  part  of  the  Gain-Spread  Theorem,  that  the  FEL- 
induced  energy-spread  is  proportional  to  the  spontaneous  emission  power  of  the 
undulator  [28]: 


((^/et)2)v  = 


d2I 


For  a  constant-period  undulator  this  results  in 
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and  for  the  symmetric  transverse  optical  klystron  [16,20]: 

[{(^/e/)2)o]0pK.  =  2(1  +cos//)[{(5c/ei(^\,))e]u„d- 


(15) 


In  these  equations  //  is  the  peak  intensity  of  the  optical  field,  and  re  the  classical 
electron  radius.  The  total  number  of  periods  of  the  undulator,  or  optical  klystron, 
is  .V.  The  harmonic  of  the  on-axis  radiation  is  n  (n  is  odd);  7  is  the  filling- 
factor,  or  electron  beam  —  optical  mode  overlap  function.  The  undulator  detuning 
parameter,  //  =  2 ~n.\Str,  and  /)  =  t~(  ^  ,Y  +  .Yj),  is  the  optical  klystron  detuning 
parameter  with  .Y ,/  a  characteristic  of  the  dispersive  section  [16,29],  The  period  of 
the  undulator  is  \0.  K  =  eB0\0/2~mc 2  is  a  dimensionless  measure  of  the  magnetic 
field  and  is  referred  to  as  the  wiggler  K  parameter,  and  £  =  nI\2/( 4  +  2A’2).  The 
resonant  energy  difference  is  &(r  =  (q  —  qr)/qr,  where 


-  go +  !*•>)■ 


(16) 


is  the  resonant  energy  of  the  linear  undulator. 

Figure  8  is  the  verification  of  the  Madey  Gain-Spread  Theorem  for  the  un¬ 
dulator  as  a  function  of  magnetic-pole-face  gap.  Figure  9  is  the  equivalent  verification 
for  the  optical  klystron.  Figures  8(a)  and  9(a)  are  the  quadratic  laser-induced  bunch- 
lengt hening  curves.  Figures  8(b)  and  9(b)  are  the  spontaneous  emission  curves  taken 
in  identical  conditions  at  the  laser-beam  wavelength  of  ol  loA.  The  excellent  agree¬ 
ment  with  the  Madey  Theorem  is  immediately  evident.  The  vertical  scales  are  the 
quadrat  ic-buncli-length  change  measured.  They  correspond  to  a  6a/a  «  GGd.  The 
horizontal  scales  represent  the  magnetic-pole-face  gap  of  the  undulator,  which  can 
be  ecpiated  to  the  magnetic  field  strength  [27],  and  the  wiggler  I\  parameter.  The 
resonant  energy  may  be  obtained  from  Eqn.  (16).  The  central  axis  corresponds  to 
magnetic-pole-face  gap  of  36.7  mm,  K  =  1.93,  and  a  resonant  energy  of  238MeV 
for  the  undulator.  For  the  optical  klystron  the  central  axis  represents  a  gap  of 
36.8mm,  l\  =  1.92.  and  a  resonant  energy  of  237 MeY. 

If  the  difference  in  the  peak  laser  intensity  and  transverse  filling  factor  is  taken 
into  account  for  the  two  curves  the  ratio  of  the  bunch  lengthening  observed  on  the 
undulator  and  the  optical  klystron  is  equal  to  the  square  of  the  total  number  of 
period-  in  each  device: 
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This  is  in  agreement  with  the  theoretical  prediction  based  on  Eqns.  (14)  and  (15). 

The  small-signal  stochastic-heating  model  of  the  SRFEL  predicts  that  for  low 
pump-laser  intensities 


Act;  oc  ((6efei)2)<,  oc  /;.  (18) 

On  ACO  the  verification  of  the  linear  response  of  the  bunch  lengthening  to  the 
pump  laser  was  done  as  follows:  Since  the  NOEL-undulator  magnetic  field  is  linear 
polarized  in  the  vertical,  y,  direction  the  electrons  oscillate  in  the  horizontal  plane. 
The  laser  is  normally  linearly  polarized  in  the  veritcal  direction.  It  is  rotated  in  the 
horizontal  plane  by  using  a  half-wave  plate.  The  initial  angle  of  the  half-wave  plate 
is  optimized  by  the  bunch-lengthening  signal,  and  bunch-lengthening  measurements 
are  then  made  as  a  function  of  the  half-wave  plate  orientation.  Angles  are  selected 
for  measurement  such  that  the  points  are  evenly  distributed  in  a  linear  space. 

Figure  10  plots  the  bunch  lengthening  as  a  function  of  the  laser  intensity.  The 
deviations  from  the  straight  line  appear  to  be  systematic  errors  since  they  were 
present  in  identical  form  over  different  experimental  runs.  They  apparently  were 
not  induced  by  either  a  transvese  displacement  of  the  laser  beam  as  the  half-wave 
plate  is  rotated,  or  the  result  of  an  elliptical  polariztion  of  the  laser  beam,  causing 
an  additional  component  with  the  proper  polarization  to  be  continually  present. 
The  elliptical-polarization  error  would  cause  all  intermediate  points  to  lie  above 
a  straight  line  whereas  the  majority  of  the  deviations  lie  below.  An  additional 
possibility  of  a  transverse-mode  deformation  was  eliminated  as  measurements  of 
the  laser  transverse  dimensions  as  a  function  of  polarization  rotation  indicate  little 
or  no  deformation.  One  possibility  remaining  is  stress-induced  birefringence  of  the 
vacuum  window,  resulting  in  a  polarization  rotation  of  the  laser  beam  as  it  enters 
the  undulator  vacuum  chamber. 


The  horizontal  uncertainty  bars  result  from  the  inaccuracies  in  the  knowledge 
of  the  power  transmission  of  the  optical  transport  system.  The  solid  line  is  the 
least-squares  fit  of  the  data  to  a  straight  line  with  x  and  y  intercepts  at  the  origin. 
Its  slope  is  1.8  ±  .1  X  10-23  sec2-cm2/W. 

The  theoretical  values  are  also  plotted  on  Figure  10.  The  small-signal  stochastic- 
heating  model  predicts  a  slope  of  2.07±.05  X  10-23  sec2-cm2/\V  under  the  operating 
conditions  of  the  experiment  and  is  represented  on  Figure  10  by  the  dashed  line. 
The  uncertainty  in  the  theoretical  value  is  the  cumulative  result  of  the  inaccuracies 
of  measurement  in  the  filling  factor  in  the  undulator,  and  the  other  storage  ring 
parameters.  The  solution  of  the  Fokker-Planck  Equation  at  the  power  levels  of 
these  data  is  indistinguishable  from  the  small-signal  stochastic-heating  model  [13]. 
The  'X’  symbols  represent  heating  values  predicted  by  numerical  simulation  of  the 
AC'O-SRFEL  amplifier  [21].  The  basic  experimental  results  show  good  agreement 
with  FEL  theory  in  both  functional  scaling  and  absolute  magnitude,  demonstrating 
that  the  basic  interaction  between  the  electrons  and  the  laser  is  well  characterized. 

Conclusion 

The  experimental  results  suggest  several  important  conclusions  concerning  the 
beam-heating  aspects  of  the  SRFEL.  The  SRFEL  is  a  stable  system.  The  induced 
heating,  resulting  from  the  FEL  on  the  storage  ring,  can  be  compensated  for  by 
the  incoherent  synchrotron-radiation  damping.  At  low  current,  where  interactions 
between  the  electrons  and  the  envirnoment  are  small  and  no  collective  effects  are 
present,  the  experimental  results  of  the  FEL  heating  are  completely  consistent  with 
basic  theory  of  the  SRFEL.  The  argeement  with  theory  is  not  only  in  the  functional 
dependence  of  the  scaling  relations,  but  in  the  absolute  magnitudes  of  the  heating 
predicted  by  theory  and  numerical  simulations.  In  addition,  on  ACO,  the  data  show 
that  the  small-signal  heating  model  gives  results  that  are  adequate  for  explaining 
behavior  and  predicting  the  magnitude  of  the  laser-induced  electron-beam  heating. 
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The  differences  of  the  small-signal  heating  model  and  and  solution  of  the  Fokker- 
Pla  nek  equation  and  numerical  simulations  [13,15,21]  will  only  become  important 
at  optical  field  densities  much  higher  than  those  accessible  on  the  ACO-SRFEL 
system. 

With  the  energy  aperature  of  the  storage  ring  dictating  the  output  power  limit 
and  efficiency  of  the  SRFEL,  improvements  must  be  made  in  the  energy  acceptance 
of  storage  rings.  For  SRFEL  systems  with  linear  wigglers  [13] 

Plater  <JtmxIPeync-  (19) 

Tlie  output  power  of  the  laser  is  P;a,er,  the  total  synchrotron-radiation  power  is 
Patjnc,  and  the  maximum  allowable  energy  spread  of  the  storage  ring  is  <r£m„.  In 
order  to  avoid  catastrophic  reduction  in  the  the  electron  beam  lifetime,  the  half¬ 
width  energy  aperature  of  the  storage  ring  Ac  should  be  at  least  six  times  the  energy 
spread  of  the  electron  beam,  At  >  6 <r£moI  [10].  This  limit  includes  the  energy- 
spread  contribution  resulting  from  FEL-indueed  beam  heating.  If  the  wiggler  energy 
envelope  is  non-C.aussian,  as  is  the  case  for  a  linear  wiggler,  the  total  energy  spread 
can  approach  ^Ac  [13,30],  Storage-ring  technology  is  such  that  Acmox  «=.  8%  as  in 
the  case  of  AD  ONE  for  the  LEDA  project  [31].  To  give  an  idea  the  limit  imposed  by 
the  aperature  on  the  output  power,  consider  the  ACO-SRFEL  [5,6]  and  the  proposed 
Stanford-SRFEL  [32].  The  maximum  allowable  energy  spread  for  the  ACO-SRFEL 
is  ~lrc,  and  for  the  Stanford-SRFEL  design  it  will  be  as  large  as  ~4  — ►  GEd. 
The  Stanford-SRFEL  has  a  proposed  total  average  current  of  1000  mA,  and  for 
this  example  consider  ACO  to  have  a  total  average  current  of  100  mA.  Assume 
that  both  SRFEL’s  will  be  storage-ring  acceptance  limited.  The  maximum  average 
power  of  the  laser  on  ACO  at  lGGMeV  (240 MeV)  over  such  conditions  would  be 
~100mW  (~300mW);  and  for  the  Stanford-SRFEL,  which  will  operate  at  lGeV, 
the  corresponding  average  power  limit  is  ~500  — ►  700  W. 

For  short-wavelength  applications  in  the  ultraviolet,  and  x-ray  regions  of  the 
electro-magnetic  spectrum  the  limit  imposed  by  the  energy  aperature  will  not  be  a 
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major  concern,  and  present  storage-ring  design  technology  should  be  sufficient  for 
sucessful  operation  [32].  Even  for  visable  wavelengths  the  energy  aperature  limit 
still  allows  the  SRFEL  to  be  a  very  attractive  high-power  source  of  tunable  coherent 
radiation.  However,  for  high-power  operation,  or  increased  high  efficiencies  for  the 
laser  the  energy  acceptance  must  be  increased.  The  energ}-  aperature  limit  will 
remain  a  central  concern  in  the  full  development  of  high-efficiencv  SRFEL’s.  Either 
storage  ring  technolog}-  must  succeed  in  improving  energ}-  acceptances  or  alternative 
systems,  such  as  the  isochronous  storage  ring  FEL  [17],  must  be  developed. 
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Appendix  -  Storage  Ring  Parameter  Measurement 

In  this  appendix  we  list  the  measurement  techniques  used  for  various  storage 
ring  parameters  necessary  for  comparing  the  laser-induced  bunch-lengthening  data 
with  theory.  In  order  to  obtain  a  quantitative  measurement  of  the  laser-induced 
bunch  lengthening  several  control  measurements  are  neccessary.  In  addition  to 
the  characterization  of  the  external  laser,  the  mean  energy  of  the  storage  ring,  the 
relative  energy  spread,  the  transverse  dimensions  of  the  electron  beam,  the  betatron 
tunes,  the  synchrotron  frequency,  the  momentum-compaction  factor,  and  average 
stored  current  must  be  measured.  Since  the  magnetic  optics  and  the  RF  cavity 


remain  constant  on  the  storage  ring,  it  is  necessary  only  to  measure  these  quantities 
once,  then  monitor  the  peak  RF-cavity  voltage  and  the  stored  current. 

The  mean  energy  is  normally  obtained  by  measuring  the  current  being  drawn 
by  the  bending  magnets  in  the  storage  ring  and  calculating  the  corresponding 
magnet  field,  and  deducing  the  energy  of  the  ring  with  a  knowledge  of  the  RF- 
cavity  frequency  [10].  On  ACO  this  method  of  measurement  is  subject  to  errors  of 
several  percent.  This  is  because  of  the  temperature-induced  voltage  variations  in 
the  shunts  used  to  monitor  the  current  supplied  to  the  bending  magnets,  and  from 
non-standard  cycling  of  the  magnets.  The  most  accurate  method  of  calibrating 
the  mean  energy,  well  suited  for  permanent-magnet  SRFEL’s,  is  that  of  measuring 
the  center  wavelength  of  the  spontaneous  emission.  The  central  wavelength  of  the 
spontaneous  emission,  either  of  the  undulator  or  the  optical  klystron,  and  the  known 
magnetic  field  given  [27]  as  a  function  of  the  magnetic-pole-face  gap  determine  the 
mean  energy  of  the  storage  ring  to  great  accuracy.  For  NOEL  on  ACO  as  a  function 
of  pole-face  gap  the  on-axis  magnetic  field  is  a  simply  an  exponential  [27].  The 
resonant  energy  of  the  undulator,  which  is  also  the  mean  energy,  is  given  by  Eqn. 
(15). 

The  two  principal  techniques  normally  employed  for  measuring  the  energy 
spread  on  storage  rings  are  not  readily  employable  on  ACO.  The  method  of  measur¬ 
ing  the  energy  spread  of  the  core  of  the  electron  distribution  by  imaging  [10,33]  is 
not  realizable  since  at  no  point  on  ACO,  where  the  beam  is  imaged,  is  the  ring 
sufficiently  dispersive  to  make  a  good  quantitative  measurment.  Scrapers,  usually 
used  for  measuring  the  tails  of  the  distribution  [33],  were  removed  from  ACO.  A 
third  method  attempted  on  ACO  for  measuring  energy  spread  was  that  of  recording 
the  betatron  tune  envelopes  when  modulated  by  the  synchrotron  oscillations  [31]. 
The  principal  problem  on  ACO  is  that  the  horizontal  tune  is  the  most  sensitive 
to  these  fluctuations,  but  there  arc  no  electrodes  on  the  storage  ring  capable  of 
effective  horizontal  excitation  or  detection. 

The  method  which  proved  to  be  more  successful  and  precise  monitored  the 


depth  of  modulation  of  the  spontaneous  emission  of  the  optical  klystron  in  or¬ 
der  to  measure  the  energy'  spread.  This  technique  permits  the  energy  spread  to 
be  measured  simultaneously  with  oscillator  operation  by  passing  the  spontaneous 
emission  through  a  monochromater  and  situating  a  photomultiplier  at  one  of  the 
minima  of  the  spontaneous-emission  curve.  Changes  in  the  energy  spread  result  in 
a  decrease  in  the  modulation  depth  and  a  subsequent  increase  in  the  detected  signal 
(16.20.35]. 

Figure  11  is  a  plot  of  the  energy  spread  detected  by  this  method  while  scanning 
RF  voltage  in  a  current  regime  where  turbulent  bunch-lengthening  thresholds  are 
present.  The  thresholds  at  ~17kV,  14 kV,  and  at  ~12kV  are  clearly  visible. 

The  transverse,  electron-beam  dimensions  were  measured  by  two  techniques. 
The  most  accurate  and  reproducible  method  was  that  of  imaging  the  synchrotron 
radiation  from  a  bending  magnet  with  a  convergent  lens  onto  a  Reticon  CCD 
[Charge  Coupled  Device],  linear  diode  array.  The  signal  from  the  CCD  was  then 
read  on  an  oscilloscope  at  a  known  rate.  The  use  of  a  digital  oscilloscope  with  a  GPIB 
interface  allowed  the  transfer  of  the  data  for  fitting.  Alignment  was  accomplished 
by  using  an  external  signal  generator  and  exciting  one  of  the  betatron  tunes.  This 
would  cause  the  beam  to  enlarge  in  one  of  its  transverse  dimensions.  Then,  the 
CCD  would  be  rotated  until  the  widest  pulse  width  was  located.  The  other  tune 
would  then  be  excited  establishing  the  other  dimension.  This  assured  that  the  two 
orthogonal  axes  of  the  electron  beam  were  being  measured. 

Betatron  tunes  were  measured  principally  by  external  excitation  and  by  ob¬ 
serving  the  transverse  beam  profile  on  one  of  the  video  cameras  which  imaged  the 
beam  in  the  storage  ring.  The  frequency  of  excitation  then  was  input  to  a  mixer 
with  the  orbit  frequency  which  gave  an  output  signal  at  the  frequency  of  that  par¬ 
ticular  tune.  This  method  suffered  from  two  shortcomings.  The  excitation  required 
for  the  two  separate  tunes  was  disproportionate  in  amplitude.  The  disproportionate 
excitation  again  was  a  result  of  there  being  no  electrode  plates  on  the  storage  ring 
capable  of  exciting  the  horizontal  tune.  Because  of  this  it  was  easy  to  lose  all  of 


the  stored  beam  by  accidently  exciting  the  strongly  coupled  tune  while  seeking  the 
weaker  tune.  Additionally,  this  method  could  be  time  consuming  if  an  accurate 
measurement  of  the  tunes  was  needed,  since  the  resonances  varied  in  their  width 
and  definition  from  run  to  run. 

To  avoid  the  problem  of  the  external-excitation  detection  of  the  betatron  tunes, 
a  technique  using  a  broadband  noise  generator  and  the  RF  spectrum  analyzer 
connected  to  an  electrode  on  the  storage  ring  was  attempted.  However,  the  tunes 
themselves  were  only  about  ~5dB  above  the  noise  level,  and  were  ~70dB  down 
in  amplitude  from  the  orbit  frequency  signal.  The  value  of  this  method  was  also 
reduced  because  of  the  unsuitable  electrodes  on  ACO.  However,  it  is  the  most 
typical  method  of  measurement  of  the  betatron  tunes  employed  on  the  principal 
storage  rings  in  operation. 

The  spectrum  analyzer  was  employed  for  the  measurement  of  the  synchrotron 
frequency.  It  is  by  far  the  most  effective  and  direct  method  of  measurement.  The 
synchrotron  frequency  could  be  measured  by  either  modulating  the  RF-accelerating 
cavity  voltage  with  an  external  low-frequency  signal  generator,  1— ►20kHz,  and 
detecting  the  level  of  coupling  between  the  excitation  and  the  electron  beam.  Or, 
a  second  technique,  which  is  non-invasive,  consists  of  recording  one  of  the  orbit 
harmonic's  natural  sideband  spectrum.  The  peaks  are  subsequently  fit  to  determine 
the  synchrotron  frequency. 

The  momentum-compaction  factor  is  the  constant  of  proportionality  between 
the  orbit  path-length  differences  and  the  energy  of  the  electrons  [10].  It  is  most  easily 
measured  by  changing  the  frequency  of  the  RF  accelerating  cavity  and  recording 
the  shift  in  the  mean  energy  of  the  stored  electrons.  The  most  sensitive  and  direct 
method  manner  of  making  this  measurement  is  to  use  the  spontaneous  emission  of 
the  optical  klystron. 

Figure  12  is  a  composite  set  of  spontaneous  emission  spectra  of  the  optical 
klystron  taken  as  a  function  of  magnetic-pole-face  gap  at  various  RF-cavity  fre¬ 
quencies.  The  resonant  energies  can  be  obtained  from  Eqn.  (16).  The  shifts  in  the 
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maxima  and  minima  were  measured  and  used  to  determine  the  induced  change  in 
energy  of  the  electrons  resulting  from  the  change  in  RF  frequency.  The  changes 
in  energy  were  in  turn  least-spuares  fit  to  a  straight  line  in  order  to  obtain  the 
momentum-compaction  factor.  The  measured  value  for  a  on  ACO  was  0.028  at 
2  10  MeY. 

The  technique  used  for  measuring  the  stored  current  on  ACO  was  to  monitor 
the  average  signal  current  from  a  phot  on.  alt  iplier  illuminated  by  synchrotron  radia¬ 
tion.  Additionally,  the  signal  amplitude  of  an  electrode  on  an  oscilloscope  was  used. 
Both  of  these  techniques  were  calibrated  with  a  DC  transformer  developed  by  M. 
Bergher  [36],  Both  allowed  very  good  relative  current  measurements,  but  limitations 
on  both  techniques  reduced  the  absolute  precision  of  stored-current  measurements 
An  impedence  mismatch  caused  the  signal  on  the  oscilloscope  to  be  inconsistent 
from  one  vertical  scale  to  the  next.  The  photomultiplier  tube  was  sensitive  to  stray 
magnetic  fields  and  changes  in  energy  of  the  storage  ring.  The  signal  from  the 
photomultiplier  tube  varied  from  run  to  run  because  of  changes  in  cycling  of  the 
ring  magnets  and  operating  point. 
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TABLE  I 


Basic  Experimental  Parameters 


Electron  Energy 

E 

238  MeV 

Momentum  Compaction 

a 

0.028 

Synchrotron  Frequency 

ft/2ir 

6190  Hz 

Electron  Bunch  Length 

<7 

200  psec 

Energy  Damping  Time 

T 

65  msec 

Electron  Beam  Width  (Horiz.) 

<?x 

0.3 1  mm 

Electron  Beam  Width  (Vert.) 

0.33  mm 

Undulator  K  Value 

K 

1.93 

Undulator  Period 

K 

7.78  cm 

Laser  Wavelength 

X 

51 15  A 

Laser  Peak  Intensity 

h 

175  W/em2 

Laser  Beam  Waist 

w 

1.2  mm 

Orbit  Frequency 

u0/2tt 

13.618MHz 

Detection  Frequency 

oj/2k 

1062.2  MHz* 

Number  of  Bunches 

1 

*78^  Harmonic 
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Figure  2  Schematic  diagram  of  the  iaser-induced  bunch-lengthening 
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Figure  3  Spectrum  Analyzer  output  showing  the  comb  spectrum  of  the 

photodiode  signal  from  synchrotron  radiation  from  ACO.  The  horizon¬ 
tal  scale  is  lOOMHz/div;  the  vertical  is  dB5/div.  The  first  line  is 
the  zero  frequency  reference  line  and  subsequent  lines  are  spaced  at 
13.018  MHz  -  the  orbit  frequency  of  the  storage  ring.  The  spectral 
envelope  corresponds  to  a  compensated  temporal  bunch  length  of 
663  picosecond  with  a  mean-squared  correlation  of  99. 2°c.  The 
spectrum  was  taken  at  an  electron  energy  of  536 MeV,  a  current  of 
90  mA,  and  an  RF  accelerating  voltage  of  17.5  kV. 
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Figure  4  Plot  of  log/qfcj)  —  logF/(u,’)  as  a  function  of  detection  frequency 


measured  as  harmonics  of  the  orbit  frequency  for  a  single  bunch¬ 
lengthening  measurement.  The  solid  line  is  a  least-squares  fit  of 
the  data  to  a  parabola,  the  response  expected  for  a  Gaussian  beam. 


The  error  bars  are  measurement  fluctuations  of  ±l<r. 


Figure  6  FEL-induced  quadratic  bunch-length  change  as  a  function  of  RF 
accelerating  cavity  voltage.  The  solid  line  is  the  least-squares  fit 
of  the  data  to  1  /Vrf.  The  statistical  errors  are  smaller  than  the 
symbols  representing  the  data  points  The  data  should  scale  as  1  /V'r/ 
as  long  as  the  FEL  interaction  does  not  effect  the  scaling  between 
the  energy  spread  and  bunch  length. 


FEL-induced  quadratic  bunch-length  change  as  a  function  of  trans¬ 
verse  laser  alignment.  The  solid  line  is  the  weighted  least-squares 
fit  to  a  Gaussian  with  a  a  =  0.75  mm,  in  good  agreement  with  the 
theoretical  value  of  0.72  mm  for  the  colinear  filling  factor.  The  error 
bars  are  measurement  fluctuations  of  ±1<t. 
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Figure  8  (a)  FEL-induced  bunch  lengthening  as  a  function  of  permanent- 

magnet  unduiator  magnetic-pole-face  gap.  The  maximum  of  the 
signal  corresponds  to  a  bunch  lengthening  of  at  a  gap  of 

36.7  mm.  equivalently  a  K  =  1.93,  or  a  resonant  energy  of  238  MeV. 
The  error  bars  are  measurement  fluctuations  of  ±I<x.  (b)  Spontaneous 
emission  of  the  permanent  magnet  unduiator  as  a  function  of  magnetic 
pole-face  gap  at  a  fixed  wavelength  (5145A). 
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Figure  9  (a)  FEL-induced  bunch  lengthening  as  a  function  of  permanent- 

magnet  optical  klystron  magnetic-pole-face  gap.  The  central  axis 
corresponds  to  a  gap  of  36.8  mm,  a  I\  =  1.92,  and  a  resonant 
energy  of  237 MeV.  When  difference  in  incident  laser  intensity  is 
accounted  for  the  maximum  observed  lengthening  is  a  factor  of  1.5 
less  than  for  the  undulator  in  agreement  with  theory.  The  2cr  error 
bars  shown  are  calculated  from  the  statistical  fluctuations  of  the 
signal  before  the  scan  was  taken,  (b)  Spontaneous  emission  of  the 
permanent-magnet  optical  klystron  as  a  function  of  magnetic-pole- 
face  gap  at  a  fixed  wavelength  (5145A). 
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Figure  10  Quadratic  laser-induced  bunch  lengthening  as  a  function  of  peak 
laser  intensity  obtained  by  polarization  rotation.  The  solid  line  is 
least-squares  fit  of  to  a  straight  line  through  the  origin,  and  has 
a  slope  of  1.8  ±  .1  X  10~23  sec2-cm2/W.  The  dashed  line  is  the 
value  predicted  by  the  stochastic  small-signal  heating  model  based 
on  the  operating  conditions  present  has  a  slope  of  2.07  ±  .05  X 
10-23  sec2-cm2/\V.  The  vertical  error  bars  are  are  measurement 
fluctuations  of  ±1(7.  The  horizontal  error  bars  are  the  result  of 
laser-power  transmission  uncertainties. 
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Abstract.  To  improve  the  gain  in  the  Orsay  storage  ring  Free  Electron  Laser  (FEL) 
experiment,  the  17  period  permanent  magnet  undulator  has  been  modified  to  form  an 
optical  klystron  (OK).  VVe  report  the  measurement  of  spontaneous  emission  and  the  effects 
on  it  of  energy  spread  and  angular  spread.  Gain  and  laser  induced  bunch  lengthening 
measurements  with  the  OK  are  also  reported  and  are  in  very  good  agreement  with  the  FEL 
classical  theory.  The  spontaneous  emission  spectrum  which  is  easy  to  measure  with  good 
signal  to  noise  ratio,  turns  out  to  be  a  very  good  diagnostic  tool  for  energy  spread  and  angular 
spread  measurements  on  storage  rings.  The  factor  of  four  increase  in  the  small  gain  obtained 
by  converting  the  undulator  NOEL  into  an  OK  was  the  critical  factor  in  the  recent 
operation  of  the  ACO  storage  ring  laser  above  threshold. 

PACS:  42.60,  42.55 


The  main  purpose  of  the  Orsay  experiment  is  to  prove 
the  feasability  and  test  the  theories  of  a  storage  ring  free 
electron  laser  (FEL)  in  the  visible  range.  A  permanent 
magnet  undulator  has  recently  been  built  [1,2]  and 
has  successfully  operated  on  the  electron  storage  ring 
ACO  in  the  energy  range  150-540MeV  [3],  However 
ACO  is  not  optimized  for  FEL  studies.  The  relatively 
low  electron  density  (for  storage  rings)  and  the  short 
length  of  the  available  straight  section  have  conspired 
to  limit  the  gain  available  with  an  undulator  to  1  to 
2  10“ 4  per  pass  at  a  wavelength  /.  =  6300  A  and 
240  McV  electron  energy.  Even  with  state-of-the-art 
mirrors,  laser  operation  is  impossible  at  this  level  of 
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gain.  Several  directions  have  been  followed  to  improve 
the  gain.  The  most  successful  approach  has  been  the 
modification  of  the  undulator  into  an  optical  klystron 
(OK).  In  this  paper  we  shall  discuss  our  experience  with 
the  OK.  The  OK  originally  proposed  by  Vinokurov 
and  Skrinsky  [4]  consists  of  two  identical  undulators 
separated  by  a  dispersive  section  forcing  the  electron 
into  a  single  large  wiggle  (Fig.  1).  This  configuration 
has  a  higher  gain  than  an  undulator  of  the  same  total 
length.  Such  a  device  can  be  used  to  advantage  on 
electron  beams  with  energy  spread  and  it  allows  the 
maximization  of  the  gain  in  an  interaction  region  of 
fixed  length. 

Early  results  have  already  been  reported  [3],  In  this 
paper  we  report  on  the  dispersive  section  optimization 
(Sect.  1),  spontaneous  emission  measurements  (Sect.  2), 
gain  measurements  (Sect.  3)  and  laser  induced  bunch 
lengthening  (Sect.  4).  The  spontaneous  emission  results 
are  very  detailed  because  they  arc  easy  to  measure  with 
a  good  signal  to  noise  ratio  and  give  much  information 
on  the  OK  FEL  behaviour. 

We  use  the  notation  of  the  theoretical  description  of 
optical  klystrons  of  [5]. 
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away  from  the  axis  of  the  electron  beam  (Fig.  2).  The 
calculated  trajectory  is  shown  in  Fig.  1  which  shows  a 
large  wiggle  of  about  2  mm  amplitude  for  an  electron 
energy  of  240  MeV  and  a  gap  of  33  mm. 


2.  Spontaneous  Emission 


2.1.  General  Features  of  the  Fundamental 

As  was  done  for  the  undulator,  the  electron  beam  is 
initially  aligned  by  optical  means  into  the  axis  of  the 
vacuum  chamber  within  ±  1  mm  [3]. 

The  spontaneous  emission  pattern  of  the  optical 
klystron  looks  the  same  to  the  eye  as  that  of  an 
undulator  [7,  6,  3],  namely  a  series  of  concentric 
coloured  rings  with  a  similar  pattern  produced  by  each 
harmonic  at  progressively  larger  opening  angles.  The 
big  difference  appears  in  the  spectrum. 

The  experimental  set  up  used  to  measure  the 
spontaneous  emission  spectrum  has  already  been 
described  [6],  It  consists  of  a  1750mm  focal  length 
spherical  mirror  placed  at  about  6.5  m  from  the  center 
of  the  optical  klystron. 

The  light  transmitted  by  a  75  pm  pinhole  placed  at  the 
focal  distance  from  the  mirror  is  sent  through  a  lens 
into  an  M20  uv  Jobin  Yvon  monochromator  (4  A 
resolution  with  0.1mm  slits).  The  output  of  the 
monochromator  is  then  sent  into  a  Hamamatsu  R  928 
photomultiplier.  Fig.  3  shows  a  spontaneous  emission 
spectrum  of  the  fundamental  at  an  energy  of  240  MeV 
for  a  gap  of  34.4  mm  at  low  current  in  the  ring. 

We  have  compared  the  envelope  of  the  oscillations 
with  that  of  the  emission  spectrum  of  a  perfect 
undulator  having  exactly  N  sinusoidal  periods,  the 
fields  outside  these  periods  being  exactly  zero.  The 
emission  spectrum  dl/df.  of  such  an  undulator  is: 


dl  1 

dr.  /. 


(3) 


with  <5  =  «rr/V  (1  —  /„//. ) , 

where  /.R  is  the  resonant  wavelength  and  n  is  the 
harmonic  number.  Fitting  the  envelope  to  curves  given 
by  (3)  gives  TV  =  8. 1  ±0.1  instead  of  7.  This  discrepancy 
is  probably  due  to  the  dispersive  section  field  which 
could  be  partly  resonant  with  the  other  7  periods.  It  is 
certainly  not  due  to  the  fringe  field  of  the  half  periods 
which  slightly  decreases  the  effective  number  of 
periods.  The  envelope  of  oscillations  also  presents  a 
long,  short  wavelength  tail  with  the  secondary 
maximum  amplitude  lower  than  expected.  This  effect 
was  also  seen  in  the  17  period  undulator  emission 
curves  and  is  due  to  the  parasitic  tail  of  the  ACO 
bending  magnets  fringe  field  [2],  Similar  curves  of 
emission  were  obtained  as  function  of  energy  or  field  in 
the  dispersive  section  (by  changing  the  gap). 


VERTICAL  MAGNETIC  FIELD 


Fig.  1.  Vertical  magnetic  field  calculated  for  the  Orsay  optical 
klystron  (gap:  33mm)  and  the  corresponding  calculated 
horizontal  electron  trajectory  at  an  energy  of  240  MeV 


Fig.  2.  Dispersive  section  permanent  magnet  configuration 
optimizing  the  low  field  gain 


Fig.  3.  Spontaneous  emission  spectrum  dlld/.dQ  measured  for  an 
electron  energy  of  238  MeV  and  a  magnetic  field  parameter  of 
K  =2.09  at  low  current  where  the  modulation  is  almost  total.  The 
current  decay  /(I)  is  superimposed 


Plotting  the  Nj  of  each  maximum  as  a  function  of  1/A, 
we  verified  a  high  linearity  (0.9999  correlation 
coefficient  in  the  range  0.3  <  A  <0.7  [pm]).  From  the 
slope,  we  can  calculate  the  experimental  value  of  Nd 
which  is  iVrf  =65.3+0.2  at  ymc1  =  240  MeV  and 
A  =  6238  A  for  a  dispersive  section  gap  of  35  mm.  This 
value  is  close  to  the  value  A7d  =  68  predicted  from  the 
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Fig.  4.  Spontaneous  emission  spectrum 
ill  il/.tld  measured  for  an  electron 
energy  of  240  McV.  a  current  of  40mA 
and  a  magnetic  gap  of  35  mm.  The  non 
complete  modulation  is  due  to 
fractional  energy  spread  of  1.2  x  10“ 3 


Fig.  5.  Spontaneous  emission  spectrum 
dI  J/.JQ  measured  at  low  curent  and 
240  MeV  The  modulation  is  reduced 
below  that  of  figure  3  by  inducing  large 
angular  motion  in  the  beam  with  a 
broad  band  noise  source  applied  to  an 
electrode  in  the  vacuum  chamber.  Note 
the  non-symmetric  modulation 
characteristic  of  the  angular  spread  as 
opposed  to  the  symmetric  modulation 
observed  in  the  case  of  a  dominant 
energy  spread  (Fig.  4) 


attributed  entirely  to  the  energy  spread  in  this  case  for 
two  reasons.  All  the  other  contributions  are  predicted 
to  be  negligible  as  we  will  see  in  the  next  two  sections, 
and  the  deduced  energy  spread  [5],  namely 


y 


1 

4n(  N  +  Nj) 


[/  — 21og/=  1.2  x  10  •' 


is  consistent  with  the  value  1.4  x  10  3  +  0.2  predicted 
from  the  measured  bunch  length  assuming  constant 
synchrotron  frequency. 

The  theory  also  predicts  that  for  a  Gaussian  energy 
spread  /~cxp(  —  er/2)  with  ox  1/A.  Calculating  for 
the  eleven  largest  fine  structure  peaks  one  verifies  the 
proportionality  of  a  versus  I j/.  with  a  correlation 
coefficient  of 0.95.  Such  a  test  is  not  powerful  since  1 !'/. 
does  not  change  much  inside  the  bandwidth  of  the 
fundamental. 


We  have  also  measured  /  and  deduced  a  for  the  first 
three  harmonics  recorded  at  the  same  energy,  current 
and  gap.  This  data  can  be  used  for  a  more  precise  check 
of  the  /.  dependence.  Table  4  summarizes  the 
measurements  after  deconvolution  from  the 
monochromator  response. 

In  Table  4  the  measured  <r  is  compared  to  the  predicted 
values  from  three  different  energy  distribution  shapes 
all  normalized  the  fundamental.  We  observe  a  good 
agreement  with  the  Gaussian  energy  spread  predicted 
from  storage  ring  theory  [9],  The  accuracy  of  this  test 
on  the  energy  distribution  depends  on  the  assumption 
that  for  a  monoenergetic  filament  beam  the 
modulation  rate  /  is  equal  to  1.  Residual  modulation 
can  be  produced  either  by  field  errors  in  the  OK,  or  an 
error  in  the  pinhole  positioning  with  respect  to  the 
mirror  focal  point.  We  have  measured  a  0.935 
modulation  depth  at  very  low  current,  subtracted  all 
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Fig.  6.  Effect  on  the  modulation  rate  of  the  horizontal  angle  of 
observation  with  respect  to  the  electron  trajectory  axis.  The  curve 
plotted  is  the  least  squares  fit  of  -  /.’log  (mduiation  rate)  to  a 
second  order  polynomial  of  the  angle.  The  fit  gives  the  horizontal 
angular  spread  aa  0.12 ±0.01  mrad 


where  ax  and  ay  are  the  horizontal  and  vertical  rms 
transverse  spreads.  Qx  and  Qy  are  coefficients 
connected  to  the  field  gradients  in  the  dispersive 
section.  They  only  depend  on  the  field  geometry  and 
are  equal  to  zero  if  the  field  is  exactly  uniform.  Qy  =  3.4 
x  10  4 mm  2  was  deduced  from  field  measurement; 
this  value  is  just  10%  lower  than  anticipated  from  the 
field  calculation.  These  calculations  also  predicted  Qx 
=  8.4x  10' 5  mm  2,  a  value  which  has  been  roughly 
confirmed  by  the  field  measurements.  Taking 
<Tx  =  <T,.  =  0.35mm  (low  current)  and  ox  =  oy  =  0.5  mm 
(30  mA  of  current  at  240  MeV)  on  has  /  =  0.9985  (low 
current)  to  /  =0.997  (high  current).  Such  a  tiny  effect 
was  covered  experimentally  by  the  energy  spread  and 
the  angular  spread.  In  some  cases,  verification  may  by 
possible  by  injecting  the  electron  beam  at  a  distance 
(x0,  v0)  from  the  dispersive  section  axis.  Equation  (12) 
still  applies  if  one  replaces  tT4(<r4)  by  cr4 +  2.Xq<t2 
(zt4  +  2i'(2a2).  Such  an  experiment  was  not  possible  on 
ACO  because  of  the  change  in  the  beam  focussing 
which  would  have  occurred  at  the  same  time.  The 
modulation  rate  only  depends  weakly  on  the  injection 
point  in  the  dispersive  section.  As  in  previous 
subsection  one  can  define  the  apertures:  horizontal 
aperture:  ± 20 mm;  vertical  aperture:  ±4.7mm. 

2. 3.  Harmonics 

Figures  8  and  9  show  the  second  and  third  harmonics 
recorded  at  240  MeV  and  a  gap  of  32.60  mm  at  0.3  mA 


-XJ  logit) 


Fig.  7.  Effect  on  the  measured  modulation  rate  of  the  vertical 
angle  of  observation  with  respect  to  the  electron  trajectory  axis. 
The  fit  gives  the  vertical  angular  spread  <r„  =  0.1 2  +  0.005  mrad 


and  5  mA  of  current,  respectively.  As  for  a  regular 
undulator,  an  ideal  optical  klystron  does  not  emit  any 
2nd  harmonic  on  axis  and  in  fact,  the  measured  peak 
intensity  in  Fig.  8  was  only  8%  of  that  of  the 
fundamental.  This  residual  intensity  is  due  to  the 
trajectory  effect  of  the  fringing  fields  of  the  two 
adjacent  bending  magnets  of  ACO.  The  third 
harmonic  has  a  shape  similar  to  the  first  harmonic. 
Fitting  for  N  to  the  envelope  curve  of  (3)  gives  A7  =  3.8 
instead  of  7.  We  know  from  the  3rd  harmonic 
measurements  on  the  17  periods  undulator  that 
inhomogeneous  effects  are  not  responsible  for  this 
broadening  [3],  The  dispersive  section  must  therefore 
be  responsible  for  this  effect  either  via  an  imperfect 
compensation  or  through  a  destructive  interference 
introduced  by  the  dispersive  section  at  the  undulator 
wavelength. 

2-4.  Dispersive’  Section  Spontaneous  Emission 

The  dispersive  section  is  equivalent  to  a  three  pole 
wiggler  and  therefore,  has  a  broad  emission  spectrum- 
If  one  is  just  interested  in  the  FEL  gain  improvement, 
one  can  ignore  the  dispersive  section  intrinsi  emission. 
However,  apart  from  the  OK  fringes  present  across  all 
the  spectra  wc  have  observed  (0.2  to  0.7  pm)  and  apart 
from  the  broadening  or  narrowing  of  the  emission 
curves  (which  we  have  already  discussed),  there  arc 
some  unique  features  which  appear  in  the  spontaneous 
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Fig.  1 1  a  and  b.  Transverse  horizontal  beam  profiles  of  the  light 
emitted  by  the  electrons  in  the  optical  klystron.  In  (a)  the 
interference  filter  center  wavelength  falls  inside  the  undulator 
fundamental  resonance;  the  peak  is  mainly  due  to  the  emission  by 
electrons  in  the  two  undulators.  The  left  hand  shoulder  is  mainly 
due  to  the  emission  in  the  dispersive  section.  It  is  still  present  in 
(b)  where  the  gap  was  decreased  to  suppress  completely  the 
undulator  emission 


and  not  in  the  vertical.  The  peak  of  Fig.  lib  shows 
maximum  definition  when  the  CCD  is  positioned  in 
the  image  plane  of  the  dispersive  section.  This  peak 
corresponds  to  light  emitted  from  a  source  displaced 
about  1 .5  mm  from  the  undulator  axis  as  compared  to 
the  theoretical  2  mm  off  axis  excursion  of  the  electrons 
in  the  dispersive  section. 

Optical  Klystron  Emission  at  Large  Gap.  At  large  gap 
the  undulator  field  vanishes  exponentially,  much  faster 
than  the  dispersive  section  field.  The  usual  coloured 
ring  [7,  6,  2]  due  to  the  interference  of  light  from  the 
undulator  periods  vanish.  However  one  still  sees  a 
double  interference  structure  centered  about  the 
undulator  axis.  A  black  and  white  reproduction  is 
given  in  Fig.  13  for  the  gap  series  {142, 121, 182  mm}  at 
240  MeV.  The  dispersive  section  maximum  field  is 
{660.  1150,  1570Gauss},  more  than  20  times  larger 
than  the  undulator  peak  field.  This  emission  is  not  due 
to  the  undulator  but  is  produced  by  the  dispersive 
section  with  some  contribution  from  the  storage  ring 
bending  magnets  fringe  field. 

The  emission  produced  in  the  dispersive  section  is 
contained  within  the  maximal  tangent  directions 
labeled  A  and  B  in  Fig.  12a.  At  large  distances  the 
radiation  pattern  is  symmetric  about  the  axis  0. 
Schematically,  due  to  the  small  I/y  emission  cone,  an 
observer  in  the  electron  orbit  plane  between  0  and  A 
will  see  the  emission  of  the  f,  — 12  part  of  trajectory 
shown  in  Fig.  12b.  In  the  0  direction  with  respect  to  A, 
the  emission  is  mainly  due  to  small  regions  around 
tangent  points  f,  and  t2  where  the  curvature  is 
important  (between  t ,  and  t2  the  curvature  is  very 
small  and  the  emission  of  this  part  of  trajectory  can  be 
neglected).  Both  points  sources  (,  and  f2  interfere 
constructively  when  the  difference  in  electron  and 
photon  transit  times  is  an  integral  number  of  optical 


Fig.  13.  Three  photographs  filtered  at  5500  A  of  the  interference 
structure  produced  in  the  optical  klystron  at  large  gap.  From  top 
to  bottom,  the  magnetic  field  is  increased  by  closing  the  gap.  The 
beam  axis  is  marked  by  the  black  cross  at  the  center  of  each 
pattern,  and  the  points  A  and  B  calculated  from  the  magnetic 
field  strength  in  the  OK  arc  marked  with  white  crosses 


periods: 
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where  d  is  practically  half  of  dispersive  section  length 
and  n  is  an  integer. 

For  monochromatic  light  the  interference  maxima 
occur  at  a  constant  angle  0  tracing  out  circles  around  A 
in  the  observation  plane.  Of  course,  the  circles  are 
incomplete  because  the  angle  of  emission  goes  outside 
the  electrons  l/y  forward  cone,  where  the  intensity 
drops  to  zero.  The  same  interference  pattern  occurs 
around  B  direction,  and  one  secs  a  double  colored 
interference  structure  centered  on  A  and  B  directions 
as  illustrated  in  Fig.  13,  for  /  =  5500A  and 
£  =  240McV.  Positions  in  fringe  are  in  a  rather  good 
agreement  with  theoretical  values  calculated  from  ( 1 3), 
where  =  8  to  10  cm,  the  angle  of  A  direction  being 
3.6mrad,  6mrad  and  7.9mrad,  respectively,  for  each 
photographs. 

Paticrns  like  Fig.  13  have  been  recorded  at  240, 
540  MeV,  and  several  different  gaps.  As  the  magnetic 
field  in  the  dispersive  section  is  increased  by  reducing 
the  gap,  the  A -Bangle  increases,  and  more  interference 
bands  become  visible.  At  small  gap,  the  usual 
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Nd  =  7l 


Fig.  15.  The  spontaneous  emission 
spectrum  as  a  function  of  wavelength 
and  undulator  gap  or  magnetic  field 
are  compared.  The  rapidly  oscillating 
measured  curves  fit  well  into  the 
theoretically  calculated  envelope  curves 
with  some  small  deviations  in  the  case 
of  the  wavelength  scan.  The  slower 
dependence  of  the  dispersion  parameter 
Nd  on  gap  is  visible  in  the  lower  curve. 
The  gain  as  a  function  of  gap  is 
approximately  the  derivative  of  the  gap 
scan,  not  the  wavelength  scan 


where  Gok/G17  is  the  OK  gain  improvement  with 
respect  to  the  original  17  period  undulator  for  a 
monoenergetic  filament  beam.  It  is  given  by  (4);  see  also 
Table  3. 

As  before,  /  is  the  modulation  rate  and  G17  is  the 
17  period  undulator  peak  gain  which  is  equal  to  [17] 

G,,  =  1.5x10  12^K2[JJ]2^,F/,  (15) 

where  ?.0  is  the  undulator  period  in  cm,  [JJ] 2  is  the 
Bessel  function  factor  [13].  K  =  eB/.J2nmc  is  the  peak 
field  magnetic  parameter,  qc  the  peak  electron  density 
in  cm  3,  and  F {  is  the  filling  factor  which  can  be 
calculated  from  the  overlap  integral  between  the 
electron  and  laser  beams.  For  coaxial  weakly  diverging 
beams  [14], 


where  »v0  is  the  laser  beam  waist  and  ah  and  <r„  the 
electron  beam  transverse  rms  sizes. 

The  theoretical  peak  gain  predicted  for  the 
measurement  shown  on  Fig.  14  is 

( G.,k ),„  =  4.0( G17)th  =  6.6 x  10  4±30%  (17) 

with  Gok/G, 7  =  5.9  (Table  3),  f  =  0,67  (measured  on 
the  spontaneous  emission)  and  (G,  7)lh  calculated  from 
(15,  16).  The  30%  error  bars  originate  primarily  from 


the  uncertainty  in  the  /)  functions  in  the  undulator.  The 
measured  peak  gain  is  (Fig.  14): 

<Gok)mea5  =  7  10-4(±0.3).  (18) 

The  agreement  between  the  theoretically  predicted  (17) 
and  measured  (18)  values  is  remarkably  good,  and 
certainly  the  best  achieved  to  date.  This  was  made 
possible  by  improvements  in  the  alignment  technique, 
and  by  thorough  characterisation  of  the  laser  mode. 
Not  only  does  the  values  of  he  gain  reported  here 
confirm  the  validity  of  the  theory,  but  it  also  represents 
a  large  increase  over  the  gain  measured  with  the 
undulator  NOEL. 


4.  Bunch  Lengthening 

Bunch  lengthening  experiments  have  been  performed 
by  superposing  an  argon  laser  colinear  with  the 
electron  beam  inside  the  optical  klystron.  The 
experimental  set  up  has  already  been  described  [10, 3], 
good  agreement  was  observed  with  stochastic  heating 
models  [15]  at  very  low  current  while  high  current 
regimes  are  dominated  by  anomalous  bunch 
lengthening  effects.  Figure  16  shows  the  bunch 
lengthening  and  spontaneous  emission  produced  by 
the  optical  klystron.  The  vertical  line  indicates  an 
undulator  parameter  K  of  1 .92.  The  results  are  in  good 
agreement  with  the  theorem  demonstrated  by  Madey 
that  the  mean  squared  energy  spread  is  proportional  to 
the  spontaneous  power  spectrum  [8]. 
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ABSTRACT 

We  have  simultaneously  measured  [1]  the  spectra  of  the  gain 
and  the  forward  spontaneous  radiation  in  the  Orsav  free  electron 
laser  so  that  we  can  extract  the  relative  shift  of  these  two 
curves.  Our  results  are  in  agreement  with  the  Madey  theorem  in 
the  low  divergence  limit,  and  demonstrate  the  violation  of  the 
theorem  produced  as  the  Rayleigh  range  of  the  wave  becomes 
comparable  to  or  smaller  than  the  length  of  the  interaction 
region . 
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Introduction 


The  Madey  theorem  [2]  has  proved  to  be  a  valuable  guide  to 
research  with  free  electron  lasers.  We  will  consider  the  form  of 
the  theorem  with  greatest  experimental  value:  the  gain  spectrum 
is  related  to  the  derivative  of  the  forward  spontaneous  emission 
spectrum. 


G(y) 
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The  Madey  theorem  simplifies  the  calculation  of  the  gain  and  the 
induced  energy  spread  in  an  arbitrary  magnetic  structure  [2,3,4], 
and  shows  that  the  easily  measured  spontaneous  spectrum  is  a 
reliable  diagnostic  of  the  small  signal  gain  [5],  which  is  a 
difficult  quantity  to  measure.  In  view  of  the  value  of  this 
theorem,  it  is  important  to  identify  its  limitations  and  where 
possible  confirm  them  experimentally. 


The  derivations  of  the  theorem  [2, 3, 6, 7, 8]  all  assume  low 
gain,  an  unsaturated  electron-photon  interaction,  and  a  plane 
wave  electromagnetic  field.  We  have  shown  experimentally  [9]  that 
the  theorem  is  obeyed  when  the  above  conditions  are  satisfied.  A 
number  of  authors  (see  for  example  references  [10],  and  [11]) 
have  derived  the  theoretical  result  that  the  gain  spectrum 
distorts  as  the  saturated  or  the  high  gain  regime  is  approached, 
although  these  results  have  not  been  experimentally  confirmed. 
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Two  papers,  [12]  and  [13],  have  analyzed  the  effects  on  the  small 
signal  gain  of  a  divergent  optical  wave,  and  found  two  deviations 
from  the  Madey  theorem  result:  a  shift  of  the  gain  curve  combined 
with  a  distortion  of  the  spectrum  for  large  enough  optical 
divergence.  This  paper  summarizes  the  results  of  the  first 
measurement  [1]  of  the  shift  in  resonance  parameter  between  the 
gain  spectrum  observed  on  a  diverging  wave  and  the  forward 
spontaneous  spectrum. 

To  first  approximation,  the  effect  of  the  divergence  of  the 
wave  is  that  a  phase  shift  appears  on  axis  as  the  wave  passes 
through  its  focus.  The  part  of  this  phase  which  is  linear  in  z 
produces  a  local  change  in  the  wavelength  of  the  probe  beam  which 
shifts  the  spectrum  of  the  gain  with  respect  to  the  Madey  result. 
The  nonlinear  part  of  the  phase  shift  distorts  the  shape  of  the 
curve.  For  large  divergence,  the  gain  curve  has  no  simple 
relation  to  the  spontaneous  spectrum;  the  Madey  theorem  is 
maximally  violated.  We  have  been  able  to  measure  the  resonance 
shift  induced  by  a  moderately  divergent  Gaussian  TEMoo  mode  under 
conditions  where  the  quadratic  and  higher  order  terms  in  the 
phase  shift  are  negligible. 

We  have  found  the  spontaneous  and  the  gain  spectra  for  the 
optical  klystron  including  the  effects  of  the  divergence  of  the 
wave  and  the  finite  size  of  the  electron  beam  in  [1].  Using  the 
notation  of  Elleaume  [14], 
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where  the  shift  in  the  resonance  parameter  2rrs  is  defined  in 
terms  of  the  lengths  of  the  first  undulator  L  and  the  dispersive 
section  d,  the  Rayleigh  range  zo  ,  and  the  standard  filling  factor 
Ff  [16] 
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The  gain  and  spontaneous  spectra  are  measured  as  a  function  of 
gap  through  filters  of  different  wavelengths.  If  we  define  6g  as 
the  absolute  shift  of  a  gain  peak  from  its  corresponding 
spontaneous  peak,  and  Ag  as  the  local  peak  to  peak  change  in  the 
gap,  the  relative  fringe  shift  5g/Ag  is 


If  -  J  +  s  +  -T  Ml 

The  first  term  of  this  expression  is  the  Madey  theorem  result. 

The  two  terms  included  in  the  shift  s  are  the  on-axis  phase  shift 
and  the  correction  term  due  to  the  reduced  phase  shift  of  the 
off-axis  electrons.  The  last  includes  the  effect  of  the  different 


measurement  wavelengths  for  the  two  spectra. 
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The  Experiment 

The  principal  difficulty  in  measuring  the  resonance  shift 
of  the  spontaneous  and  the  gain  spectra  is  the  small  size  of  the 
effect  compared  to  the  1/2N  relative  width  of  the  spectra.  The 
shift  of  an  optimized  TEMoo  mode  is  difficult  to  identify  in  the 
presence  of  measurement  noise.  This  difficulty  is  alleviated  if 
the  experiment  is  done  with  an  optical  klystron  since  the  fine 
structure  introduced  onto  the  spectra  increase  the  precision  of 
the  measurement.  Our  experiment  was  done  with  the  optical 
klystron  modification  [17]  of  the  undulator  NOEL. 

The  measurement  of  the  gain  [18]  requires  the  alignment  of 
an  external  laser  coaxial  with  the  electron  beam.  It  is 
preferable  to  measure  the  spontaneous  spectrum  simultaneously  in 
order  to  eliminate  the  various  sources  of  drift  in  the 
experiment.  However,  the  laser  beam  which  emerges  perfectly 
aligned  with  the  spontaneous  radiation  must  be  filtered  out  [1] 
to  allow  the  variations  in  the  weak  spontaneous  power  to  be 
observed.  This  is  done  by  shifting  the  analyzing  monochrometer 
away  from  the  laser  wavelength.  Although  the  difference  in 
wavelength  will  introduce  another  term  in  the  resonance  shift  of 
the  spectra,  the  effect  of  this  extra  term  can  be  subtracted 
during  the  analysis  of  the  data. 

The  resonance  shift  was  measured  for  two  values  of  the 
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probe  laser  beam  divergence.  The  first  experiment  was  performed 
after  focussing  the  laser  into  the  center  of  the  optical  klystron 
with  a  moderately  large  divergence  characterized  by  the  Rayleigh 
range  zo  =  ttw5/\  =  1.1  ±  .1  m.  For  the  second  experiment,  the 
laser  mode  was  configured  essentially  as  a  plane  wave  with  a 
considerably  smaller  divergence  and  a  Rayleigh  range  zo  =  2.3  ± 

.1  m . 


For  each  scan  the  resonance  shift  is  determined  from  the 
locations  of  the  peaks  and  the  valleys  by  multiple  linear 
regression  to  two  parallel  lines  in  gap  vs.  peak  number.  The 
shift  is  normalized  as  in  (4)  to  the  mean  peak  separation.  In  the 
high  divergence  experiment,  four  scans  were  taken,  two  at  5130  A 
which  showed  a  shift  of  .09  and  .11,  and  two  at  5160  A  with  a 
shift  of  .65,  and  .66.  The  mean  of  these  values,  and  therefore 
the  shift  at  the  laser  wavelength  5145  A,  is  . 38  ±  .01.  The 
theoretical  value  for  this  divergence  is  .32  ±  .02,  which  is 
composed  of  the  following  parts:  .32  =  .25  (Madey  Theorem)  +  .11 
±  .01  (on  axis  divergence  effect)  -  .04  ±  .01  (finite  beam  size 
correction  of  divergence  effect),  where  the  divergence  effects 
are  calculated  using  the  measured  values  of  the  electron  beam 
size  cr  2:  320  n  ±  20%.  The  agreement  between  the  theoretical 
prediction  and  the  experimental  value  is  reasonably  close  (see 
the  discussion  of  the  systematic  error  below) . 


The  wavelength  detuning  induced  resonance  shift  which 
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MEASURED  THEORETICAL 
resonance  result 
parameter 
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Table  I 


appears  in  the  above  data  is  6s  =  1.8xio~26X.  This  is  identical 

to  the  value  deduced  theoretically  from  the  measured  value  of  N^ 

=  90:  6s  =  (N+N^SX/X  =  1.9xlO~26X  ±  5%. 
a 

In  the  low  divergence  experiment,  four  successive 
measurements  were  taken  at  one  wavelength  detuning.  The  mean 
value  of  the  measured  shift  is  .09  ±  .01.  This  is  to  be  compared 
to  a  theoretical  value  of  .13  ±  .02  =  .25  -  .15  ±  .01  (wavelength 
detuning)  +  .05  ±  .01  (divergence  effect)  -.02  (beam  size 
correction)  where  the  electron  beam  size  is  u  2  350  n  ±  203s. 


The  small  deviations  from  the  theory  in  the  above  results 
are  significant  at  two  standard  deviations,  and  arise  from  some 
systematic  effect.  We  think  the  problem  arises  from  the  imperfect 
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characterization  of  the  phase  front  curvature  of  the  probe  laser 

beam.  The  Strehl  ratio  was  found  to  be  S  .  9  by  measuring  the 

laser  beam  waist  at  several  positions.  Since  a  light  beam  is 

considered  diffraction  limited  for  S  >  .8,  this  beam  is  of  very 

high  quality.  The  RMS  wavefront  deviation  c  is  given  from  the 

w 

Strehl  ratio,  for  small  a^/X,  by  [19] 
aw  I - 

2*-f  =  v!  1-S  (5) 

which  yields  <r  /X  =*  .05  in  our  case.  The  actual  evolution  of  the 
electromagnetic  phase  along  the  trajectory  of  an  electron 
therefore  contains  an  additional  component  due  to  the  evolution 
of  the  wavefront  distortion.  If  the  laser  beam  is  not  too  badly 
aberrated,  the  phase  front  distortion  is  slowly  varying  in  the 
transverse  dimensions,  and  all  electrons  see  approximately  the 
same  effect.  We  expect  under  these  circumstances  that  the  optical 
wave  will  shift  phase  by  an  additional  amount  on  the  order  of 
equation  (5)  as  the  probe  bear,  goes  through  its  focus,  which 
means  that  the  shift  s  in  (3)  contains  an  unknown  additional 
component  on  the  order  of  ±cr^/X.  This  sets  a  fundamental  linit  on 
the  precision  of  this  type  of  experiment.  The  errors  observed 
between  the  theoretical  and  experimental  values  of  Table  I  are  on 
the  order  of  .05,  in  agreement  with  our  hypothesis. 

Although  the  high  divergence  maximal  violation  of  the  Madey 
theorem  has  been  beyond  our  reach  for  this  experiment,  we  have 
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been  able  to  observe  the  linear  resonance  shift  which  occurs  like 
a  warning  signal  just  before  the  strong  distortion  of  the  gain 
curve.  We  have  shown  experimentally  that  the  Madey  theorem  is 
violated  for  divergent  optical  modes. 
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ABSTRACT 

We  have  simultaneously  measured  the  spectra  of  the  gain  an 
the  forward  spontaneous  radiation  in  the  Orsav  free  electron 
laser  so  that  we  can  extract  the  relative  shift  of  these  two 
curves.  Our  results  are  in  agreement  with  the  Madey  theorem  in 
the  low  divergence  limit,  and  demonstrate  the  violation  of  the 
theorem  produced  as  the  Rayleigh  range  of  the  wave  becomes 
comparable  to  or  smaller  than  the  length  of  the  interaction 
region.  We  extend  the  single  particle  theory  to  account  for  the 
divergence  of  the  wave  to  first  order.  The  predictions  from  this 
linearized  model  are  in  reasonable  agreement  with  the 
experimental  results  in  view  of  the  fundamental  limitations  on 
this  kind  of  experiment. 
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The  Xadey  theorem  '1]  has  proved  to  be  a  valuable  guide  to 
research  with  free  electron  lasers.  Although  the  theorem  can  be 
divided  into  two  parts  and  used  in  different  contents,  we  will 
consider  the  form  of  the  theorem  with  greatest  experimental 
value:  the  gain  spectrum  is  related  to  the  derivative  of  the 
forward  soonfaneous  emission  soectrum . 
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where  dW { y  )  /dfidu  is  the  energy  radiated  per  unit  frequency 
interval  per  unit  solid  angle  in  the  forward  direction  by  one 
electron  in  a  pass  through  a  linearly  polarized  magnetic 
structure.  The  Madey  theorem  simplifies  the  calculation  of  the 
gain  and  the  induced  energy  spread  in  an  arbitrary  magnetic 
structure  [1,2,3],  and  shows  that  the  easily  measured  spontaneous 
spectrum  is  a  reliable  diagnostic  of  the  small  signal  gain  [4], 
which  is  a  difficult  quantity  to  measure.  In  view  of  the  value  of 
this  theorem,  it  is  important  to  identify  its  limitations  and 
where  oossible  confirm  them  exoerimentallv. 


The  derivations  of  the  theorem  [1,2, 5, 6, 7]  all  assume  low 
gain,  an  unsaturatec  electron-photor.  interaction,  and  a  plane 
wave  electromagnetic  field.  We  have  shown  experimentally  [3]  thai 
the  theorem  is  obeyed  when  the  above  conditions  are  satisfied.  A 
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number  of  authors  !see  for  example  references  [3],  and  1 1C’)  have 
derived  the  theoretical  result  that  the  gain  spectrum  distorts  as 
the  saturated  or  the  high  gain  regime  is  approached,  although 
these  results  have  not  been  experimentally  confirmed.  Two  papers, 
[11]  and  [12],  have  analyzed  the  effects  on  the  small  signal  gain 
of  a  divergent  optical  wave,  and  found  two  deviations  from  the 
Madev  theorem  result:  a  shift  of  the  gain  curve,  and  a  distortion 
of  the  spectrum  for  large  enough  optical  divergence.  Our  paper 
reports  the  results  of  the  first  measurement  of  the  shift  in 
resonance  parameter  between  the  gain  spectrum  observed  on  a 
diverging  wave  and  the  forward  spontaneous  spectrum.  Since  the 
date  of  this  work,  the  effects  of  divergence  on  the  spectrum  of 
the  Crsay  oscillator  have  also  been  observed  and  reported  [13,7], 
This  latter  work,  however,  suffers  from  the  difficulty  that  the 
input  mode  content  is  not  measurable  so  that  a  numerical 
comparison  with  the  theory  is  not  possible. 

To  first  approximation ,  the  effect  of  the  divergence  of  the 
wave  is  that  a  phase  shift  appears  on  axis  as  the  wave  passes 
through  its  focus.  The  part  of  this  phase  which  is  linear  in  z 
produces  a  local  change  in  the  wavelength  of  the  probe  beam  which 
shifts  the  spectrum  of  the  gain  with  respect  to  the  Madey  result. 
The  nonlinear  part  of  the  phase  shift  distorts  the  shape  of  the 
curve.  For  large  divergence,  the  gain  curve  has  no  simple 
relation  tc  the  spontaneous  spectrum;  the  Madey  theorem  is 
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=  2~m  -  - 


The  gain  and  spontaneous  spectra  are  measured  as  a  function  ">* 


gap  through  fi! 

(G)  (S) 

gm  -g.-x  as 


:ers  of  different  wavelengths.  If  we  define  5g  - 
:he  absolute  shift  of  a  aain  oeak  from  its 


corresponding  spontaneous  peak,  Ag  =  g.-i~gjr.+  :  as  the  local  peak  to 
peak  change  in  the  gap,  and  3X  =  X'  '  as  the  difference  in 

the  measurement  wavelengths  of  the  spontaneous  from  the  gain 
spectra,  the  relative  fringe  shift  Sg/Ag  can  be  found  by 
linearizing  equations  (15) 


(M-Nd)  ~ 


The  first  term  of  this  expression  is  the  Madey  theorem  result. 

The  two  terms  included  in  the  shift  s  are  the  on-a:-:is  phase  shift 
produced  in  a  Gaussian  mode,  and  the  correction  term  which 
accounts  for  the  reduced  phase  shift  experienced  by  the  off-axis 
electrons.  The  last  term  accounts  for  the  fact  that  the 
spontaneous  and  the  gain  spectra  are  measured  at  different 
wavelengths . 
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of  the  spontaneous  and  the  gain  spectra  is  the  small  size  of  the 
effect  cor.parec  to  the  1/2”  relative  width  of  the  spectra.  The 
shift  of  an  optimized  TEMoo  mode  is  difficult  to  identify  in  the 
presence  of  measurement  noise.  This  difficulty  is  alleviated  if 
the  experiment  is  done  with  an  optical  klystron  since  the  fine 
structure  introduced  onto  the  spectra  increase  the  precision  of 
the  measurement.  Our  experiment  was  done  with  the  optical 
klystron  modification  IIS]  of  the  undulator  NOEL. 


The  comparison  of  two  subsequent  scans  of  the  gain  or  the 
spontaneous  emission  spectra  is  complicated  by  the  number  of 
extraneous  parameters  which  must  be  held  constant,  and  in 
principle  measured.  Small  resonance  shifts  of  the  spectrum  can 
easily  occur  if  the  beam  energy  changes  between  scans  (the  magnet 
current  control  shunts  are  not  temperature  stabilized  on  ACC);  if 
the  electron  orbit  changes  (due  to  energy  or  magnetic  focussing 
or  RF  frequency  drifts);  or  if  there  is  play  in  the  mechanical 
gap  adjustment  system.  It  is  preferable  to  measure  the  two 
spectra  simultaneously. 


?he  measurement  of  the  cain  '19’  recuires  the  alignment  of 


an  external  laser  coaxial  with  the  electron  beam.  If  the 
spontaneous  emission  measurement  is  to  be  oer formed 


simultaneously ,  tr.e  laser  Dear,  emerges  perrect_y  a^zgr.ec  w ;• 


spontaneous  radiation  and  may  interfere  with  the  measurement 
its  spectrum.  We  use  spectral  filtering  (figure  1)  to  make 
oossible  the  simultaneous  measurements. 


A  semi  transparent  mirror  splits  off  some  of  the  radiatio 
for  analysis  by  the  gain  electronics,  which  performs  its  funct 
without  modification.  The  remainder  of  the  light  is  directed  t 


the  spontaneous  emission  diaar.ostic: 


the  mor.ochromato 


is  detuned  from  tr.e  Argon  laser  line  by  a  sufficient  margin,  i 
is  possible  to  filter  out  the  laser  light  while  allowing  the 
variations  in  the  relatively  weak  spontaneous  power  to  be 
observed.  Although  the  difference  in  wavelength  will  introduce 
another  term  in  the  resonance  shift  of  the  spectra,  the  effect  of 
this  extra  term  can  oe  subtracted  during  the  analysis  of  the  data 


The  gain  measurement  is  set  up  as  usual  on  the  5145  A  line 
of  the  Argon  laser.  The  monochromator  for  analyzing  the 
spontaneous  radiation  is  detuned  alternately  to  the  short  and  the 
long  wavelength  side.  Scans  are  taken  as  a  function  of  the 
magnetic  gap,  overshooting  enough  on  each  side  to  eliminate 
mechanical  backlash.  The  data  is  recorded,  the  locations  of  the 


means  ana 


va__eys  are  ceterr.ir.ad,  c.nc  tr.e  relative  shift  o: 


r.e  gain  spectrum,  is  calculated  oy  a  linear  regression  analysis 


Since  the  vaveler.cth  is  disnlaced  bv  the  same  amount  in  both 
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directions,  the  resonance  parameter  shift  at  the  laser  wavelength 
can  be  determined  simoly  as  the  mean  of  the  two  shifts  at  the 
displaced  wavelengths  (see  figure  2). 

The  laser  mode  matching  was  performed  by  measuring  the 

laser  mode  size  at  several  positions  with  the  aid  of  a  Ret i con 

diode  array.  Alignment  of  the  electron  beam  was  obtained  to  a 

precision  of  about  SC  urn  by  optimizing  the  measured  gain  with  the 

aid  of  four  horizontal  and  -vertical  steering  coils.  The 

experiment  was  performed  at  relatively  low  currents  so  that  the 

product  of  the  modulation  depth  of  the  spontaneous  emission  and 

the  signal  level  of  the  gain  would  be  a  maximum.  The  electron 

beam  size  was  also  measured  with  the  Reticon.  The  optical 

klvstron  tarameter  IT.  was  measured  from  a  wavelength  scan  of  the 

a 

spontaneous  emission  at  the  central  value  of  the  magnetic  gap. 

The  fringe- to- fringe  spacing  AX  determines  this  quantity  through 

**  .  *T  _  *>  '  A  > 


The  resonance  shift  was  measured  for  two  values  of  the 
probe  laser  beam  divergence.  The  first  experiment  was  performed 
after  focussing  the  laser  into  the  center  of  the  optical  klystro 
with  a  moderately  large  divergence  characterized  by  the  Rayleigh 
range  z;  a  ~w 8/X  =  l.l  -  .1  m.  For  the  second  experiment,  the 
laser  mode  was  configured  essentially  as  a  plane  wave  with  a 
considerable*  smaller  divercence  and  a  Ravleigh  ranee  zc  =  2.3  ± 


resonance 

parameter 


shift 


f. 05  =.25  -.27 


3  3  ^ 


1.66 


j  ( 5 1 30A  ) 

I 

1.59  =. 25  +. 27  +.11 

(516CA) 


1.1  n  320  cm 


- .  04 


,  C9 


-.  15 


3  m  3  5  0  cm 


(  5 1 3  7  A  ) 


Table  I 


Tor  each  scan  the  resonance  shift  is  determined  from  the 
locations  of  the  peaks  and  the  valleys  by  multiple  linear 
regression  to  two  parallel  lines  in  gap  vs.  peak  r.umoer.  The 
shift  is  normalized  to  the  mean  peak  separation.  In  the  high 
divergence  experiment,  four  scans  were  taken,  two  at  5130  A  which 
showed  a  shift  of  .09  and  .11,  and  two  at  5160  A  with  a  shift  of 


.65,  and  .66.  The  mean  of  these  values,  and  therefore  the  shift 
at  the  laser  wavelength  5145  A,  is  .33  =  .01.  The  theoretical 


value  for  this  divergence  is  .32  ±  .02,  which  is 
following  parts:  .32  =  .25  (Madey  Theorem)  -  .11 
divergence  effect)  -  .04  =  .01  (finite  beam  size 
divergence  effect),  where  the  divergence  effects 
using  the  measured  values  of  the  electron  beam  s 


composed  of  the 
-  .01  (on  axis 
correction  of 
are  calculated 
ze  r  ^  320  _m  = 


I 
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2C"o.  The  agreement  between  the  theoretical  predict 
experimental  value  is  reasonably  close  (see  the  discussion  of  the 
svstematic  error  beicwi . 


.  r.e  wave.; 


ununc  rnuucea  resonance  sruft  wr.icn 


aooears  in  the  above  data  is  -3s  =  1.3x10' 


hors  is  icsntica. 


to  the  value  deduced  theoretical Iv  from  the  measured  value  of  N . 

a 

=  90:  3s  =  (::->!.}  S\/X  =  1.9xl0“-5\  ±  5?<5. 
c ' 


In  the  lev;  divergence  experiment,  which  was  performed 
almost  one  year  following  the  first  measurements,  it  was  no 
longer  possible  to  take  simultaneous  measurements.  Because  of  the 
now  large  angular  fluctuations  of  the  probe  laser,  the  noise 
level  contributed  by  the  laser  in  the  spontaneous  measurement 
completely  dominated  the  spontaneous  signal.  The  best  which  could 
be  done  was  to  take  a  series  of  measurements  contiguous  in  time 
and  alternating  between  gain  and  spontaneous  emission.  Seven  such 
measurements  were  taken  at  low  stored  current  so  that  the  beam 
lifetime  would  be  long.  After  normalizing  for  the  varying  scan 
rates  of  the  optical  klystron  magnetic  gap,  we  found  that  these 
seven  measurements  are  consistent  internally.  This  indicates  that 
the  uncontrolled  external  drifts  were  small  during  the  time  scale 
of  the  measurements,  and  we  have  been  able  to  extract  reliable 
information  even  in  a  set  of  non-simul taneous  experiments. 
Although  these  measurements  were  taken  at  only  one  wavelength 


cettninc.  the  behavior  of  the  orevious  results  verifies  that  we 
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know  how  to  calculate  the  resonance  parameter  shift  due  to  a 
wavelength  detuning.  Tour  successive  measurements  of  the  shift 
give  a  mean  of  .C9  =  .31.  This  is  to  be  compared  to  a  theoretics 
value  of  .13  ±  .02  =  .25  -  .15  =  .01  (wavelength  detuning)  -r  .05 
±  .01  (divergence  effect)  -.02  (beam  sire  correction)  where  the 
electron  beam  sire  is  r  ^  350  -m  =  20%. 

Again,  the  agreement  is  reasonable.  We  note,  however,  that 
there  is  a  small  deviation  of  the  measurements  towards  the 
extremes  which  may  be  significant.  The  high  divergence 
measurements  show  a  slightly  larger  shift  and  the  low  divergar.ee 
results  shew  a  smaller  shift  than  the  theory  predicts. 


These  deviations 

from  the  theory 

may  be  caused 

by 

the 

imperfect  character i rat 

ior.  of  the  phase 

front  cur vat 

ure 

of  the 

probe  laser  beam.  The  3 

trehl  ratio  was 

found  to  be  S 

^  , 

.9  by 

measuring  the  laser  beam  waist  at  several  positions.  Since  a 

light  beam  is  considered  diffraction  limited  for  S  b  .3,  this 

bean  is  of  very  high  cualitv.  The  RMS  wavefront  deviation  r  is 

w 

given  from  the  Strehl  ratio,  for  small  c  X,  bv  [20’ 

w 


which  yields  c  'X  ^  .05  in  our  case.  The  actual  evolution  of  the 
electromagnetic  phase  along  the  trajectory  of  an  electron  is  the: 
given  by  (5)  with  an  additional,  component  due  to  the  evolution  o 
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the  wavefront  distortion.  If  the  laser  beam  is  not  too  badly 


distorted , 

the  distortion 

is  slowly  varying  in  the 

trar.s' 

verse 

dimensions 

,  and  all  elect 

ror.s  see  approximately  the 

same 

effect 

We  expect 

under  these  cir 

cunstar.ces  that  the  linear 

part 

of  the 

distorted 

comoonent  will 

shift  chase  bv  an  amount  o 

v.  .  -  tr 

order 

equation  (17)  as  the  probe  beam  goes  through  its  focus,  which 

means  that  the  shift  s  in  (13)  will  have  an  unknown  component  on 

the  order  of  =c  This  sets  a  fundamental  li/nit  on  the 

w 

precision  of  this  type  of  experiment.  Indeed,  the  errors 
observec  between  the  theoretical  and  experimental  values  of  Table 
I  are  on  the  order  of  .35. 


In  summary,  we  have  found  in  this  experiment  that: 

a)  within  the  anticipated  error  bars,  the  divergence  of  the 
optical  wave  produces  a  resonance  parameter  shift  of  sign  ar.d 
magnitude  predicted  by  the  linearized  small  signal  theory;  and 

b)  a  sufficiently  low  divergence  beam  shows  no  shift,  in 
agreement  with  the  Maaey  theorem. 


Although  the  high  divergence  maximal  violation  of  the  Made\ 


heorem  has  been  beyond  our  reach  for  this  experiment,  we  have 


been  able  to  observe  the  linear  resonance  shift  which  occurs  like 
a  warning  signal  just  before  the  strong  distortion  of  the  gain 
curve.  We  have  shown  experimental ly  that  the  Maaey  theorem  is 
violated  for  divergent  optical  modes.  Any  application  of  “he 
theorem  must  test  for  low  divergence 
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gure  Captions 


gure 


Scfier.atic  diagram  of  the  experiment.  The  gain  spectrum 
[8,133  -s  measured  simultaneously  with  the  forward 
spontaneous  emission  spectrum  [15].  In  the  spontaneous 
measurements,  the  laser  power  is  filtered  out  by 
setting  the  bandpass  of  the  monochrometer  at  a 
wavelength  displaced  slightly  from  the  Ar  line  (see 
Table  I ) . 


gure  2 


The  gain  spectrum  shown  in  the  central  trace  is 
clearly  displaced  from  the  two  spontaneous  curves.  The 
peak  positions  of  a  spontaneous  spectrum  taken  at  5145 
A  would  fall  midway  between  those  of  the  top  and 
bottom  spontaneous  curves.  The  shift  of  the  gain  peaks 
in  a  low  divergence  experiment  is  almost  exactly  .25 
of  the  fringe  spacing  as  predicted  by  the  Xadey 
theorem,  once  the  wavelength  shift  is  factored  out. 
This  figure  shows  an  example  of  the  increase  of  the 
shift  (to  .38)  as  a  result  of  the  violation  of  the 
Madey  theorem  induced  by  the  divergence  of  the  input 
=  1.1  m  in  this  case) . 
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Abstract.  We  derive  the  most  general  equations  of  motion  for  the  electrons  and  the 
electromagnetic  field  in  a  free-electron  laser  including  the  effects  of  diffraction  and  pulse 
propagation.  The  field  evolution  is  expressed  in  terms  of  the  amplitudes  and  phases  of  a 
complete  set  of  transverse  modes.  The  analytic  solution  is  given  in  the  small-signal  regime, 
where  the  theory  is  shown  to  be  in  excellent  agreement  with  a  recent  experiment  at  Orsay. 

PACS:  42.60,  42.20,  42.55 


Stimulated  by  the  original  free-electron-laser  experi¬ 
ments  [1]  in  1977,  a  number  of  authors  have  contri¬ 
buted  to  the  development  of  a  purely  classical  theory 
for  the  electron  dynamics  and  the  electromagnetic 
wave  growth  in  these  devices.  The  initial  work  as¬ 
sumed  the  light  could  be  represented  by  a  single¬ 
frequency  plane  wave  [2,  3],  The  first  generalization 
was  required  to  explain  the  extremely  short  pulse 
phenomena  observed  at  Stanford  [4,  5],  The  inclusion 
of  the  longitudinal  modes  in  the  theory  [6-8]  per¬ 
mitted  the  explanation  of  the  cavity  detuning  curve, 
and  predicted  a  range  of  phenomena  in  the  pulse 
structure  which  have  yet  to  be  observed.  More  re 
cently,  the  theory  has  been  broadened  to  include  the 
transverse-mode  structure  of  the  optical  beam  [9-15]. 
Until  our  work  at  Orsay  [16],  no  experimental  infor¬ 
mation  has  been  available  to  test  the  validity  of  these 
so-called  3D  theories. 

In  this  paper,  we  present  a  new  approach  to  calculating 
the  three-dimensional  effects  operative  in  free-electron 
lasers.  The  previously  mentioned  approaches  consider 
the  growth  of  the  field  r?(r,  r)  along  the  propagation  or 
r  axis  by  evaluating  its  change  at  each  point  (  v,  y),  and 
integrating  numerically  through  the  interaction  region 
in  the  time  domain  [9-14]  or  in  the  frequency  domain 
[15].  These  techniques  all  demand  long  computer  runs 
if  they  are  to  be  applied  to  a  real  experimental 
situation.  Our  approach  decomposes  the  problem  into 


the  minimum  number  of  physically  observable  quan¬ 
tities  :  the  transverse  optical  modes  of  the  system.  The 
field  evolution  is  expressed  in  terms  of  a  complete  set  of 
orthogonal  transverse  modes;  equations  are  developed 
for  the  propagation  of  the  amplitude  and  phase  of  each 
mode.  In  physical  systems  which  operate  on  a  few  of 
the  lowest-order  modes,  this  approach  greatly  in¬ 
creases  the  accuracy,  and  may  reduce  the  required 
computer  time  for  the  calculation  by  working  in  a 
vector  space  well  matched  to  the  solution  of  the 
problem.  For  the  oscillator  case,  the  appropriate 
choice  of  modes  is  the  set  of  eigenmodes  of  the  cavity. 
For  the  amplifier,  the  vector  space  of  modes  is  de¬ 
termined  by  the  characteristics  of  the  input  mode, 
which  is  presumably  a  TEM00  Gaussian  mode.  In 
either  device,  an  optimum  design  would  result  in  the 
excitation  of  as  few  of  the  higher-order  modes  as 
possible.  The  modal  decomposition  method  is  there¬ 
fore  well  adapted  to  the  prediction  and  optimization  of 
the  operation  of  the  free-electron  laser  (FEL). 

In  the  first  section,  we  derive,  in  their  most  general 
form,  the  equations  governing  the  dynamics  of  the 
complex  mode  amplitudes.  The  subsequent  sections 
reduce  these  equations  to  the  familiar  case  of  the  small 
signal,  low-gain  result  (Sect.  2).  Here,  the  problem 
becomes  linear,  the  mode  evolution  can  be  described 
by  a  matrix  transformation,  and  we  retrieve  the  well 
known  gain  equation  complete  with  filling  factor.  The 
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theory  is  then  applied  to  the  case  of  the  Orsay  experi¬ 
ment,  where  the  results  are  in  excellent  agreement  with 
an  experiment  [16]  performed  recently  which  exhibits 
the  off-diagonal  terms  of  the  gain  matrix. 


I.  Theoretical  Development 
of  the  Fundamental  Equations 

The  FEL  system  is  properly  described  by  the  coupled 
Maxwell  and  Lorentz  force  equations.  From  these,  we 
shall  derive  a  self-consistent  set  of  equations  describing 
the  electron  and  the  transverse  optical  mode  dynamics. 
We  use  the  dimensionless  notation  originally  deve¬ 
loped  by  Colson  (in  fact  this  work  is  a  generalization  of 
Colson's  work  to  include  transverse  modes  and  we 
shall  stay  as  close  as  possible  to  his  original  notation). 
Let  us  recall  his  main  equations  describing  the  field 
and  electron  dynamics  in  the  slowly  varying  phase  and 
amplitude  approximation  [18]: 

dv 

—  =a  cos(s  +<j>),  (1) 

at 


(2) 


dd_ 

dr 


(3) 


where 


C(t)  =  (k  +  k0)r(f)-cur 

is  the  dimensionless  electron  phase. 

v{t)  =  L[(k  +  k0)fl:(t)-k] 

the  dimensionless  resonance  parameter, 

_  ct 
T=  L 

the  dimensionless  interaction  time. 


a'(z.t)  = 


4neN  LKE(z.  f)e'*<=  l 


the  dimensionless  complex  field  amplitude,  and 

sun,.*™??*;**”. 


the  dimensionless  gain  parameter. 


(4) 

(5) 

(6) 

(7) 

(8) 


Here  we  consider  an  N-period  helical  undulator  of 
length  L,  magnetic  period  /.r)  =  2 n  k0,  peak  magnetic 
field  B.  and  deflection  parameter  K  =93.4  B 
[Gauss] [cm],  An  electron  beam  of  energy  ymc*, 
and  number  density  g  travels  along  the  axis  of  the 
undulator;  an  individual  electron  has  longitudinal 
coordinate  r(t)  and  longitudinal  velocity  c/L(t)  at  time 


r.  A  helically  polarized  plane  ware  of  wavelength 
/.  =  2n;k,  frequency  to.  and  electric  field  £(:,  r)  =  E(:,  t) 
exp{i[kz-ajr  +  <fi(z,t)]}  interacts  with  the  electrons.  In 
(3),  <  >-oVo  is  the  average  over  the  initial  phase  C0  and 
resonance  parameter  v0  of  the  electron  population  at 
the  position  ;. 

Equations  (1,2)  are  derived  directly  from  the  Lorentz 
force  equation  and  describe  the  effect  of  the  radiation 
field  on  the  electrons.  The  work  done  by  the  longitu¬ 
dinal  field  on  the  electrons  is  neglected  here,  which  is  a 
good  approximation  provided  that  the  modes  are  not 
too  divergent  [14]  /J(o0<^2it2K.\  y.  Equation  (3)  is 
derived  from  the  Maxwell  equations  and  describes  the 
effect  of  the  electron  on  the  radiation  field.  The  set  (1), 
(2).  and  (3)  is  self-consistent.  Indeed,  those  equations 
are  very  close  to  being  the  most  general  classical 
equations  describing  the  FEL  dynamics.  They  apply  to 
high  and  low  gain  devices  (r'£>  1  or  r  <?  1),  high  field 
and  low  field  cases  (a'?>l  or  «' <g  1  >,  and  include  the 
effects  of  multiple  longitudinal  modes  (laser  lethargy 
effects)  through  the  £  dependence  of  r\  E.  and  <p.  Slight 
modifications  allow  their  extension  to  the  cases  of: 

-  the  planar  undulator  [19], 

-  the  tapered  undulator  [18], 

-  the  optical  klystron  [14],  and 

-  space  charge  effects  [19]. 

However,  the  plane-wave  approximation  cannot  ac¬ 
curately  describe  the  transverse  effects  produced  by  the 
finite  transverse  extent  of  the  optical  mode  and  the 
electron  beam.  A  filling  factor  calculated  with  an 
ad-hoc  overlap  integral  can  be  added  to  the  results  of 
this  calculation,  and  gives  satisfactory  results  in  the 
small  signal  regimt  only  so  long  as  one  is  not  interested 
in  the  exact  transverse  field  profile. 

To  relieve  this  last  restriction  on  the  theory,  we  assume 
the  field  to  be  described  in  free  space  by  the  paraxial 
wave  equation  [20] : 

/  f2  f2  e  t 

U-5  +  —  -2ifc—  £<r.f)e‘*''-"  =  0.  (9) 

\cx  cy-  czj 

This  equation  is  derived  from  the  wave  equation 
(V2  —  c_2t,2/«V)<f  =0  in  the  slowly  varying  amplitude 
and  phase  approximation,  and  has  been  widely  used  in 
laser  field  calculations  [17,  20].  The  general  solution  of 
(9)  can  be  expressed  as  a  linear  combination  of  a 
complete  set  of  orthogonal  modes.  If  we  define  these 
modes  by  the  complex  amplitude  £mexp(it/’m).  where 
Em  is  real  and  m  is  the  generalized  index  of  the  mode  (in 
the  two-dimensional  transverse  space  we  consider,  m 
represents  two  integer  numbers),  the  most  general 
expression  for  the  field  is 

E(r.  f  (c’*""  =  X  rjf )  £m(r)e,v~"' ,  (10) 

m 

where  cm  is  complex  and  time-independent  in  free 
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space.  The  orthogonality  relation  reads 

.JxJv  =  (Ill 

t  j  /»i  n  mn 

n\\  o 

where  we  have  chosen  a  convenient  normalization 
which  makes  the  E„  dimensionless.  The  modes  can  be 
chosen  in  a  variety  of  symmetries,  but  it  is  useful  to 
exhibit  their  specific  form  in  cylindrical  symmetry: 


~p\  v\,j  ,  r  1  Jcos/0| 
\  (/>  +  /)!  ( I  -r- <>,„ I  «(:)’  w(rr  ) sin /t/ ( 


ll 


Ir 


(Id) 


(13) 


'■Wv’l.-l' 

u'rlir)=  iW)  -(2/,  +  /+1)tan  '  r- 

where  r  is  the  radial  and  0  is  the  azimuthal  coordinate. 
wn  is  the  beam  waist.  L‘[2r2  w’l  is  the  associated 
Laguerre  poly  nomial,  and 


w2(r)  =  w^(l  +  —  j, 

(14) 

\  -O'' 

R(r)  =  --(t+  72.). 

(15) 

7tW2 

-0  ~  ;  • 

(16) 

These  modes  are  very  useful  for  the  case  of  a  cylindri¬ 
cal  electron  beam  aligned  to  the  axis  of  the  light  beam. 
For  an  ellipsoidal  electron  beam  profile,  or  off-axis 
electron  injection,  the  rectangular  eigenmodes  are 
more  appropriate.  Although  we  will  use  the  cylindrical 
modes  in  the  examples,  we  proceed  with  the  general 
theoretical  development  which  makes  no  assumptions 
on  the  specific  form  of  the  modes. 

In  FEL.  the  coefficients  in  ( 10)  become  time  dependent. 
We  wish  to  calculate  the  evolution  of  the  amplitude 
and  phase  of  these  mode  coefficients.  Proceeding 
through  the  derivation  of  (1—3).  making  only  the  slowly 
varying  amplitude  and  phase  approximation,  but  now 
using  (10)  and  (11).  we  find 


—  =  X  l‘U  £r*  cos<; + 1 + </>  j , 

^  ^  m 

cX=v 

cx 

(  r  7TW  q 

where  we  have  made  the  new  definitions 
4ne\LK 

8  n2e2SL2K2 

r{rj)=  — — t - <?(r  f), 

y  me* 

CJ:-  fl  =  k'J"  fll  . 


(17) 

(18) 

(19) 

(20) 

(21) 

(22) 


As  before.  (17)  and  (18)  describe  the  effect  of  the 
radiation  field  on  the  electrons,  and  (19)  describes  the 
growth  or  decay  of  the  radiation  field  due  to  its 
interaction  with  the  electrons.  The  change  in  (17)  is 
quite  straightforward.  Equation  (19l  shows  clearly  the 
fact  that  the  growth  in  the  m,h  mode  amplitude  and 
phase  is  given  by  the  operlap  integral  of  the  inphase 
and  out-of-phase  components  of  the  charge  density 
w  ith  the  complex  conjugate  of  that  mode,  as  one  w  ould 
expect.  We  note  that  the  only  assumptions  made  on 
the  modes  £?m exp(ivm)  used  in  ( 17-19)  are  orthogo¬ 
nality  and  completeness.  This  means  these  equations 
are  also  valid  for  the  cases  of  waveguide  modes  and 
dielectrically  loaded  cavities.  In  this  case,  w0  is  no 
longer  the  mode  waist  in  the  usual  Gaussian  sense,  but 
is  defined  by  (11).  As  before,  these  equations  are  self- 
consistent.  An  example  of  this  fact  is  the  energy 
conservation  equation 


(23) 


which  is  derived  from  (17)  and  (19).  The  left-hand  side 
of  (23).  the  total  energy  gained  by  all  the  modes,  is 
equal  to  the  energy  loss  integrated  over  all  of  the 
electrons  in  the  beam. 

Equations  (17-19)  retain  all  of  the  generality  of  (1-3). 
They  are  valid  for  high  and  low  fields,  and  high  and 
low  gain  systems.  They  take  into  account  the  evolution 
of  the  transverse  modes  explicitly,  and  the  evolution  of 
the  longitudinal  modes  implicitly,  by  keeping  track  of 
the  i  dependence  of  the  charge  density  r(r,  t)  and  of  the 
mode  amplitudes  a(z,  t).  For  simplicity  in  the  follow  ing 
development,  we  drop  the  explicit  z  dependence  which 
has  been  thoroughly  discussed  by  Colson  [18],  and 
concentrate  on  the  transverse  phenomena. 

As  discussed  in  [19],  the  generalization  to  the  case  of 
the  planar  undulator  is  no  more  than  a  change  in  the 
definition  of  the  two  parameters 


IneNLKUJJ 
m  yhne1  m 


..  4n2e2NL2K2[Jjy 
r  =  ,  ,  t?(r.  f ) . 

yJmc‘ 


(24) 

(25) 


where 


[JT]  =  J0 


4  +  2  K2I 


K2 


4  +  2 K2  ' 


(26) 


Equations  (17-19)  can  be  integrated  numerically  to 
find  the  evolution  of  the  optical  wave  in  any  Compton 
regime  FEL.  In  a  high-field  experiment.  (18)  and  (19) 
are  nonlinear  in  a,  and  the  wave  evolution  can  only  be 
obtained  numerically.  In  this  case,  (17-19)  provide  a 
precise  and  efficient  technique  for  solving  the  general 
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problem.  In  a  low-field  situation  such  as  we  find  at 
Orsay.  however,  the  problem  becomes  linear,  and  can 
be  solved  analytically.  We  proceed  with  the  low  field 
case  in  the  next  section 


2.  The  Low-Field  Solution 

The  low-field  case  is  defined  by  |<jp|  1  for  every  mode. 

In  other  words,  the  electrons  do  not  become  over¬ 
bunched.  Experiments  which  operate  in  this  domain 
include  the  low -field  amplifier  experiments,  and  sto¬ 
rage  ring  FEL  oscillators  which  saturate  by  mecha¬ 
nisms  other  than  overbunching.  The  ignition  of  any 
FEL  oscillator  also  occurs  in  this  domain. 


2.1.  The  Gain  Matrix 


Equations  (17  19)  can  be  solved  b>  integrating  (17) 
and  (18)  to  lowest  order  m  the  fields  am  and  inserting 
the  result  for  '  into  (191.  If  the  electrons  are  uniformly 
distributed  initially  in  phase,  we  find 

=  f  dx  jj  t/r" Mm„ir.  r")it„lt") .  (27) 

n  o 

where 


•v^»l  r.O=yJ  ‘‘—I-  '1  v.  y)  EJ  x.  t .  r )  Ej  x.  y.  r") 

.  e  t)  -  v»U.  V.  i»nl«  f  'y  pfl) 


Equation  (27)  describes  a  linear  evolution  of  the  mode 
amplitudes,  and  upon  integration,  gives  the  relation 


<Ut=l)  =  (/+G)m„u„(t=0), 


(29) 


where  I  is  the  identity  matrix,  and  G.  which  is  generally 
not  Hermitian.  has  elements 

9mn  =  i^'i^",V/m„(T.T") 

0  0  0 

+  I  <1*1  j  dr 2  Jf  dtyUJ r,,Tj) 

0  0  0 

7  dri  / t/r,  'jdT^MJty  rj  +  ... .  (30) 

0  0  0 

The  higher-order  terms  in  ijmn  are  proportional  to  r 2 
and  higher  powers  of  r.  and  are  negligible  in  the  low- 
gain  case. 

Evidently  this  matrix  is  of  great  interest  since  multiple 
passes  of  the  electron  beam  will  result  in  multiple 
products  of  this  matrix,  greatly  simplifying  the  calcu¬ 
lation  of  the  modes'  growth.  We  shall  discuss  the 
consequences  for  an  oscillation  experiment  in 
Sect.  2.2. 


Let  us  note  that  this  gain  matrix  is  generally  complex 
and  defines  the  growth  of  the  amplitude  of  the  field. 
Sometimes  people  speak  of  the  gain  in  a  mode  "m"  as 
the  energy  gained  by  this  mode  in  a  pass  through  the 
undulator.  This  gam  is  simply  2  Re|l/mmi  Of 

course,  one  must  keep  in  mind  that  energy  is  radiated 
into  other  modes,  and  (hat  cross  terms  will  mix  a 
multiple  mode  input.  If  the  input  beam  is  truly  mo¬ 
nomode.  the  power  radiated  into  the  n,h  mode  is  lower 
than  that  into  the  m'h  mode  by  the  ratio 
-  R*-’  if/mm!  which  is  small  for  low  gain  ( r 1 ) 
systems.  It  is  only  in  this  case  that  it  makes  sense  to 
speak  of  the  gain  of  a  mode.  In  high-gain  systems, 
however,  the  off-diagonal  terms  can  lead  to  substantial 
emission  of  energy  into  the  higher-order  transverse 
modes.  If  the  input  beam  is  multimode,  of  course, 
mode  mixing  occurs  at  all  power  levels. 

Let  us  calculate  iimn  in  the  simple  case  of  experimental 
interest  where  the  electron  beam  is  cylindrical,  and  a 
good  choice  of  modes  is  the  cylindrical  cavity  eigen- 
modes  (12)  and  (13).  We  restrict  ourselves  to  the 
weakly  diverging  case  nw^P/.L  where  the  gain  takes 
on  its  most  familiar  form.  The  mode  amplitudes  and 
phases  in  (28)  become  independent  of  r,  and  we  can 
integrate  the  first  term  in  (30)  to  find  the  gain.  The 
average  over  the  resonance  parameter  in  (28)  becomes, 
under  the  assumption  of  a  Gaussian  distribution  of 
centroid  v  and  deviation  <rv 


<e 


r  ’>  =  e  - 

/  Vo 


(31) 


Under  the  weak-divergence  approximation,  and  as¬ 
suming  negligible  pulse  slippage  effects  (long  electron 
bunch  length  a ,  g>  /VA),  the  only  time-dependence  in 
(28)  is  that  of  (31).  For  small  spread  crv<t  1  the  integral 
gives  the  well  known  gain  spectrum 


.dxd\ 


—f  **.>•>  EJx,y)  E„(x,  v)  e'^-» 


1  -  cos  vr  -  -y  sin  vf  -  ^  -  j  cos  i;  +  sin  v 

- - . - +  i  — - - - 1 - 


(32 


In  the  usual  experimental  case  (unfortunately).  |ymn|  <?  1 
and  the  energy  gain  Gm  on  the  mode  m  becomes 


dxdy 
— ~  2r£' 
n"- 


-cosvc—  2  sm  >, 


(33! 


This  is  exactly  the  gain  one  calculates  by  using  the 
filling  factor  obtained  oy  integrating  the  mode  profile 
overlap  with  the  gain  profile.  Specializing  to  the 
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TEM(10  case  with  a  Gaussian  electron  beam  of  w  idth  <t, 
we  find 

v 

1  —  cos  vc  —  ~  sin  v . 


complete  with  the  familiar  filling  factor. 

The  v  dependence  of  G„  is  the  well  known  spectral 
dependence.  The  imaginary  part  of  is  not  new.  It 
describes  the  phase  shift  of  the  radiation  field  as 
described  by  Colson  [18].  The  inhomogeneous 
broadening  term  in  (31)  clearly  distorts  and  reduces  the 
magnitude  of  the  gain  spectrum  if  it  is  present  in  the 
integral  of  (30). 

The  effect  of  the  divergence  of  the  beam  on  the 
diagonal  terms  in  G  is.  to  first  order,  and  for  a 
filamentary  electron  beam,  the  addition  of  a  time- 
varying  phase  which  shifts  the  resonance  curve  in  (32) 
by  a  constant  depending  on  the  mode 


Gaussian  TF.VI,,,  modes,  and  the  phase  shift  per  round 
trip  ot,  becomes,  for  it  =  2  reflections  per  amplifica¬ 
tion  and  identical  radius  of  curvature  mirrors. 

=  4<2p  +  /  +  1 )  tan  ' (Z.,  2j„).  where  Lc  is  the  optical 
cavity  length. 

The  matrix  (36)  is  the  fundamental  matrix  of  the 
problem.  Its  diagonalization  allows  the  calculation  of 
the  mode  evolution  up  to  the  onset  of  saturation: 

(I  +  G)C  =  PAP .  (38) 

[(/  +  G)C]m  =  P.rP~  1 .  (39) 

where  the  columns  of  P  are  composed  of  the  eigenvec¬ 
tors  of  (I  +  G)C.  and  .t  is  diagonal.  The  fastest  mode 
growth  will  be  obtained  with  the  eigenmode  having  the 
highest  eigenvalue  modulus.  Optimization  of  the  FEL 
oscillator  will  then  consist  of  maximizing  the  desired 
eigenvalue  of  [I  +  G)C. 

If  the  gain  is  low  (as  it  is.  unfortunately,  for  our  system 
on  ACO).  one  can  diagonalize  the  evolution  matrix 
(36) 


v  — *  v - y(2p  +  /+  1).  (35) 

Equation  (35)  means  that  the  gain  curves  of  the  modes 
are  shifted  with  respect  to  each  other.  This  effect  has 
been  calculated  for  the  fundamental  TEM00  mode  in 
the  energy  loss  approximation  [14],  and  has  recently 
been  observed  experimentally  at  Orsay  [21].  It  should 
be  noted  that  for  many  practical  situations  where  the 
cavity  is  optimized  for  gain  on  the  TEM0(,  mode,  this 
expression  is  valid  for  only  the  lowest-order  mode.  The 
higher  modes  become  distorted  in  form  as  well  as 
simply  shifted  in  resonance  parameter  by  (35). 

2.2.  The  Low-Field  Oscillator 

We  now  discuss  some  consequence  of  the  linearity  of 
the  low-field  problem  on  the  optimization  of  an  optical 
cavity  for  an  FEL  oscillator  experiment.  In  such  an 
experiment  the  light  pulses  reflect  n  times  on  the  cavity 
mirrors  (n^2)  between  interactions  with  electrons  in 
the  undulator.  The  matrix  governing  the  mode  evolu¬ 
tion  from  one  amplification  to  the  next  is 

( 1  +  G)C .  (36) 

where  G  in  the  gain  matrix  defined  previously,  and  C  is 
the  cavity  matrix  describing  the  n  reflections  on  the 
mirrors.  In  a  set  of  cavity  eigenmodes.  C  is  diagonal 

C.x'  =  .  (37) 

with  c)y  =  /(>  =  g/exp(ia/),  where  1—  q]  are  the  total 
losses  on  the  n  reflections,  including  transmission, 
absorption,  scattering,  and  diffraction.  If  diffrac¬ 
tion  is  negligible,  the  eigenvectors  v1  become  the 


(/  +  G)C.‘ = 


(40) 


using  the  cavity  eigenmodes  V  as  the  basis  for  a 
perturbation  expansion  of  the  new  modes  To  first 
order  in  the  non-degenerate  case,  the  result  is 


G  =  (?>ei,'(l+^). 


/)  t! 

-'=*>  +  V  E-t-  -  n  \m . 
,,  e11'  -  o  e"-  '  m> 


(41) 

(42) 


Under  these  conditions,  the  FEL  design  is  optimized 
by  maximizing  the  diagonal  term  corresponding  to 
the  desired  mode.  From  (33)  and  (12)  it  is  clear  that  the 
beam  size  »v0  must  be  reduced  dow  n  to  the  order  of  the 
electron-beam  size  in  order  to  optimize  the  coupling, 
but  if  the  mode  becomes  too  divergent,  the  time 
dependent  terms  in  Em  and  E„  of  (28)  begin  to  reduce 
the  gain.  The  optimal  situation  lies  between  these  two 
extremes,  and  has  been  calculated  in  detail  (using  the 
energy  loss  approximation)  by  Colson  and  Elleaume 
[14]. 

The  optimization  procedure  must  also  be  limited  by 
the  stability  condition  [17]  on  the  cavity.  For  the 
Orsay  experiment,  the  radius  of  curvature  chosen  to 
optimize  the  small-signal  gain  was  R  =  3  m,  which  is 
acceptably  close  to  the  stability  limit  of  2.75  m.  There 
are  cavity  designs  in  which  C  is  degenerate  for  which 
the  optimization  procedure  is  not  necessary.  For  these 
designs.  x;  is  a  constant  independent  of  the  index.  In 
these  cavities,  any  combination  of  modes  reproduces 
itself  after  n  reflections.  If  n  =  2  as  in  the  Stanford  and 
the  Orsay  experiments,  the  concentric  and  the  plane- 
parallel  cavities  are  degenerate,  and  the  confocal  cavity 
is  degenerate  on  the  p  modes  iquasi-degenerate).  These 
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Fig.  I.  Simplified  schematic  diagram  of  Ihe  gam-measurement  ap¬ 
paratus  [23]  showing  the  argon-laser  focussing  system,  the  collimat¬ 
ing  iris,  and  the  double  demodulation  detection  system 


Fig.  2.  The  measured  gam  as  a  function  of  the  iris  diameter  [16] 
normalized  to  the  measured  beam  waist  at  the  iris.  The  solid  points 
were  taken  closing  the  iris  and  the  open  points  while  opening  it.  The 
error  bars  are  the  one  sigma  statistical  errors.  All  points  have  the 
same  horizontal  error  bar  which  is  shown  for  the  point  at  2.7.  The 
solid  curve  is  calculated  using  the  measured  values  for  the  electron 
and  laser  beam  sizes.  The  effect  of  each  higher-order  mode  is  shown 
by  the  dashed  curves 


cavity  designs,  however,  are  useless  since  they  stand 
critically  on  the  stability  boundary.  The  tolerance  on 
the  mirror  radius  of  curvature  is  on  the  order  of  \g00\ 
(obtained  from  the  perturbation  expansion)  which  is 
difficult  to  meet  if  the  gain  is  low.  For  two-mirror 
devices  where  n>2,  such  as  the  Novosibirsk  experi¬ 
ment  w  here  n  =  8,  in  general  for  mirrors  of  equal  radius 
of  curvature  there  exist  ni 2+  1  cavity  designs  with  a 
degenerate  C  matrix: 


2  tan  — 

it 

and  n  2  with  a  quasi-degenerate  C  matrix  in  which  the 
odd  /  modes  change  sign  on  every  amplification.  Only- 
two  of  the  degenerate  and  one  of  the  quasi-degenerate 
cavities  correspond  to  the  unstable  cavities:  the  others 
are  potentially  useable  in  an  experiment.  The  value  of  a 
degenerate  C  matrix  is  that  the  eigenvectors  of  the 
amplifier  plus  cavity  matrix  (36)  arc  equal  to  those  of 
the  gain  matrix  alone,  multiplied  by  a  constant.  This 
degeneracy  allows  the  cavity  to  oscillate  on  the  most 
favorable  combination  of  modes  which  best  fits  the 
electron  beam  shape.  In  this  manner,  the  gain  can  be 
increased  by  factors  of  two  or  three  over  the  gain  of  an 
optimized  TEM„„  mode,  particularly  if  the  electron 
beam  size  is  smaller  than  the  TtMn„  mode.  The 
tolerance  on  the  mirror  radius  for  the  degeneracy  of  C 


will  still  be  tight,  and  the  experimental  utility  of  these 
cavities  remains  to  be  investigated. 

3.  Application  to  Gain-vs-Aperture  Experiment 

3.1.  Description  of  the  Experiment 

The  gain  of  the  Orsay  FEL  has  recently  been  measured 
with  the  optical  klystron  in  place  [22]  in  an  amplifier 
experiment  using  an  external  argon  ion  laser  to  pro¬ 
vide  the  coherent  mode.  A  detailed  description  of  the 
apparatus  can  be  found  in  [23],  and  a  schematic 
description  is  given  in  Fig.  1.  The  laser  beam  is  anal¬ 
yzed  at  a  distance  d  from  the  optical  klystron  after 
passing  through  an  adjustable  collimating  iris  (Fig.  1), 
which  is  centered  on  the  laser  mode  emerging  from  the 
interaction  region.  The  gain  is  measured  as  the  ratio  of 
the  power  detected  in  phase  with  both  the  electron 
repetition  frequency  and  the  chopper  frequency  (the 
amplified  power)  divided  by  the  power  in  phase  with 
the  chopper  alone  (the  incident  laser  power). 
Calibration  is  performed  as  in  [23]. 

The  again  is  recorded  as  a  function  of  the  iris  aperture, 
and  large  variations  are  observed  [16].  One  set  of  data 
points  is  reproduced  in  Fig.  2.  where  the  gain  is 
normalized  to  its  value  for  the  iris  completely  open, 
and  ihe  iris  diameter  is  normalized  to  the  measured 
beam  waist  at  the  iris.  The  data  is  taken  at  maximum 
gain,  which  means  vf%0  for  the  optical  klystron,  and 
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f-ig  .V  Calculated  curves  for  the  gain  as  a  function  of  iris  diameter 
under  the  conditions  of  the  Orsay  experiment  [!(>].  The  electron 
beam  dimension  i'  =  rr|  rr  /. L  is  varied  to  show  the  effects  of  the 
beam  size  on  the  excitation  of  the  higher  order  modes.  The  value 
i  =  0  ’6  corresponds  to  rr  =  (I  TS  mm  which  is  very  close  to  the  value 
at  which  the  experimental  points  were  recorded 


rig.  4.  Calculated  gain  as  a  function  of  iris  diameter  for  several  iris 
positions  cl.  under  the  conditions  of  the  Orsay  experiment  [16].  The 
ratio  of  the  iris  to  optical  klystron  distance  if  divided  by  the  optical 
klystron  length  L  is  varied  through  the  range  0.5  to  x. .  The 
experimental  points  of  Fig.  2  were  taken  for  if  L  =  9 


the  laser  beam  was  carefully  aligned  to  within  about 
0.05  mm  of  the  axis  of  the  electron  beam.  The  change  in 
the  measured  gain  as  the  iris  is  closed  means  that  the 
laser  is  not  uniformly  amplified  in  its  transverse  profile. 
In  fact,  this  experiment  provides  a  very  sensitive 
technique  for  measuring  the  power  emitted  into  the 
higher-order  modes  even  in  the  small  gain  limit  and  for 
a  monomode  input  beam.  Clearly  a  calculation  of  the 
gm„  is  necessary  in  order  to  explain  these  results.  In  the 
next  section,  we  apply  the  theory  we  have  developed  to 
the  case  at  hand,  and  in  Sect.  3.3,  precise  comparison  is 
made  between  the  experimental  and  the  theoretical 
results. 


3.2.  Multimode  Emission  in  a  Single-Mode 
Amplifier  Experiment 

In  this  subsection,  we  assume  the  incident  wave  is  a 
single  mode  TEM00  beam  with  a  weak  field  (|«0|  <  1), 
and  perfectly  aligned  onto  the  electron  beam.  As 
discussed  previously,  we  take  the  cylindrical  eigen- 
modes  based  on  the  form  of  the  input  beam.  Using 
the  notation  of  Sect.  1,  the  input  laser  field  reads 

£f(  r )  =  t-0  rfe”'  ‘’,r’ .  (44) 

where  the  subscript  0  refers  to  the  TFM„0  mode  of  (12) 


and  (13).  From  (29).  the  output  field  £s(r)  becomes 
E5  =  E'  +  c0  t  </,o£y<r>eiVj,r’.  (45) 

J  =  0 

Assuming  low  gain,  the  output  power  passing  through 
the  iris  aperture  is 

8^P  =  J(AS|£'|2  J dSE20  +  2c i  Re 

jl  gjoidSEoEjC*™- '■*'},  (46) 

where  Ji/S  covers  the  iris  aperture.  The  gain  is 
therefore 


2  Re 1 1  0,0 |</S£0£, 


f  dSE* 

For  purely  cylindrical  1  =  0  modes,  G  can  be  written 

x  .,  iLj  £„(.v)c~xc/.v 

G  =  2 ReJc/oo  +  X  </P„e'-P,‘n  =<>% -  .  (48) 

fe'(/x 

0 

where  is  the  Rayleigh  range  of  the  laser  mode.  :  is 
the  distance  between  the  iris  and  the  laser  beam  waist, 
and  A'  =r"/2w2(r)  where  r0  is  the  iris  diameter.  There 
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are  two  interesting  limiting  eases: 

G  =  2Re{</u0}  .V— x  (iris  open).  (49) 

G  =  2  Re |(t/00)+  t  </l,,)ci'"'3n  '-J 

JV  — ►  0 (iris  closed).  (50) 

It  is  obvious  from  (48  50)  that  the  gain  changes  with 
iris  diameter  in  a  way  which  depends  on  the  magni¬ 
tudes  of  the  off-diagonal  terms  in  the  gain  matrix. 

The  generalization  is  straightforward  to  the  case  of  the 
multimode  input  beam,  and  to  imperfect  alignment  of 
the  laser  and  electron  beams,  although  the  calculation 
becomes  more  difficult.  This  calculation  also  applies  to 
high-power  input  laser  beams  (u0^>l)  provided  one 
keeps  in  mind  that  the  cip0  are  functions  of  a0. 

3.3.  Application  to  the  Orsay  Experiment 

The  experimental  points  shown  on  Fig.  2  were  taken 
under  the  following  approximate  conditions 

laser 

beam:  measured  beam  waist  u0  =0.67 mm 

wavelength  =  5145  A 

-  measured  beam  waist  at  iris  u(r)  =  2.7  mm 

-  distance  from  optical  klystron  to  iris 
d  =  1 1.6m. 

electron 

beam:  -  Gaussian  and  cylindrical  with 

<t^0.32  mm. 

optical 

klystron:  -  Nd  =  80  [19.  22.  24] 

-  resonance  parameter  corresponding  to 
maximum  gain  with  iris  open. 

The  solid  curve  of  Fig.  2  has  been  calculated  using  (27, 
28.  and  50)  for  the  planar  configuration  (24  and  25), 
taking  into  account  the  10  lowest  order  /  =  0  modes. 
The  dashed  curves  of  Fig.  2  show  the  contribution  of 
each  individual  mode.  These  curves  are  the  same 
whether  an  undulator  or  optical  klystron  is  used.  Very 
similar  curves  (not  shown)  were  calculated  for  other 
resonance  parameters  indicating  that  as  expected,  the 
diffraction  effects  do  not  change  much  as  a  function  of 
detuning  parameter  for  modes  with  low  divergence. 
Figure  3  shows  the  calculated  effect  for  several  dimen¬ 
sionless  electron  beam  transverse  dimensions 
Z=  |  an  /L  =0.4.0.76.  1.5,  3.  where  L  is  the  length  of 
the  magnetic  interaction  region.  27  =  0.76  corresponds 
to  the  value  a  =0.35  mm.  close  to  that  used  in  Fig.  2. 
The  flattening  of  the  curves  as  I  is  increased  to  27  =  3  is 
due  to  the  vanishing  of  </0>  </00  for  j  +  0  as  <r 
increases. 


Figure  4  shows  the  calculated  effect  for  various  iris 
distances  d  from  the  optical  klystron  center,  normal¬ 
ized  to  the  optical  klystron  length:  d  L  =  0.5,  1.  1.8.  2.5. 
and  x.  The  experimental  points  of  Fig.  2  were  obtain¬ 
ed  for  d  L  =  9.  The  inversion  of  the  effect  is  due 
primarily  to  the  term  exp[i2ptan“  '(;,;„)]  with  p=  1  in 
(50)  which  switches  from  +  1  to  -  1  as  r  goes  from  zero 
to  infinity  (the  mode  p=l  gives  the  predominant 
effect).  At  short  distances  the  TKMln  mode  interferes 
constructively  on  axis  I.  and  at  long  distances,  it 
changes  sign. 
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A  storage-ring  free-electron  laser  oscillator  has  been  operated  above  threshold  at  a 
visible  wavelength  A  =  6500  A. 

PACS  numbers:  42.60.-v 


A  free-electron  laser  (FEL)  generates  stimu¬ 
lated  radiation  by  interaction  of  relativistic  free 
electrons  with  a  spatially  oscillating  magnetic 
field  (undulator).  The  wavelength,  a,  of  the  stim¬ 
ulated  radiation  emitted  along  the  axis  is  given 
approximately  by  (in  SI  units) 


(1) 


where  e  and  m  are  the  electron  charge  and  mass, 
c  is  the  speed  of  light  in  vacuum,  B0  and  A0  are 
the  undulator  peak  magnetic  field  and  period,  and 
y  is  the  total  electron  energy  divided  by  me2.  The 
first  operation1  of  a  FEL  has  been  demonstrated 
by  the  Stanford  group  at  A  =3.4  pm  using  43-MeV 
electrons  produced  by  a  linear  accelerator.  In 
any  practical  design,  A0-e  3  cm;  therefore  it  is 
obvious  from  Eq.  (1)  that  operation  of  a  visible  or 
ultraviolet  FEL  needs  a  high-energy  accelerator 
of  ymc2z  100  MeV.  In  this  energy  range,  storage 
rings  are  superior  to  linear  accelerators  in  terms 
of  electron  density,  energy  spread,  and  emit- 
tance,  thus  resulting  in  a  higher  optical  gain. 

In  this  Letter,  we  report  the  first  operation  of  a 
storage-ring  free-electron  laser  oscillator  above 
threshold  at  A  =*  6500  A. 

We  have  used  the  Orsay  storage  ring  ACO,  an 


“older”  machine,  first  operated  in  1965.  Its 
characteristics  for  this  experiment  are  given  in 
Table  I.  A  permanent  magnetic  undulator  of  17 
periods  optimized  for  electrons  of  240-MeV  ener¬ 
gy  has  been  built.2  However,  the  limited  straight- 
section  length  (1.3  m)  available  on  ACO  left  us 
with  a  low  gain  of  a  few  parts  in  10* 4  per  pass  in 
earlier  work.3  Laser  oscillation  was  therefore 
impossible.  Now,  a  gain  enhancement  by  a  fac¬ 
tor  2  to  7  has  been  achieved  by  modifying  the  un¬ 
dulator  into  an  optical  klystron4  (OK).  This  was 
done  by  replacing  the  three  central  periods  by  a 
three-pole  wiggler  (dispersive  section).  This 
section  strongly  enhances  the  bunching  originat¬ 
ing  from  electron  and  radiation  interaction  in  the 
first  undulator  section  and  therefore  gives  a  larg¬ 
er  energy  exchange  between  electrons  and  radia¬ 
tion  field  in  the  second  undulator  section.4 

Mirror  reflectivity  degradation  has  been  ob¬ 
served3  at  240-MeV  electron  energy.  This  de¬ 
gradation  is  mainly  due  to  the  uv  part  of  the  spon¬ 
taneous  emission  from  the  klystron.  This  prob¬ 
lem  was  overcome  by  operating  the  OK  at  mini¬ 
mum  K  (by  increasing  the  magnetic  gap)  to  mini¬ 
mize  the  harmonic  content  of  the  spontaneous 
emission,  and  the  storage  ring  at  the  minimum 
electron  energy  to  also  diminish  the  flux.  It  is 


TABLE  I.  ACO  characteristics  in  the  FEL  experiment. 


Energy 

160  -  166  MeV 

Circumference 

22  m 

Bunch  to  bunch  distance 

11  m 

Electron  beam  current  for  oscillation 

16  to  100  mA 

rms  bunch  length  o, 

0.5  to  1  ns 

rms  bunch  transverse  dimensions  o,,  oy 

0.3  to  0.5  mm 

rms  angular  spread  »/,  ov' 

0.1  to  0.2  mrad 

rms  relative  energy  spread 

(0.9  to  1.3)  >  10‘3 

Electron  beam  lifetime 

60  to  90  min 

1652 


tv)  1983  The  American  Physical  Society 


■Vlv^rlvlv: 


'  '-A'*-* 


\'<>l  I  Ml  SI,  Nl  M  HI  K  IK 


PHVSICAf  REVIEW  LETTERS 


31  ()<  iohik  198.3 


TABLE  II.  Optical  cavity  characteristics. 


Length 

5.5  m 

Mirror  radius  of  curvature 

3  m 

Rayleigh  range 

1  m 

Wavelength  of  maximum  Q 

620  to  680  nm 

Average  mirror  reflectivity  at  6328  A 

99.965T 

Round  trip  cavity  losses  at  6328  A 

7  x  10‘4 

Mirror  transmission 

3  x  10'5 

possible  by  decreasing  both  K  and  y  to  keep  a 
constant  around  6500  A  which  gives  the  maximum 
reflectivity  for  the  mirrors  used  in  this  experi¬ 
ment.  The  OK  was  finally  operated  at  A'  =1.1  to 
1.2  with  a  dispersive  section  characterized  by5 
Nd  -  95  where  Sd  is  the  number  of  laser  optical 
wavelengths  passing  over  an  electron  in  the  dis¬ 
persive  section.  The  optical  cavity  characteris¬ 
tics  are  given  in  Table  II. 

Laser  oscillation  was  obtained  after  a  careful 
alignment  of  the  electron  beam  on  the  cavity  axis 
(within  0.1  mm  all  along  the  1.3-m  OK  length) 
and  maximization  of  emission  as  a  function  of  the 
storage-ring  radio-frequency  cavity  voltage  and 
the  optical-cavity  length.  Laser  operation  re¬ 
quires  very  precise  synchronism  between  light- 
pulse  reflections  and  electron-bunch  revolution 


FIG.  1.  Normalized  laser  power  as  a  function  of  rf 
frequency  and  equivalent  mirror  displacement.  Curve 
a  U  recorded  close  to  the  maximum  gain/loss  ratio 
•od  curve  b  close  to  the  laser  threshold.  The  shift 
between  the  two  curves  is  probably  due  to  some  slight 
cavity-length  drift. 


frequency.  To  avoid  backlash  and  mirror  mis¬ 
alignment,  fine  tuning  of  the  cavity  length  was 
performed  by  slightly  changing  the  frequency  in¬ 
stead  of  translating  the  mirrors.  Laser  oscilla¬ 
tion  lasted  typically  1  h  after  each  electron  injec¬ 
tion. 

Figure  1  shows  two  “detuning  curves”  of  laser 
power  (normalized  to  the  maximum)  versus  fre¬ 
quency  variation  and  equivalent  mirror  displace¬ 
ment.  Curve  (a)  has  a  3.4  pm  full  width  at  half 
maximum  of  equivalent  mirror  displacement; 
curve  ( b ),  as  recorded  much  closer  to  laser 
threshold,  has  only  a  1.6  pm  full  width  at  half 
maximum.  The  shift  in  displacement  between 
curves  (a)  and  (b)  is  probably  due  to  some  slow 
cavity  length  drift  [  e.g.,  a  temperature  drift  of 
(0.02  °C)/(30  min)  would  be  sufficient].  In  this  ex¬ 
periment,  the  gain  was  not  proportional  to  the 
electron  current,  mainly  because  of  the  anoma¬ 
lous  bunch  lengthening  and  energy  spreading  at 
high  current;  these  effects  make  the  gain  versus 
ring  current  reach  a  maximum  and  then  decrease. 
One  consequence  of  this  is  that  the  ratio  of  gain 
to  cavity  losses  always  remained  just  above  1 
during  laser  operation.  Wider  detuning  curves 
are  expected  for  higher  gain/loss  ratios. 

Figure  2  presents  the  horizontal  and  vertical 
transverse  profiles  of  the  laser  output.  They  are 
in  a  very  good  agreement  with  the  expected  cavi¬ 
ty  TEMqo  profile.  The  slight  discrepancies  might 
arise  from  some  residual  instability  of  the  laser 
too  close  to  threshold  or  nonuniform  mirror  re¬ 
flectivity. 

Figure  3  shows  the  laser  time  structure  in  a 
200-ms  total  time  scale.  Quasirandom  peaks  ap¬ 
pear  in  curve  («).  Curve  (b)  is  obtained  by  in¬ 
creasing  the  detector  sensitivity  and  shows  some 
substructure  in  the  peaks.  Although  not  seen  in 
Fig.  3,  each  subpeak  is  also  usually  more  or  less 
modulated  at  13  kHz,  very  close  to  the  theoreti¬ 
cal  rf  synchrotron  frequency,  and  also  shows  a 
typical  rise  time  around  200  ps  corresponding  to 
a  gain  minus  cavity  losses  of  2x  10"* 'pass.  It 
should  be  noted  that  this  time  structure  highly  de- 
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FIG.  2.  Experimental  laser  horizontal  {solid  curve), 
vertical  (dashed  curve)  transverse  profiles,  and  the 
calculated  (dotted  curve)  cavity  TEM00  profile. 


pends  on  how  far  the  laser  is  from  threshold. 
Quasiperiodic  peaks  have  sometimes  been  seen 
near  a  40-Hz  frequency.  And  the  laser  variation 
also  reproduces  the  pulsed  27.2-MHz  structure 
of  the  electron  beam. 

Figure  4  presents  two  spectra:  (a)  is  recorded 
without  amplification  (optical  cavity  completely 
detuned)  and  (b)  is  recorded  at  laser  operation 
(cavity  tuned).  For  case  (b),  the  laser  oscillates 
at  three  wavelengths,  with  the  strongly  dominant 
one  being  at  A  =6476  A;  each  wavelength  is  locat¬ 
ed  at  a  maximum  of  the  gain  versus  wavelength 
curve.5  Smaller  gain/loss  ratios  would  restrict 
the  number  of  laser  wavelengths  to  two  or  one. 
Typical  laser  lines  are  Gaussian  if  averaged  over 
a  long  time  scale  of  ^  1  sec,  with  2  to  4  A  full  width 
at  half  maximum.  Figure  4  also  shows  an  enlarge¬ 
ment  of  the  main  laser-line  spectrum  recorded 
by  using  a  one-dimensional  charge-coupled-de¬ 
vice  (CCD)  detector  instead  of  the  usual  mono¬ 
chromator  exit  slit.  The  aperture  time  is  3  ms. 
Each  narrow  square  peak  in  this  curve  is  record¬ 
ed  on  only  one  CCD  element  and  corresponds  to 
a  0.3  A  spectral  width  which  is  of  the  same  order 
as  the  monochromator  resolution.  We  conclude 
that  there  is  a  residual  inhomogeneous  contribu¬ 
tion  to  the  laser  linewidth  probably  connected 
with  the  long-time-scale  laser-pulse  structure 
(see  Fig.  4).  The  central  wavelength  of  any  line 
is  always  equal  to  the  wavelength  of  maximum 
emission  of  spontaneous  emission  with  the  cavity 
completely  detuned  (no  amplification)  plus  0.15 


M 


FIG.  3.  Laser  time  structure  over  a  200-ms  interval. 
Curve  b  is  recorded  in  the  same  condition  as  a,  but 
with  a  higher  detector  sensitivity. 

of  the  wavelength  interfringe  distance  (see  Fig.  4) 
instead  of  0.25  as  predicted  from  Madey’s  theo¬ 
rem.0  This  discrepancy  is  probably  due  to  the 
transverse  multimode  content  of  the  spontaneous 
emission  stored  in  the  cavity;  laser  operation  is 
only  achieved  on  the  TEMm  mode.  Laser  tunabil- 
ity  was  obtained  between  640  and  655  nm  by 
changing  the  magnetic  gap  [equivalent  to  a  change 
of  K  in  Eq.  (1)1.  The  range  of  tunability  is  in  fact 
limited  for  the  moment  by  the  mirror  reflectivity. 

A  typical  75-pW  average  output  power  has  been 
recorded  at  50-mA  current  of  166-MeV  electrons. 
This  corresponds  to  a  typical  60 -mW  output  peak 


FIG.  4.  Spectra  of  the  cavity  output  radiation  under 
two  conditions :  curves,  cavity  detuned  (no  amplifica¬ 
tion)  and  curve  b,  cavity  tuned  (laser  on). 
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power  over  the  1-ns  electron  bunch  length  (in¬ 
cluding  effects  of  both  the  long-  and  short-time¬ 
scale  pulse  structure  of  the  laser),  and  a  2-kW 
intracavity  peak  power.  75  pW  give  a  2.4*  10' 5 
efficiency  with  respect  to  the  3.1  W  of  total  syn¬ 
chrotron  radiation  power  generated  in  the  whole 
storage  ring  at  this  energy  and  current.  This  ef¬ 
ficiency  is  about  0.4  the  Renieri  typical  efficien¬ 
cy7 

V  =(j£/£)77  ml[  , 

where 

mirror  transmission 

^m'r  round  trip  cavity  losses 

and  <se/E  =1.2x10' 3  is  the  RMS  relative  energy 
spread.  In  this  experiment  t)ml,  =0.043.  This 
low  mirror  efficiency  is  mainly  due  to  the  high 
absorption  in  the  mirror  dielectric  compared  to 
its  transmission.  The  rather  low  Renieri  aver¬ 
age  output  power  is  therefore  the  result  of  high 
sensitivity  of  the  OK  to  energy  spread,  poor  mir¬ 
ror  efficiency,  and  weak  total  synchrotron  radia¬ 
tion  power  at  166  MeV  (power  is  proportional  to 
the  fourth  power  of  the  electron  energy).  The 
laser  output  power  always  decreases  with  the 
electron  beam  current  and,  except  in  regions  too 
close  to  threshold,  is  usually  proportional  to  the 
current  as  predicted  by  the  same  model.  How¬ 
ever,  important  discrepancies  appear  in  the 
structure  for  long  time  scales  (see  Fig.  3)  which 
is  predicted  to  be  constant.  Moreover,  bunch 
shortening  instead  of  bunch  lengthening  has  been 
seen  at  laser  turnon;  such  an  effect  has  already 
been  seen  at  much  lower  current  by  simulating 
the  FEL  with  an  argon  laser.8 

In  summary,  this  work  demonstrates  for  the 
first  time  the  feasibility  of  the  storage-ring  free- 
electron  laser.  Future  experiments  will  continue 
to  analyze  the  saturation  mechanism,  and  the 
long-  and  short-time-scale  structure.  Higher 
gain/loss  ratios  are  expected  from  the  use  of 
better  mirrors,  a  second  rf  cavity,  and  smaller 
beam  transverse  dimensions  obtained  by  using 
positrons  instead  of  electrons. 
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ABSTRACT 

We  report  the  first  measurement  of  the  off-diagonal  terms  in 
the  transverse  gain  matrix  in  a  free  electron  laser.  We  show  that 
the  higher  order  transverse  modes  stimulated  in  the  FEL 
interaction  diminish  in  amplitude  as  the  electron  beam  size  is 
increased,  that  the  far  field  and  the  near  field  effects  of  the 
excited  modes  are  not  equivalent,  and  that  in  this  low  divergence 
case,  the  mode  mixing  is  insensitive  to  the  resonance  parameter. 
Remarkably  close  agreement  is  demonstrated  with  the  theory  [1]. 

The  effects  of  a  multiple  mode  input  beam  are  also  observed,  and 
the  theoretical  expression  is  derived. 


X.  INTRODUCTION 


Beginning  in  1980  a  series  of  amplifier  measurements  have 
been  performed  of  the  stimulated  emission  of  the  Orsay  storage 
ring  free  electron  laser.  These  measurements  were  part  of  the 
successful  effort  to  achieve  oscillation  using  the  ACO  storage 
ring  system  [2,3].  The  original  gain  experiments  [4,5]  on  the 
superconducting  undulator,  and  a  later  series  of  measurements 
[6,7]  on  the  permanent  magnet  undulator  NOEL  and  its  optical 
klystron  configuration  demonstrated  that  the  amplitude  and  the 
spectrum  of  the  gain  could  be  correctly  calculated  by  the 
classical  independent  particle  theory  [8].  These  experiments  were 
sensitive  to  the  total  energy  emitted  by  the  stimulated  emission 
process,  and  did  not  provide  any  information  on  the  spatial 
structure  of  the  amplification  process.  In  this  paper  we  report  on 
an  extension  of  these  measurements  which  probes  the  three 
dimensional  nature  of  the  interaction. 

By  means  of  an  adjustable  aperture  placed  in  the  amplified 
beam,  we  have  been  able  to  measure  the  excitation  of  higher  order 


transverse  modes  induced  by  the  nonuniform  transverse  profile  of 
the  electron  beam.  This  information  is  mapped  into  the  spatial 
distribution  of  the  electric  field  at  the  analyzing  aperture.  The 
ratio  of  the  transmitted  stimulated  power  to  the  transmitted  power 
of  the  external  laser  is  measured  as  a  function  of  the  aperture 
diameter.  The  functional  dependence  of  the  ratio  on  the  aperture 
diameter  is  a  diagnostic  on  the  transverse  mode  mixing  which  has 
occurred  during  the  amplification  process. 

If  the  spatial  distribution  of  the  stimulated  radiation  is 
identical  to  that  of  the  incident  laser,  their  diffraction 
characteristics  are  identical,  and  the  ratio  of  the  powers 
transmitted  through  the  aperture  will  be  independent  of  its 
diameter.  A  dependence  will  be  seen  only  if  the  stimulated 
radiation  contains  a  different  combination  of  transverse  modes 
than  the  incident  radiation.  The  aperture  dependence  of  this  ratio 
is  therefore  a  good  diagnostic  of  mode  mixing  in  the  amplifier. 

In  a  free  electron  laser,  the  input  laser  mode  is  perfectly 
replicated  only  under  the  conditions  that  the  electron  beam  is 
uniform  across  the  optical  mode,  the  electric  field  strength 
remains  below  the  saturated  level,  and  the  divergence  of  the  laser 
is  negligible  in  the  interaction  region.  If  the  electron  beam 


dimensions  are  comparable  to  those  of  the  optical  mode,  if 
saturation  occurs,  or  if  the  divergence  is  significant,  higher 
order  modes  will  be  excited.  In  most  configurations  of  practical 
interest,  one  or  more  of  the  latter  conditions  will  hold. 

The  mode  content  of  the  optical  field  propagating  through  the 
FEL  can  be  conveniently  thought  of  as  a  vector  whose  elements  are 
the  amplitudes  and  phases  of  the  modes.  As  the  beam  propagates 
through  free  space,  the  phases  evolve  in  a  fashion  characteristic 
of  each  mode,  but  the  amplitudes  remain  constant  (in  the  absence 
of  losses).  The  amplifier  contributes  an  additional  vector  of  complex 
field  amplitudes  which  describes  the  stimulated  radiation,  and 
which  is  in  general  not  parallel  to  the  input  vector  [1],  The 
field  vector  at  the  output  of  the  FEL  amplifier  is  the  sum  of  the 
input  vector  propagated  to  the  output  and  the  stimulated  emission 
vector.  The  magnitude  squared  of  the  resultant  will  be  larger  than 
that  of  the  input  by  a  factor  which  is  the  net  amplification,  and 
it  will  be  rotated  in  the  multi-dimensional  mode  space.  If  the 
gain  is  large  and  mode  mixing  occurs,  the  output  field  vector  will 
be  strongly  rotated.  Under  these  circumstances,  the  laser  can  be 
expected  to  run  multi-mode,  or  even  to  couple  the  modes  so  strongly 
that  the  modal  decomposition  becomes  useless  as  in  the  guided  wave 
case.  A  self-consistent  calculation  of  the  high  gain  field  dist¬ 
ribution  must  be  performed  to  identify  the  steady  state 
field  vector  and  the  gain  it  experiences.  In  a  low  gain 


oscillator,  however,  the  same  level  of  mode  mixing  will  produce  no 
effects  because  the  rotation  of  the  field  vector  by  the  gain  is  negligible. 


The  gain  of  the  Or say  FEL  is  very  low,  on  the  order  of  10 
[4-7].  Even  though  the  conditions  for  mode  mixing  are  present  (the 
dimensions  of  laser  and  electron  beam  are  comparable),  the 
rotation  of  the  field  vector  is  minuscule.  As  could  be  expected 
under  these  conditions,  the  oscillator  produced  a  nearly  pure 
TEMoo  Gaussian  mode  [3],  (Even  the  Stanford  oscillator  which  has  a 
gain  of  about  10%  [9]  produces  a  Gaussian  mode  [10]).  The  newer 
FEL  systems  now  under  construction  at  Stanford,  Orsay,  and 
elsewhere  are  projected  to  have  gain  on  the  order  of  unity.  The 
mixing  phenomena  of  the  transverse  modes  will  play  an  important 
role  in  the  operation  of  these  new  devices.  Our  measurement 
technique  permits  us  to  obtain  some  information  on  the  three 
dimensional  amplification  problem  even  in  a  devicewith  such  a  low 
gain  that  no  effects  can  be  observed  in  the  modes  of  the  oscillator. 

Two  series  of  experiments  have  been  performed,  both  with  the 
undulator  NOEL  modified  into  the  optical  klystron  configuration 
[7],  The  first  series  of  experiments  (which  we  will  identify  as 
series  I)  served  to  identify  the  phenomenon,  and  stimulated  our 
theoretical  effort  [l]  on  the  problem.  A  second  series  (series  II) 
was  planned  in  order  to  extend  the  comparison  between  theory  and 
experiment  by  changing  the  placement  of  the  analyzing  aperture.  We 
report  here  on  the  results  of  these  two  experiments. 
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7 r.  MODE  mixing  theory  applied  to  the  gain -vs.- aperture  experiment 


Elleaume  and  Deacon  [l]  have  identified  the  equations  of 
motion  for  the  excitation  of  the  higher  order  transverse  modes  in 
an  undulator  and  solved  them  in  the  low  field  case.  Their 
principal  low  field  results  are  that  the  complex  mode  amplitudes 
evolve  according  to 

c4  = 


The  linear  relationship  between  input  and  output  mode  amplitudes 
holds  even  in  the  high  gain  region  [l],  although  the 
transformation  becomes  difficult  to  calculate.  In  what  follows,  we 
restrict  our  examples  to  the  small  gain  limit.  In  this  limit,  for 
low  divergence,  and  for  small  beam  slippage  parameter  (as  in  the 
Orsay  experiment),  the  elements  of  the  gain  matrix  £  are 


E.J M)  e 


We  note  that  this  result  does  not  include  the  effects  of  the 
divergence  of  the  optical  modes  or  the  electron  beam,  a  valid 
approximation  in  the  case  of  the  Or say  experiment,  where  the 
Rayleigh  length  is  1.7  times  the  length  of  the  optical  klystron, 
and  the  beta  functions  are  1.5  to  3  times  longer. 

In  the  large  divergence  case,  the  time  integrals  in  equation 
(30)  of  reference  [1]  depend  on  both  the  mode  numbers  and  the 
transverse  position.  The  spectral  term  between  the  curly  brackets 
of  (2)  which  is  generated  by  this  time  integral  is  therefore 
intimately  tied  to  the  transverse  integral  and  cannot  be  factored 
out  of  the  gain  matrix.  Under  these  conditions,  the  mode  mixing 
phenomena  may  have  a  strong  dependence  on  the  resonance  parameter 
A  full  treatment  of  the  problem  including  divergence  would 
probably  involve  a  numerical  solution. 

We  can  apply  the  low  divergence  results  to  the  case  of  the 
optical  klystron  (or  the  tapered  wiggler)  since  there  is  no 
essential  change  in  the  physics  which  goes  into  the  derivation  of 
(2).  Within  the  same  approximations,  the  only  change  is  in  the 
longitudinal  dimension.  The  time  integrals  which  lead  to  (2) 
become  more  complex  so  that  the  spectral  term  changes  its  form. 
Since  this  term  factors  out  in  the  low  divergence  case,  for  a 


fixed  frequency  experiment  such  as  our  amplifier  work,  the 
spectral  term  is  optimized  once  and  remains  constant  during  the 
experiments.  As  long  as  the  longitudinal  dynamics  are  separable 
mathematically,  the  mode  mixing  phenomena  will  depend  only  on  the 
beam  dimensions,  the  resonator  configuration,  and  the  gain. 

Let  us  apply  equations  (1)  and  (2),  following  [1],  to  the 
case  in  which  the  incident  beam  contains  an  arbitrary  number  of 
modes.  The  input  field  is 


where  we  use  the  summation  convention:  repeated  indices  are 
summed.  The  functions  E^(x,y),  and  i^(x,y)  are  defined  in  [1] 
as  the  normalized  amplitude  and  the  phase  of  the  mode.  The 


output  field  is,  from  (1), 


The  gain  measured  in  this  experiment,  which  is  the  ratio  of 
the  stimulated  power  to  the  incident  power  transmitted  through  the 
aperture,  is,  for  small  gain  systems, 


where  the  integral  is  over  the  surface  of  the  aperture  opening. 
For  simplicity  of  notation,  we  define  a  normalized  mode 
integral 


which  defines  the  elements  of  a  Hermitian  matrix  T_.  Since  the 
integral  is  over  the  surface  out  to  the  aperture,  it  is  clear  that 
these  quantities  depend  on  the  diameter.  Due  to  the  different  mode 
structures,  each  t^  has  its  own  aperture  dependence,  and  its  own 
unique  phase.  Substituting  into  (5),  the  gain  becomes 


This  relation  can  be  expressed  more  compactly  in  matrix 
notation.  Adopting  C  as  the  vector  of  complex  mode  amplitudes,  we 
find  the  result 


<tJT*c> 

where  the  superscript  H  means  Hermitian  con]ugate.  Physically 
speaking,  this  result  is  not  hard  to  understand.  The  numerator  is 
the  stimulated  power  term.  Each  input  mode  cj  produces  a  set  of 
stimulated  modes  via  the  gain  matrix  g;j  :  the  vector  of  stimulated 
modes  is  G|c).  Each  stimulated  mode  "beats"  with  another  mode  c^ 
when  it  passes  through  the  aperture;  this  other  mode  is  either  an 
incident  mode  or  an  amplified  mode  ^C|G,  .  Each  combination  of 

modes  produces  a  nonzero  contribution  to  the  transmitted  power 
when  the  aperture  diameter  is  finite.  This  contribution  is 
proportional  to  a  particular  mode  integral,  whose  amplitude  is 
given  by  2(jJ  •  The  measured  gain  is  the  ratio  of  the  stimulated 
power  to  the  incident  power,  and  the  denominator  takes  care  of  the 
effects  of  the  iris  on  the  incident  beam. 


To  make  the  situation  more  clear,  let  us  look  at  the  explicit 
sums  in  the  case  of  two  input  modes,  "o”  and  "p”.  We  assume  that 
the  input  beam  is  predominantly  "o”,  and  linearize  in  Cp.  Equation 
(8)  reduces  to 


where  we  have  rationalized  the  denominator,  and  collected  the 
off-diagonal  terms  in  the  sums. 


We  note  here  that  for  single  mode  input  (Cp  =  0),  this  result 
reduces  to  that  of  equation  (48)  in  [1].  The  additional  terms  are 
the  various  beating  terms  of  the  'stimulated  emission  with  the 
input  modes,  integrated  across  the  iris.  The  last  term,  for 
example,  is  the  product  of  the  mth  stimulated  mode  (produced  via 


from  the  pth  input  mode)  and  the  0th  input  mode, 


integrated 


across  the  analyzing  aperture,  and  normalized  to  the  input 


intensity.  It  is  the  last  three  sums  which  are  apparently 
responsible  for  the  difference  between  the  experimental  points  in 
figures  7  and  8  and  the  single  input  mode  theory  of  [l]. 


If  the  electron  beam  is  very  large  compared  to  the  laser 
beam,  the  off-diagonal  terms  (2)  of  G.  go  to  zero  by  the 
orthogonality  of  the  modes.  The  diagonal  terms  all  have  the  same 
magnitude  due  to  the  mode  normalization.  Equation  (9)  reduces, 
under  these  circumstances,  to 

G  =  2  Re^„  j-  (large  electron  beam)  (10) 

which  means  that  the  gain  for  a  combination  of  two  modes  is  the 
same  as  that  of  one  of  them  alone.  If  we  redefine  the  sum  of  two 
modes  as  our  new  starting  mode,  it  is  clear  that  an  arbitrary 
third  mode  amplitude  can  be  added  without  changing  the  result.  It 
follows  by  induction  that  for  an  arbitrary  combination  of  modes, 
the  direction  of  the  mode  vector  is  not  changed  by  the  FEL 
amplification  process  for  a  large  beam. 

A  measurement  of  the  gain  under  these  circumstances  will  show 
no  variation  as  a  function  of  the  transverse  position  in  the  beam 
(as  shown,  (10)  is  independent  of  iris  diameter).  Any  variation  of 
the  gain  in  the  transverse  plane  is  proof  of  mode  mixing.  Although 
the  data  of  figures  7  and  8  were  taken  with  an  unknown  multiple 
mode  input,  the  deviation  from  unity  of  the  results  are  a  clear 
indication  of  the  presence  of  mode  mixing.  The  finite  beam  size 
allows  off-diagonal  terms  in  (2),  and  thereby  produces  the  mode 
mixing  observed  in  this  experiment. 


7TT .  EXPERIMENTAL  METHOD 


The  apparatus  we  used  was  essentially  identical  to  that  reported 
in  the  original  gain  measurements  [4-7].  Those  unfamiliar  with  the 
details  of  the  setup  should  refer  to  the  schematic  diagram  in  [4] 
along  with  the  associated  description  of  the  experiment  and  the 
calibration  proceedure.  A  chopped  external  argon  laser  is 
focussed  coaxially  onto  the  electron  beam  in  the  optical  klystron 
where  it  is  amplified  by  the  FEL  interaction.  The  output  beam  is 
passed  through  an  adjustible  aperture,  filtered  in  a 
monochromator,  and  detected  in  a  resonant  detector.  The  amplified 
power  is  extracted  in  a  double  demodulation  scheme:  a  13  Mhz 
lock-in  detects  the  sum  of  the  spontaneous  and  the  stimulated 
power,  and  a  subsequent  lock-in  operating  at  the  frequency  of  the 
chopper  extracts  a  signal  proportional  to  the  stimulated  power  at 
the  detector.  The  FEL  gain  is  given  by  the  ratio  of  the  stimulated 
power  to  the  incident  laser  power. 

In  this  experiment,  the  gain  is  measured  as  a  function  of  the 
opening  of  an  iris  aligned  on  the  amplified  laser.  Both  the 
stimulated  power  and  the  incident  laser  power  are  measured  for 
several  apertures;  their  ratio  is  the  "gain"  as  a  function  of  the 
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iris  diameter.  If  both  beams  were  pure  TEMoo  modes,  the  gain 
measured  this  way  would  be  independent  of  the  iris.  If,  on  the 
other  hand,  the  stimulated  radiation  is  a  combination  of  modes, 
the  measured  ratio  will  show  some  variation  as  the  aperture  is 
varied,  since  closing  the  aperture  breaks  the  orthogonality  of  the 
modes.  A  measurement  of  the  gain  as  a  function  of  the  iris  is 
therefore  a  source  of  information  on  the  mode  content  of  the 
stimulated  radiation.  As  explained  below,  we  have  measured  the 
mode  amplitudes  stimulated  from  a  single  mode  input  beam,  which 
provides  direct  information  on  the  cross  terms  in  the  transverse 
gain  matrix. 

In  setting  up  for  the  experiment,  the  undulator  gap  is 
scanned  to  set  a  maximum  of  the  optical  klystron  gain  spectrum  [7] 
at  the  laser  wavelength.  Alignment  of  the  incident  laser  on  the 
electron  beam  is  performed  by  optimizing  the  amplitude  of  the 
stimulated  power  when  adjusting  the  angle  and  position  of  the 
electron  beam.  A  half-wave  plate  in  the  laser  input  beam  is  also 
adjusted  to  ensure  that  the  polarization  of  the  input  mode  is 
parallel  to  that  of  the  undulator  emission. 

One  analyzing  aperture  was  set  up  as  close  as  possible  to  the 
interaction  region:  just  outside  the  exit  window  at  a  distance  of 


288  cm  from  the  center  of  the  optical  klystron.  This  unit  was  a 
set  of  fixed  apertures  of  diameters  .5,  1.0,  1.5,  2.0,  and  2.5  mm, 
mounted  on  a  rotating  exchanger  block.  The  selected  aperture  could 
be  rotated  into  position  and  aligned  precisely  in  the  plane 
transverse  to  the  laser.  Before  each  measurement,  the  new  aperture 
was  aligned  by  recentering  the  diffracted  beam  downstream  of  the 
aperture  onto  the  original  axis.  A  second  aperture  was  placed  far 
from  the  interaction  region  at  a  distance  of  14.7  m.  Since  the 
laser  beam  was  larger  here,  a  minimum  aperture  of  2  mm  was 
sufficient,  and  an  adjustible  iris  was  used,  eliminating  the  need 
for  realignment  on  each  change  in  diameter.  The  powers  of  the 
transmitted  laser  and  of  the  stimulated  emission  were  measured  in 
separate  channels  and  recorded  for  about  100  seconds.  The  laser 
was  then  blocked  in  order  to  permit  a  measurement  of  the  noise 
background;  the  aperture  diameter  was  readjusted;  the  aperture  was 
realigned  (if  necessary);  and  the  two  power  measurements  were 
repeated.  The  curves  of  figures  4-8  were  taken  by  this  method. 

In  order  to  verify  that  no  detection  artifacts  were 
masquerading  as  real  effects,  we  verified  first  the  linearity  of 
the  vacuum  photodiode  detector,  and  then  the  position  sensitivity 
of  its  high  and  low  frequency  response.  By  inserting  a  series  of 
broadband  filters  up  to  a  neutral  density  of  .8,  we  simulated  the 


reduction  in  power  produced  by  closing  the  aperture  diameter 
without  changing  the  ratio  of  laser  to  stimulated  power.  The 
measured  ratio  remained  within  a  fraction  of  one  standard 
deviation  from  unity.  The  detector  was  also  moved  from  side  to 
side  with  a  focussed  input  to  check  that  spatial  variations  in  its 
sensitivity  were  not  responsible  for  the  signal.  The  same  negative 
result  was  obtained:  the  detection  system  is  linear  and  sensitive 
only  to  the  power  transmitted  through  the  analyzing  iris. 

Since  our  goal  is  to  measure  the  dependence  of  the  stimulated 
power  on  the  aperture,  it  is  important  to  be  sure  that  no  portion 
of  the  laser  output  mode  is  scraped  off  at  any  point  in  the 
optical  transport  system.  (In  confirmation  of  this  fact,  we  note 
that  the  aperture  dependance  of  the  gain  disappeared  when  the 
monochromator  slits  were  closed  down,  thereby  scraping  off  a 
portion  of  the  focal  area).  During  set-up,  the  laser  mode  was 
centered  on  the  transport  mirrors,  and  the  monochromator  slits 
were  opened  to  ensure  that  when  the  analyzing  operture  was  open, 
all  of  the  stimulated  power  entered  the  detector.  When  the 
analyzing  aperture  was  placed  in  its  farthest  position  from  the 
undulator,  no  problem  was  observed  in  power  transport.  When  the 
aperture  was  placed  in  its  near  position  to  the  undulator,  and  set 
at  its  minimum  opening,  the  growth  of  the  output  beam  due  to 
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diffraction  did  result  in  perhaps  a  5%  power  loss  on  one  of  the 
beam  transport  mirrors. 

The  interpretation  of  the  results  demands  an  accurate 
knowledge  of  the  electron  and  optical  beam  sizes.  In  an  attempt  to 
simplify  the  comparison  with  theory,  considerable  effort  was 
expended  to  make  the  external  Argon  laser  run  in  a  single 
transverse  mode.  For  series  I,  this  effort  was  successful.  For  the 
later  series,  the  aged  laser  plasma  tube  performed  less  well,  and 
even  with  an  iris  internal  to  the  Argon  laser  oscillator  closed 
and  optimized,  the  beam  showed  significant  higher  order  mode 
power . 

The  laser  beam  sizes  were  repeatedly  measured  before  and 
during  series  II.  The  input  laser  mode  was  measured  with  a  Reticon 
diode  array  with  a  25  micron  diode  separation.  The  beam  shape  and 
the  full  width  at  half  maximum  was  recorded  at  five  positions 
about  the  focus,  both  in  the  vertical  and  the  horizontal 
dimensions.  These  measurements  were  fit  to  a  Gaussian  TEMoo  mode 
to  extract  the  beam  waist  and  focal  position  information.  The 
fits,  shown  in  figure  1,  show  moderate  agreement  with  the  data, 
but  cannot  be  used  reliably  to  estimate  the  higher  mode 
intensities.  The  astigmatism  of  the  Argon  laser's  asymmetric  cavity 


is  apparent;  as  expected,  the  vertical  focus  is  tighter  than  the 
horizontal,  and  occurs  slightly  before  it.  Figure  2  shows  the 
transverse  profile  of  the  intensity  at  the  center  of  the  undulator 
(nearly  Gaussian),  and  4.3  m  before  the  focus,  where  the  higher 
mode  power  becomes  apparent. 

For  the  data  analysis,  we  also  require  a  knowledge  of  the 
size  of  the  laser  beam  at  the  aperture.  This  information  is  obtained 
from  the  dependence  of  the  total  detected  laser  power  on  the 
diameter.  If  the  diameter  is  known  and  the  aperture  is  well 
aligned,  the  laser  beam  size  can  be  calculated  from  the  reduction 
of  the  transmitted  laser  power,  under  the  assumption  of  a  Gaussian 
beam. 

The  electron  beam  size  is  less  well  known.  It  is  measured  as 
shown  in  figure  3.  The  synchrotron  radiation  emitted  from  the  beam 
in  the  dipole  magnet  158  cm  in  front  of  the  center  of  the  optical 
klystron  is  imaged  onto  the  Reticon.  The  vertical  and  the 
horizontal  beam  shapes  are  typically  observed  to  be  Gaussian,  with  a 
width  which  varies  according  to  the  current  in  the  ring.  These 
dimensions  also  change  with  the  tunes  of  the  ring,  so  separate 
electron  beam  measurements  are  required  on  each  injection  during 
the  experiment.  No  diagnostics  are  available  on  the  beam  size  in 


the  center  of  the  optical  klystron.  This  is  an  unfortunate 
limitation  since  the  beta  functions  of  ACO  have  not  been  measurea 
with  the  sextupoles  in  operation.  In  the  absence  of  the  trajectory 
errors  induced  by  the  sextupoles,  the  theoretical  model  of  the  ACO 
optics  predicts  that  the  beam  sizes  change  by  less  than  25%  and 
6%,  respectively  in  the  horizontal  and  the  vertical  planes,  when 
moving  from  the  electron  beam  measurement  point  to  the  center  of 
the  OK.  These  numbers  can  be  taken  as  an  estimate  of  the 
systematic  error  in  the  electron  beam  size  determination. 

For  the  second  series  of  experiments  the  extraction  mirror 
had  been  removed  for  the  repair  of  an  intracavity  iris  in  the  same 
assembly.  We  were  therefore  obliged  to  run  with  the  probe  laser 
beam  exiting  the  vacuum  system  through  the  rear  oscillator  mirror. 
This  mirror  was  tilted  slightly  so  that  the  reflected  beam  could 
be  intercepted  on  a  beam  stopper.  The  transmission  of  this  mirror, 
a  high  reflector  at  6328  A,  is  about  80%  at  the  Argon  5145  A 
wavelength.  It  acts  as  a  defocussing  lens  with  a  focal  length  of 
about  -5.2  m.  These  experiments  also  had  to  be  performed  at  low 
energy  (166  MeV)  in  order  to  avoid  damaging  [  9 ]  the  oscillator 
mirrors  with  the  UV  generated  in  the  optical  klystron.  All 
measurements  were  done  with  the  5145  A  laser  line;  in  Series  I,  Nd 
=  80,  K  =  1.89,  and  Wo  =  670  microns;  in  Series  II,  Nd  =  70,  K  = 
.9,  and  Wo  =  610  microns. 


RESULTS  OF  THE  EXPERIMENTS 


One  data  scan  taken  in  series  I  is  compared 
with  theory  in  figure  4.  This  scan  was  taken  under  optimum 
conditions  with  well  characterized  electron  and  optical  beam 
sizes.  Here,  the  input  laser  beam  is  a  pure  TEMoo  mode,  and  the 
amplified  beam  propagates  freely  for  14.3  m  before  being  analyzed 
with  the  iris.  The  error  bars  can  be  estimated  from  the  accuracy 
with  which  the  two  consecutive  scans  overlap.  The  solid  curve  is 
the  theoretical  calculation  [l]  based  on  the  measured  distances 
and  beam  sizes;  no  fitting  has  been  done.  The  agreement  is 
excellent. 

This  experiment  is  the  first  to  demonstrate  the  existence  of 
the  off-diagonal  terms  in  the  gain  matrix  (transverse  mode 
mixing),  and  the  agreement  between  theory  and  experiment 
in  figure  4  establishes  the  validity  of  our  theoretical 

understanding  of  mode  mixing  in  FELs  in  the  weak  field  limit.  In 
the  calculation  of  the  theoretical  curve,  the  amplitudes  of  the 
first  few  terms  in  the  matrix  are  found  to  bejg^/g^j  =  .46,  and 
Ig  /g  /=  .17;  the  higher  terms  are  less  than  10%.  This  means  that 
the  amplitudes  added  to  the  two  higher  order  modes  are  .46  and  .17 
of  the  amplitude  added  to  the  lowest  order  mode,  respectively. 


(The  power  added  to  the  higher  order  modes  is  very  low  in  this 

-10  -* 
experiment,  10  and  10  of  the  laser  power  vs.  4X10  for 

the  lowest  mode,  because  the  gain  is  low  and  no 

higher  order  mode  power  is  present 

in  the  input  beam.)  It  is  clear  that  significant  mode  mixing  is 
ocurring  with  an  electron  beam  to  laser  beam  size  ratio  of  ^  =  25*/l 
=  0.96.  If  a  laser  of  moderately  high  gain  is  operated  with  such  a 
beam  ratio,  the  oscillator  mode  will  differ  significantly  from  a 
Gaussian. 

The  data  of  Figures  5  and  6  were  taken  during  the  initial 
investigations  of  the  effect.  For  these  curves,  the  laser  beam 
dimensions  are  not  as  well  characterized  as  for  Figure  4,  so 
theoretical  curves  are  not  available  for  comparison.  However, 
this  data  illustrate  certain  qualitative  features  of  the  mode 
mixing  effect  such  as  the  effect  of  changing  the  electron  beam 
size  and  the  resonance  parameter. 

Figure  5  show  the  results  of  two  series  I  measurements  taken 
as  the  stored  current  decayed  by  a  factor  of  1.9;  the  resonance 
parameter  was  changed  only  a  small  amount  (5v=  .1)  between  the  two 
scans.  Significant  differences  are  apparent  between  these  curves. 
This  data  demonstrates  the  effect  of  a  reduction  of  the  electron 


beam  size. 


In  ACO  the  transverse  dimensions  are  strongly  dependent  on 
the  current  stored  in  the  ring  because  of  the  multiple  Touschek 
effect  and  ion  trapping.  Between  these  two  scans,  the  electron 
beam  dimensions  have  fallen  from  approximately  <T  =  620  microns  at  79 
ma  of  average  beam  current  to  O'  =  480  microns  at  42  ma,  which  makes 
the  beam  ratios  *|=1.9  and  1.4,  respectively.  The  direct  effect  of 
the  current  on  the  size  of  the  gain  is  removed  by  the 
normalization,  leaving  only  the  effect  of  the  beam  size.  The  curve 
taken  with  the  smallest  electron  beam  shows  the  largest  gain 
depression  at  small  aperture.  This  confirms  the  theoretical 
expectation  that  smaller  electron  beam  sizes  result  in  more 
excitation  of  the  higher  order  transverse  modes.  All  off-diagonal 
terms  in  the  gain  matrix  tend  towards  zero  as  the  electron  beam 
becomes  large  --  see  equation  (2).  The  shape  of 

each  curve  depends  on  the  detailed  structure  of  the  mode  integrals 
(6)  and  their  relative  amplitudes.  The  shape  changes  as  the  beam 
size  drops  and  modes  of  progressively  higher  order  are  excited. 

This  effect  is  also  visible  in  figure  8. 

Figure  6  shows  a  later  set  of  data  taken  after  the  lifetime 
of  the  stored  beam  had  stabilized  to  a  larger  value.  These  three 
sets  of  points  were  taken  consecutively  at  approximately  the  same 
electron  beam  size  (®*=  420,  380,  and  350  microns),  but  widely 


varying  resonance  parameters  (i?=  3.6,  3.0  ,and  -3.0),  which 
correspond  to  the  f ourth-f rom-the-center  gain  and  absorption  peaks 
on  either  side  of  the  optical  klystron  spectrum.  No  systematic 
deviation  appears  between  these  curves,  and  we  conclude  that  the 
transverse  mode  coupling  is  independent  of  the  resonance 
parameter,  as  predicted  by  the  low  divergence  theory.  It  is 
apparent  from  (2)  that  in  the  low  divergence  limit,  the  spectral 
term  has  separated  out  as  an  independent  factor.  This  implies  that 
both  the  relative  amplitudes  and  phases  of  the  gain  matrix  and 
therefore  the  shape  of  the  gain-vs.-iris  curve  are  independant  of 
the  resonance  parameter. 

This  result  is  supported  by  a  second  type  of  measurement.  For 
a  given  aperture  diameter,  a  scan  was  taken  vs.  the  resonance 
parameter.  The  peak  positions  were  identified,  and  a  new  scan  was 
taken  for  another  diameter.  No  spectral  shifts  or  changes  in  the 
shape  of  the  spectrum  were  observed  over  a  range  of  iris  openings 
between  full  open  and  .9  on  the  horizontal  scale  of  the  figures. 

The  data  of  series  II  were  taken  with  the  primary  goal  of 
verifying  the  dependence  on  the  distance  of  the  aperture  from  the 
interaction  region.  The  evidence  of  multiple  mode  input  discussed 

above  implies  that  the  full  equation  (8)  must  be  used  to 


calculate  the  aperture  dependence  of  the  gain.  Unlike  the  case  of 
figure  4  where  an  <ib  initio  calculation  is  possible,  here  several 
fitting  parameters  would  be  needed  since  the  C  vector  of  the  mode 
amplitudes  and  phases  is  unknown.  Since  these  quantities  cannot 
be  obtained  from  the  measurements  of  figures  1  and  2,  and  since 
we  judge  the  results  of  a  multi-parameter  fit  to  be  useless,  no 
theoretical  values  are  calculated  for  the  Series  II  measurements. 
The  experimental  results  alone  are  sufficient  for  the 
identification  of  a  very  important  phenomenon  which  can  occur  in 
free  electron  lasers:  active  guiding. 

Figure  7  shows  a  composite  of  three  successive  measurement 
scans  taken  in  the  near  field  at  a  distance  from  the  center  of 
the  optical  klystron  of  2.88  m,  only  13  cm  from  the  output 
mirror.  The  large  fluctuations  in  the  amplitude  are  caused  by  the 
need  to  realign  the  aperture  each  time  a  new  opening  is  selected. 
A  large  increase  in  the  gain  is  observed  at  small  aperture 
openings.  A  set  of  data  was  also  taken  in  the  far  field  at  a 
distance  of  14.3  m,  and  is  presented  in  figure  8.  The  error  bars 
for  this  set  are  comparable  to  those  of  the  other  14.3  m  data  of 
figures  4-6. 

The  data  of  figure  8  were  taken  under  the  same  conditions 
as  those  of  figure  7.  There  is  a  large  and  significant  difference 
between  the  data  sets  taken  at  the  different  distances.  In  the 
near  field,  we  observe  that  the  amplitude  of  the  stimulated 
radiation  is  concentrated  at  the  center  of  the  amplified  laser 


mode.  The  on  axis  gain  is  enhanced  by  almost  50%.  In  the  far 
field  on  the  other  hand,  the  higher  order  modes  diverge  and  shift 
phase  with  respect  to  one  another  so  that  the  enhancement 
disappears.  A  change  in  the  distance  at  which  the  aperture  is 
placed  shifts  the  phases  of  the  T_ matrix  elements  of  equation 
(6),  and  therefore  changes  the  shape  of  the  gain-vs . -iris  curve. 
For  single  mode  input,  the  stimulated  amplitude  is  actually 
reduced  on  axis  in  the  far  field  as  we  have  seen  in  figure  4. 

The  significance  of  these  observations  lie  mainly  in  the 
implications  they  have  for  the  high  gain  mode  of  operation.  If 
the  amplification  is  large,  the  propagation  of  the  growing  wave 
becomes  dominated  by  the  near  field  diffractive  behavior  of  the 
stimulated  radiation.  We  infer  from  the  results  of  figure  7  that 
the  mode  which  propagates  in  a  high  gain  system  will  be  more 
strongly  concentrated  about  the  electron  beam  than  the  freely 
diffracting  unamplified  input  wave.  If  the  gain  per  Rayleigh 
range  is  large,  this  phenomenon  can  lead  to  guided  propagation  in 
which  the  shape  of  the  growing  mode  is  unchanged  throughout  the 
interaction  region.  This  phenomenon,  in  which  the  effects  of 
diffraction  are  compensated  by  strong  axial  amplification  [11], 
is  called  "gain  focussing". 

Although  the  imaginary  part  of  the  gain  matrix  (2)  cannot 
be  seen  in  a  low  gain  experiment,  in  a  high  gain  system  it  can 
shift  the  phase  of  the  propagating  wave  significantly.  As  Tang 
and  Sprangle  [12]  have  observed,  the  optical  phase  is 


retarded  by  the  interaction  on  the  positive  gain  side  of  the 
resonance  curve.  This  implies  that  the  stimulated  wave,  which  we 
have  observed  to  be  localized  to  the  axis  of  the  system,  tends  to 
distort  the  phase  fronts  of  the  input  wave  in  such  a  way  as  to 
produce  a  real  focussing  of  the  optical  mode.  Both  the  real  and 
the  imaginary  parts  of  the  gain  conspire  to  actively  guide  the 
node  in  the  high  gain  case  as  is  discussed  in  a  recent 
theoretical  paper  by  Scharlemann  et  al .  [13]. 

The  results  of  figure  4  demonstrate  that  our  present 
understanding  of  the  three  dimensional  effects  of  amplification 
is  numerically  correct.  The  results  of  figure  7  show  that  the 
gain  is  enhanced  on  axis  in  the  near  field.  Taken  together,  these 
two  results  lend  strong  experimental  support  to  the  theoretical 
prediction  of  active  guiding  in  FEL  systems. 
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T.  CONCLUSION 


This  experiment  is  the  first  capable  of  resolving  an 
excitation  of  the  higher  order  transverse  modes  in  a  free  electron 
laser,  and  has  established  the  existence  of  off-diagonal  terms  in 
the  transverse  gain  matrix.  The  technique  we  have  developed  is  a 
useful  diagnostic  of  the  transverse  mode  mixing  effect  which  is 
invisible  in  today’s  FEL  oscillators  due  to  the  low  gain,  but 
which  will  play  a  dominant  role  in  the  operation  of  the  high  gain 
devices  now  being  constructed  or  planned. 

High  gain  free  electron  lasers  will  show  strong  effects  of 
mode  mixing  in  any  system  in  which  the  electron  and  laser  beam 
sizes  are  comparable.  We  have  demonstrated  that  mode  mixing 
effects  are  typically  large  in  free  electron  lasers  but  that  they 
have  been  masked  in  the  existing  devices  by  the  low  value  of  the 
gain.  The  results  of  this  experiment  indicate  that  the  basic 
theory  is  valid,  and  that  it  is  complete  in  the  linear  regime. 

As  discussed  in  the  r  re  section,  these  results  lend  strong 
experimental  support  to  the  theoretical  prediction  of  active 
guiding  in  high  gain  free  electron  laser  systems. 
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FIGURE  CAPTIONS 


Figure  1 

The  laser  beam  was  measured  with  a  Reticon  diode  array  to 
determine  its  mode  characteristics.  The  beam  parameter  W(z)  measured  in 
Series  II  is  plotted  and  compared  to  a  TEMoo  Gaussian  fit.  The  beam  is 
well  focussed  into  the  klystron,  whose  position  is  indicated  at 
the  center  of  the  horizontal  axis. 

Figure  2 

Photographs  of  the  transverse  laser  profiles  during  series 
II,  taken  by  the  Reticon  at  two  longitudinal  positions:  at  the 
center  of  the  optical  klystron,  and  4.3  m  in  front  of  it.  In  the 
center  of  the  laser  interaction  region,  the  laser  profiles  can  be 
fit  well  with  a  Gaussian  shape.  However,  far  from  the  focus,  the 
profiles  acquire  a  significant  non-Gaussian  pedestal  which  is 
apparent  in  the  photo  on  the  left  (the  discontinuity  in  the 
vertical  profile  is  an  artifact  caused  by  a  bad  diode).  This  data 


shows  the  laser  beam  in  its  most  Gaussian  configuration  for  the 
Series  II  measurements.  Much  cleaner  Gaussian  behavior  was 
obtained  in  Series  I. 

Figure  3 

The  electron  beam  measurement  scheme.  Synchrotron  light 
emitted  in  the  bending  magnet  before  the  klystron  is  imaged  onto 
the  Reticon. 

Figure  4 

Measured  ratio  of  stimulated  power  to  incident  laser  power 
transmitted  through  an  analyzing  aperture  as  a  function  of  its 
diameter.  The  diameter  of  the  aperture  is  normalized  to  the  laser 
beam  wai  parameter  at  the  measurement  point,  and  the  gain  ratio 
is  normalized  to  unity  at  large  aperture  opening.  The  points  were 
measured  on  two  success. ve  scans;  the  error  bars  can  be  estimated  by 
the  deviation  of  the  points  from  each  other.  The  solid  curve  is 
the  theoretical  [l]  result,  using  the  measured  beam  parameters:  Wo 
=  670  microns  for  the  laser,  and  O'  =  320  microns  for  the  electron 
beam.  The  iris  was  placed  in  the  far  field  at  14.3  m  from  the 
optical  klystron,  where  the  optical  wave  had  grown  to  wq=  3*2  mm. 


Figure  5 


The  gain  is  measured  as  a  function  of  the  iris  diameter  as  in 
figure  4,  but  now  for  two  different  electron  beam  sizes.  The 
triangles,  taken  for  the  small  electron  beam  size,  show  the 
largest  deviation  from  unity,  and  therefore  the  largest  excitation 
of  the  higher  order  modes.  As  explained  in  the  text,  smaller 
electron  beam  sizes  result  in  larger  values  for  the  off-diagonal 
terms  in  the  transverse  gain  matrix.  The  shift  in  the  resonance 
parameter  between  the  two  curves  is  inconsequential,  as  shown  in 
figure  6. 

Figure  6 

The  gain  is  measured  as  a  function  of  the  iris  diameter  for  three 
different  values  of  the  resonance  parameter  which  correspond  to 
gain  and  absorption  peaks  displaced  by  four  optical  klystron 
fringes  (more  than  half  of  the  gain  curve)  from  the  center.  The 
data  all  fall  along  the  same  curve  within  the  error  bars,  and  we 
conclude  that  the  mode  mixing  effect  does  not  depend  on  detuning 
in  this  low  divergence  experiment.  The  changes  in  the  beam  size 
are  small  (and  unavoidable  due  to  the  current  decay  in  ACO  during 
these  measurements). 


Figure  7 

In  the  second  series  of  experiments,  the  gain  is  measured  as 
a  function  of  aperture  for*a  position  as  close  to  the  output  of  the 


optical  klystron  as  possible.  The  error  bars  are  large  here  because 
of  the  constant  realignment  of  the  aperture  position.  The  beam  sizes 


are  0"^  =  4,50^.  and  <T^  =  310^<t  for  the  triangles,  and  <TH  =  340 ^  and 
(Ty  =  290 jn  for  the  remaining  points.  The  input  laser  is  operating 
multimode  (see  figure  2). 


Figure  8 

Two  scans  of  the  gain  vs.  diameter  taken  in  the  far  field  of 
the  interaction  region  under  the  same  conditions  as  those  of  figure 
7.  The  deviation  of  these  results  from  those  of  figure  7  is  due  to 
the  difference  between  the  near-field  and  the  far-field  phases  of  the 
higher  order  excited  modes.  As  was  seen  in  figure  5»  "the  scan  with 
the  smaller  electron  beam  0^  =  360^/<  and  (7^  =  330^  (the  solid  points) 
shows  a  larger  deviation  from  unity  than  that  with  the  larger  beam 
(Tjj  =  450  u  and  0^  =  420  u  (the  open  squares). 
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Resume.  —  On  a  fait  fonctionner  un  laser  a  electrons  fibres  sur  anneau  de  stockaee  a  6  500  A  de  longueur  d’onde. 
L'effet  laser  a  ete  obtenu  avec  des  electrons  de  166  MeV  circulant  dans  l’anneau  de  stockage  ACO.  La  lumiere  emise 
par  les  electrons  lorsqu’ils  passent  a  travers  un  klystron  optique (une  version  modifiee  de  1'onduleur)  est  stockee  dans 
une  cavite  optique  a  faibles  pertes.  On  analyse  tout  d'abord  les  phenomenes  sous  le  seuil.  La  seconde  partie  decrit 
les  caracteristiques  spectrales  et  macrotemporelles  du  laser.  On  montre  que  la  puissance  moyenne  est  en  accord 
avec  la  limite  imposee  par  le  chauffage  du  paquet  d'electrons  (limite  de  Renien).  Finalement  une  augmentation  d'un 
factcur  100  de  la  puissance  Crete,  a  ete  obtenue  en  asservissant  le  laser  a  un  generateur  exterieur  de  signaux  basse 
frequence. 

Abstract  —  A  storage  ring  free-electron  laser  oscillator  has  been  operated  above  threshold  at  a  visible  wavelength 
~  6  500  A.  This  laser  was  obtained  with  electrons  of  166  MeV  circulating  in  the  storage  ring  ACO.  The  light 
emitted  by  the  electrons  passing  through  an  optical  klystron  (a  modified  version  of  an  undulator)  is  stored  in  a 
high-Q  optical  cavity.  First  some  below  laser  threshold  phenomena  are  reported.  Then  the  spectral  and  macro¬ 
temporal  structure  of  the  laser  are  described.  The  average  power  is  shown  to  be  consistent  with  the  limit  imposed 
by  the  heating  of  the  electron  beam  (Renieri’s  limit).  Finally  a  factor  of  100  enhancement  in  the  peak  power,  over  the 
Renieri's  limit,  has  been  obtained  by  driving  the  laser  with  an  external  low  frequency  generator. 


In  recent  years,  much  attention  has  been  given  to 
radiation  emitted  by  relativistic  electrons.  This  is 
principally  synchrotron  radiation  in  a  dipole  field  [1], 
in  an  undulator  [I]  or  from  a  free  electron  laser.  The 
free  electron  laser  (F.E.L.)  is  a  very  promising  source 
of  coherent  radiation  in  the  spectral  range  from  the  far 
I.R.  to  the  V.U.V. ;  the  earliest  proposals  for  such  a 
device  were  made  by  Ginzburg  [2]  and  Motz  [3]. 

There  are  basically  two  classes  of  F.E.L.  One  type 
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operates  in  the  «  collective  »  regime  where  the  elec¬ 
tronic  density  and  Coulomb  interaction  are  very 
strong,  the  electron  energy  is  rather  low  (a  few  MeV) 
and  the  radiation  is  emitted  in  the  far  I.R.  [4].  The 
order  class  operates  in  the  low  density  or  «  Compton  » 
regime,  where  Coulomb  interactions  are  negligible,  the 
energy  is  high  (I0-103  MeV)  and  the  emission  comes 
at  short  wavelengths  (I.R.,  Visible,  U.V.).  Our  experi¬ 
ment  belongs  to  the  second  category.  This  type  of 
F.E.L.  has  been  studied  by  Madey  [5]  and  his  group 
at  Stanford  where  the  first  laser  action  was  achieved 
in  1977  in  I.R.  [5].  Since  that  time  many  theoretical 
and  experimental  studies  have  been  done  on  this 
subject  [6]. 

Most  F.E.L.  experiments  involve  either  linear 
accelerators  and  microtrons  [7]  or  electron  storage 
rings  [8],  One  expects  the  storage  ring  D.E.L.  to  reach 
much  shorter  wavelengths,  particularly  in  the  U.V. 
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and  V.U.V.  spectral  range.  Our  experiment  is  the  first 
example  of  an  F.E.L.  working  on  a  storage  ring.  We 
have  recently  reported  the  first  operation  of  the 
laser  [9],  In  this  paper  we  described  in  more  detail  the 
operation  of  the  F.E.L.,  its  behaviour  as  a  driven 
pulsed  laser  and  the  interpretation  of  its  macro-tem¬ 
poral  structure.  We  have  used  the  Orsay  storage  ring 
ACO,  an  « old »  machine,  first  operated  in  1965. 
Its  characteristics  for  this  experiment  are  given  in 
table  I. 


Table  I. 

ACO  characteristics  in  the  F.E.L. 
Energy 

Circumference 
Bunch  to  bunch  distance 
Electron  beam  current  for  oscilla¬ 
tion 

R.M.S.  bunch  length  <r, 

R.M.S.  bunch  transverse  dimen¬ 
sions  <tx,  cry 

R.M.S.  angular  spread  a'x ,  a'y 
R.M.S.  relative  energy  spread 
Electron  beam  life  time 
R.F.  frequency 


experiment 

160-166  MeV 
22  m 
11  m 

16  to  100  mA 
0.5  to  1  ns 

0.3  to  0.5  mm 
0. 1  to  0.2  mrad 
0.9  to  1.3  x  10'3 
60  to  90  min 
27.2361  MHz 


Table  II. 

Optical  klystron  characteristics 
Overall  length  1.3  m 

Type  SmCo5  permanent  magnets 


Undulators 

Number  of  periods  2  x‘ 

Period  78 

Transverse  pole  width  100  mm 

Minimum  magnetic  gap  33  mm 
Maximum  field  0.31  T 

K  0  to  2.3 

K  in  the  F.E.L.  experi¬ 
ments  1.1  to*  1.2 

Dispersive  section 
Length  240  mm 

Maximum  field  0.58  T 

Nd  (see  text)  70  to  100 

Nd  in  the  F.E.L.  experi- 


The  core  of  an  F.E.L.  is  an  undulator  (periodic 
transverse  magnetic  field)  on  passage  through  which 
high  energy  electrons  radiate  a  series  of  harmonics 
of  a  fundamental  line  of  wavelength.  A,  given  by  : 

X  =  J7  1  +~  +  V202]  0) 

where  y  is  the  electron  energy  normalized  to  their 
energy  at  rest,  A0  the  undulator  period  and  0  the  angle 
of  observation  with  respect  to  the  electron  trajectory 
axis.  K  is  the  deflection  parameter  of  the  undulator 
given  by  : 

€  B  / 

K  =  - — 2—2.  =  93.4  s  ;  in  s.I.  units  (2) 

2  7 rm0  c 

where  B0  is  the  peak  undulator  magnetic  field,  e  and 
m0  the  electron  charge  and  mass.  The  characteristics 
of  our  undulator/optical  klystron  are  listed  in  table  II. 

A  light  beam  of  convenient  wavelength  in  coinci¬ 
dence  with  an  electron  beam  is  amplified  when  pass¬ 
ing  through  an  undulator.  The  magnitude  of  the 
optical  gain  was  first  calculated  by  Madey  [5].  We 
have  measured  the  gain  of  the  Orsay  undulator  made 
with  17  periods  of  7.8  cm  [10]  placed  on  the  ACO 
storage  ring  working  at  240  MeV.  The  stimulating 
light  was  an  argon  laser  at  5  145  A.  The  measured 
gain  value  was  near  2  x  10" 4  [II],  This  gain  is  too 
small  to  achieve  laser  oscillation;  therefore,  we  have 
transformed  our  undulator  into  an  optical  klys¬ 
tron  [12]  (O.K.).  This  is  done  by  replacing  the  3  central 


periods  of  the  undulator  by  a  longer  period,  analogous 
to  a  3-pole  wiggler,  called  the  dispersive  section. 
The  role  of  this  section  is  to  «  bunch  »  the  electronic 
density  faster  than  a  regular  undulator  [12).  There¬ 
fore,  the  optical  klystron  gain,  Go  k  is  larger  than  the 
gain  Gunll  of  the  undulator  of  same  length  by  a  maxi¬ 
mum  factor  of  : 


where  <te/£  is  the  R.M.S.  relative  energy  spread  of  the 
electrons  circulating  in  the  ring 

The  synchrotron  light  emitted  by  the  OK.  exhibits 
a  characteristic  modulation  [13]  which  is  shown  on 
figure  la.  The  optical  gain  (Fig  lb)  which  by  virtue 
of  the  Madey’s  theorem  [14]  is  the  derivative  of  the 
spontaneous  emission  with  respect  to  energy,  also 
appears  modulated.  For  a  given  fundamental  wave¬ 
length  GokJGuluj  depends  on  the  magnetic  gap  and  on 
the  energy  spread.  The  measured  value  confirms  a 
predicted  factor  of  4  gain  enhancement  at  240  MeV 
and  1  x  10" 3  energy  spread. 
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Fig  I.  — (a)  On  axis  spectral  distribution  of  the  optical 
klystron  spontaneous  emission  measured  using  240  MeV 
electrons  and  a  magnetic  gap  of  34.4  mm.  (b)  The  corres¬ 
ponding  gain  curve  versus  magnetic  gap  measured  with  an 
external  laser. 


In  this  paper  we  present  results  obtained  with  our 
optical  klystron  working  at  a  low  electron  energy 
(166  MeV).  This  energy  range  was  chosen  in  order  to 
minimize  the  cavity  mirror  degradation  due  to  V.U.V. 
emissions  [11]  from  the  undulator  and  dispersive 
section. 

The  light  emitted  by  the  O.K.  was  stored  in  the 
highest-Q  optical  cavity  that  we  could  obtain 
(Table  III).  The  measured  average  round  trip  losses 


Table  III. 


Optical  cavity  characteristics 
Length 

Mirror  radius  of  curvature 
Rayleigh  range 
Wavelength  of  maximum  Q 
Average  mirror  reflectivity  at  6  328  A 
Round  trip  cavity  losses  at  6  328  A 
Mirror  transmission 


5.5  m 
3  m 
1  m 

620  to  680  nm 
99.965  % 

7  x  10~4 
3  x  1(T5 


for  the  mirrors  that  we  used  were  near  2  x  10“ 4  at 
6  328  A  when  the  mirrors  were  in  the  air.  However 
when  placed  in  the  ultra-vacuum  of  the  cavity  and 
irradiated  with  the  light  produced  by  the  undulator 
or  optical  klystron  the  losses  grow  rapidly  [11],  At 
240  MeV  the  degradation  drives  the  losses  almost 
instantaneously  above  1  or  2  x  10“ 3.  However  it  is 
possible  to  operate  the  O.K.  at  lower  electron  energy 
by  increasing  the  magnet  gap  (thereby  decreasing  K) 
in  order  to  keep  the  emitted  wavelength  near  6  500  A, 
wavelength  of  maximum  reflectivity  for  the  mirrors 
used  in  these  experiments.  A  lower  K  decreases  the 
harmonics  content  of  the  undulator  emission  and 
therefore  the  U.V.  flux  received  by  the  mirrors.  The 
O.K.  was  finally  operated  at  K  =  1.1  to  1.2  with  a 
dispersive  section  characterized  by  [13]  Nd  =  95 
where  (Vd  is  the  number  of  laser  optical  wavelengths 
passing  over  an  electron  in  the  dispersive  section. 

At  a  low  energy  of  150  MeV,  after  an  initial  degra¬ 
dation  the  losses  were  stabilized  near  (7  ±  1).10“4. 
At  this  energy  the  optical  gain  is  only  slightly  smaller 
than  the  losses  and  very  interesting  «  below  threshold  » 
phenomena  may  be  observed  (§  1).  At  higher  energy 
(160-166  MeV)  an  intense  amplification  of  the  stored 
emission  is  observed  (laser  effect).  The  various  charac¬ 
teristics  of  the  laser  are  detailed  in  §  2  and  discussed 
in  §  3. 


1.  Below  threshold  phenomena  (150  MeV). 

When  the  optical  cavity  mirrors  are  set  parallel  to 
each  other,  the  spontaneous  emission  of  the  O.K.  is 
stored  in  the  cavity.  The  resulting  optical  spectrum 
(Fig.  2,  lower  curve)  is  the  spontaneous  emission  spec¬ 
trum  modified  by  the  phase  shifts  of  the  various  cavity 
transverse  modes  in  which  it  is  stored  [15]  and  by  the 
output  mirror  transmission.  When  the  cavity  length 
is  exactly  matched  with  the  electron  round-trip  time 


CAVITY  TUNED 


Fig.  2.  —  Fine  structure  spectrum  of  the  stored  sponta¬ 
neous  emission  from  the  optical  klystron.  The  « cavity 
tuned  »  curve  is  obtained  for  a  perfect  time  synchronism 
between  light  and  electron  pulses  in  the  cavity.  The  «  cavity 
detuned  »  curve  is  obtained  for  complete  desynchronism. 
Gain  and  absorption  appear  on  the  higher  and  lower  wave¬ 
length  side  of  the  maximum. 


in  the  ring  the  cavity  is  considered  to  be  « tuned  ». 
Then  amplification  of  the  stored  spontaneous  emis¬ 
sion  occurs  through  the  gain  mechanism  mentioned 
above.  One  clearly  sees  in  figure  2  an  amplification  at 
A  =  6  470  A  by  a  factor  close  to  2  and  an  absorption 
at  A  ~  6  440  A  by  a  factor  0.8.  This  amplification 
occurs  at  several  wavelengths  along  the  spontaneous 
emission  curve  with  a  maximum  depending  mostly 
on  the  mirrors  reflectivity  curve.  The  spectral  distribu¬ 
tion  of  the  amplified  beam  is  no  longer  proportional 
to  the  derivative  of  the  stored  spontaneous  emission 
because  of  the  non-linearity  of  the  multipass  amplifi¬ 
cation  when  the  gain  value  is  close  to  the  cavity  losses. 
Also  the  noise  that  can  be  seen  at  the  top  of  the  amplifi¬ 
ed  curve  is  clearly  related  to  the  fact  that  we  are  very 
near  the  laser  threshold.  A  small  gain  variation  pro¬ 
duces  a  large  amplification  variation.  By  aligning  the 
optical  cavity  axis  and  the  electron  beam  trajectory 
within  0.1  mm,  a  maximum  amplification  ratio  of  3 
has  been  observed  for  a  50  mA  ring  current.  Figure  3 
shows  horizontal  and  vertical  transverse  profiles 
of  the  stored  spontaneous  emission  recorded  at 
A  6  470  A,  corresponding  to  the  most  amplified 
wavelength.  The  net  amplification  (obtained  by  sub¬ 
tracting  the  cavity  detuned  curve  from  the  cavity 
tuned  one)  has  a  profile  very  close  to  the  fundamental 
cavity  TEM00  profile. 

The  amplification  is  a  sharp  function  of  the  adjust¬ 
ment  of  the  cavity  length.  The  dependence  of  its  value 
versus  the  detuning  of  this  length  is  shown  on  figure  4. 
The  amplification  as  a  function  of  the  mirror  displa¬ 
cement  from  the  tuned  position,  5,  can  be  computed 
by  summing  the  contributions  over  the  electron  posi¬ 
tions  (characterized  by  <50)  inside  the  electron  bunch 
and  over  a  quasi-infinite  number  of  pass,  n,  in  the 
optical  cavity.  Let  F(S)  be  that  sum,  then  : 


F(6)  =  (I  -  R2) 
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Fig.  3.  —  Horizontal  and  vertical  profiles  of  the  stored 
spontaneous  emission  at  the  wavelength  of  maximum  gain. 
The  difference  between  the  «  cavity  tuned  »  and  «  cavity 
detuned  »  curves  show  a  good  agreement  with  the  cavity 
TEM00  profile. 


L0NGITU01NAL  MIRROR  TRANSLATION 

Fig  4.  —  Maximum  on  axis  amplification  ratio  versus 
cavity  length.  The  dark  curve  shows  the  theoretical  predic¬ 
tion  for  0.23  m  R.M.S.  bunch  length,  3.6  x  10"4  average 
loss  per  mirror  and  5. 1  x  10" 4  peak  gain  per  pass. 


where  R  is  the  average  reflectivity  per  mirror,  f(50) 
is  the  normalized  shape  function  of  the  electron 
bunch  and  qp  the  peak  gain  at  the  wavelength  consi¬ 
dered.  From  (4)  one  deduces  : 


F{5  =  x>)  =  1  and  F(5  =  0) 


W o)  n 
1  _ 

1  -  R 2 
(5) 


valid  for  gp  <  1  -  R2  =  p. 


For  a  uniform  distribution  of  the  electrons  inside 
the  bunch  (5)  reduces  to  : 

F(5  =  0)  =  [1  -  gJP)-'. 

This  would  yield  gJP  =  0.60  in  the  case  of  figure  4 
(amplification  by  2.5).  In  our  case  (Gaussian  electronic 
distribution)  the  integral  has  been  computed  nume¬ 
rically  using  the  measured  R.^LS.  bunch  length  of 
0.23  m.  The  best  fit  to  the  experimental  results  shown 
in  figure  4  was  obtained  for  P  =  1.2  x  10~4  (close 
to  the  measured  value  of  7.0  x  10"4)andgp//>  =  0.71. 

Very  good  agreement  has  also  been  found  in  the  tail 
of  this  detuning  curve  showing  a  significant  amplifica¬ 
tion  of  1.02  for  a  detuning  as  high  as  0.8  mm. 

2.  Above  threshold  phenomena. 

At  166  MeV  laser  oscillation  is  obtained  by  a  careful 
alignment  of  the  electron  beam  on  the  cavity  axis 
(within  0.1  mm  and  0.1  mrad)  and  maximization  of 
emission  as  a  function  of  the  storage  ring  R.F.  cavity 
voltage  and  the  optical  cavity  length.  Laser  operation 
requires  very  precise  synchronism  between  the  light 
pulse  round  trip  frequency  and  the  bunch  revolution 
frequency.  To  avoid  backlash  and  mirror  misalign¬ 
ment.  fine  tuning  was  performed  by  slightly  changing 
the  R.F.  frequency  instead  of  translating  the  mirrors. 
Laser  oscillation  lasted  typically  1  h  after  each  elec¬ 
tron  injection. 

Due  to  the  weak  gain/loss  ratio  the  cavity  length 
(or  R.F.  frequency)  is  a  very  crucial  parameter. 
Figure  5  shows  two  «  detuning  curves  »  of  laser  power 
(normalized  to  the  maximum)  versus  R.F.  frequency 
variation  and  equivalent  mirror  displacement.  Curve 
(a)  has  a  3.4  pm  F.W.H.M.  of  equivalent  mirror 
displacement;  curve  (b).  recorded  much  closer  to 
laser  threshold,  has  only  a  1.6  pm  F.W.H.M.  The  shift 
in  displacement  between  curves  (a)  and  (b)  is  probably 


RF  frequency  (Ml) 


Fig.  5.  —  Normalized  laser  power  as  a  function  of  R.F. 
frequency  and  equivalent  mirror  displacement.  Curve  (a) 
is  recorded  close  to  the  maximum  gain/loss  ratio  and  curve(b) 
close  to  the  laser  threshold.  The  shift  between  the  two  curves 
is  probably  due  to  some  slight  cavity  length  drift.  The  upper 
scale  has  to  be  multiplied  by  a  factor  2 
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due  to  a  slow  cavity  length  drift  (a  temperature  drift 
of  0.02  °C  in  30  min  would  be  sufficient).  However  the 
length  of  the  cavity  was  stable  enough  to  allow  laser 
operation  without  adjustments  during  typically  more 
than  half  an  hour. 

The  horizontal  and  vertical  transverse  profiles  of  the 
laser  mode  are  very  close  to  the  cavity  TEM00  mode 
( Fig.  6)  as  expected  from  the  below  threshold  measure¬ 
ment.  The  small  observed  discrepancies  may  arise 
from  some  residual  instability  of  the  laser  too  close  to 
threshold  or  from  non  uniform  mirror  reflectivity. 


Fig.  6.  —  Experimental  laser  horizontal  ( - )  and  verti¬ 

cal  (-  -  -)  transverse  profiles  and  the  calculated  (...)  cavity 
TE\l0o  profile 

The  laser  exhibits  a  particular  time  structure  :  (i)  it 
is  pulsed  at  the  ACO  storage  ring  frequency  (27.2  MHz 
with  two  bunches  in  the  ring),  (ii)  its  structure  on  a 
long  time  scale  (of  the  order  of  1  s)  is  quite  irregular 
as  shown  by  the  photograph  in  figure  7a,  (iii)  the  laser 
may  be  triggered  by  an  external  low-frequency  gene¬ 
rator  (Fig.  7b).  This  will  be  discussed  in  section  3. 


(b)  (a) 


Fig.  7.  —  Laser  time  structure  over  a  200  ms  total  interval, 
(a)  «  natural  »  time  structure;  (b)  time  structure  when  the 
electron  beam  transverse  position  is  modulated  (upper 
trace).  The  vertical  scale  is  1/4  the  scale  on  (a)  (see  text). 

Figure  8  presents  two  spectra  :  (a)  is  recorded 
without  amplification  (optical  cavity  completely  detun¬ 
ed),  (b)  is  recorded  at  laser  operation  (cavity  tuned). 


Fig.  8.  —  Spectra  of  the  cavity  output  radiation  under  two 
conditions  :  (a)  cavity  detuned  (no  amplification),  (b)  cavity 
tuned  (laser  on). 


For  case  (b),  the  laser  oscillates  at  three  wavelengths, 
with  a  strongly  dominant  one  at  ).  =  6  476  A;  each 
wavelength  is  located  at  a  maximum  of  the  gain  versus 
wavelength  curve  [5],  Typical  laser  lines  are  Gaussian 
if  averaged  over  a  long  time  scale  of  ^  1  s  with  2  to 
4  A  F.W.H.M. 

The  central  wavelength  of  any  line  is  always  equal  to 
the  wavelength  of  maximum  emission  of  spontaneous 
emission  with  the  cavity  completely  detuned  (no 
amplification)  plus  0.15  of  the  wavelength  interfringe 
distance  (see  Fig.  8)  instead  of  0.25  as  predicted  from 
Madey’s  theorem  [14],  This  discrepancy  is  probably 
due  to  the  transverse  multimode  content  of  the  spon¬ 
taneous  emission  stored  in  the  cavity;  laser  operation 
is  only  achieved  on  the  TEM00  mode.  Laser  tunability 
was  obtained  between  640  nm  to  655  nm  by  changing 
the  magnetic  gap  (equivalent  to  a  change  of  K  in 
Eq.  (1)).  The  range  of  tunability  is  limited  for  the 
moment  by  the  mirror  reflectivity  bandwidth. 

The  laser  spectrum  has  also  been  recorded  on  a 
Reticon  diode  array.  The  laser  was  driven  by  the  Reti- 
con  «  start »  signal  in  order  to  synchronize  the  laser 
pulses  with  the  Reticon  read  out  For  a  driving  fre¬ 
quency  of  ~  800  Hz,  as  discussed  in  section  3,  the 
laser  pulse  intensities  were  unstable  and  it  was  observ¬ 
ed  that  the  higher  the  pulse  the  narrower  its  spectral 
line.  For  the  higher  pulses  the  observed  spectral 
width  is  less  than  1  A  (Fig.  9)  well  below  the  above 
mentioned  spectral  width. 

In  the  experiments  the  gain  is  not  directly  propor¬ 
tional  to  the  electron  current,  mainly  due  to  the  ano¬ 
malous  bunch  lengthening  and  energy  spreading  at 
high  current.  These  effects  make  that  the  gain  versus 
ring  current  reaches  a  broad  maximum  near  50  mA 
and  then  decreases.  However  despite  this  effect  and 
despite  the  irregular  pulsed  time  structure  described 
above  the  laser  average  power  is  very  stable  and 
depends  almost  linearly  on  the  ring  current  as  shown 
by  the  recording  of  figure  10.  This  fact  shows  clearly 
that  the  free  electron  laser  has  saturated  and  reached 
its  equilibrium  average  power. 


Fig.  9.  —  Recording  of  the  spectral  distribution  of  a  laser 
pulse  on  a  reticon  diode  array.  The  laser  was  driven  at  a 
frequency  of  8.8  kHz  (see  text).  Each  diode  corresponds  to 
0.3  A  in  wavelength  which  is  also  the  spectral  resolution  of 
the  monochromator  used. 
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Fig:  10.  —  Recording  of  the  laser  average  output  power 
together  with  the  ring  current.  The  ratio  of  the  two  quantities 
is  a  constant  on  most  of  the  laser  operation  range. 

3.  Laser  saturation  and  peak  power  enhancement 

It  has  been  shown  that  storage  ring  F.E.L.  would 
saturate  by  laser  induced  electron  energy  spread 
which  would  decrease  the  gain  by  inhomogeneous 
broadening  [16].  The  energy  spreading  is  also  expected 
to  lengthen  the  electron  bunch  and  decrease  the  peak 
density  and  therefore  the  peak  gain  [17].  Renieri  [16] 
has  shown  that  the  average  power  (PJ  is  a  fraction  q 
of  the  total  synchrotron  radiation  (P5)  emitted  by  a 
single  bunch  in  the  whole  storage  ring.  r\  represents 
the  laser  efficiency  since  it  is  the  ratio  between  the 
laser  power  and  the  Radio  Frequency  (R.F.)  power 
given  by  the  R.F.  cavity  to  the  electron  beam  to  com¬ 
pensate  for  the  power  lost  by  electrons  by  synchrotron 
radiation  emission,  rj  has  been  calculated  for  a  weakly 
saturated  optical  klystron  to  be  (18)  : 

11  *  ,v),)L"<»o/z’> 

where 

_  output  mirror  transmission 
'■  round  trip  cavity  losses  ' 

is  the  cavity  efficiency,  /  is  the  optical  klystron 


modulation  rate  [13],  .V  is  the  number  of  periods  of 
one  of  the  two  identical  undulators  constituting  the 
optical  klystron.  g0/P  is  the  ratio  between  the  start-up 
peak-gain  and  the  round  trip  cavity  losses. 

The  laser  output  power  always  decreases  with  the 
electron  beam  current  and.  except  in  regions  too  close 
to  threshold,  is  usually  proportional  to  the  current 
(see  Fig.  10).  Such  a  result  is  in;good  agreement  with 
the  theory  since  />s  is  proportional  to  the  current  but  rj 
is  only  weakly  current  dependent.  A  typical  75  pW 
average  output  power  has  been  recorded  at  50  mA 
current  of  166  MeV  electrons;  this  corresponds  to  an 
efficiency  rj  =4.8  x  10" 5  with  respect  to  the  1.55  W 
of  single  bunch  total  synchrotron  radiation.  Taking 
the  following  measured  quantities,  rjc  —  0.047  (due 
to  the  high  absorption  losses  in  the  dielectric  layers  of 
the  mirrors  compared  to  the  transmission),  /  =  0.7, 
N  +  Nd  ~  100  and  assuming  g0/P  =  3,  one  calcu¬ 
lates  rj  =  5.7  x  10" 5. 

Considering  the  errors  in  the  calibration  of  ring 
current  and  power  meter  and  the  uncertainty  in  the 
determination  of  g0/P ,  one  concludes  in  the  good 
agreement  of  the  measured  average  power  with  the 
one  predicted  by  theory  [18].  This  theory  does  not 
consider  the  anomalous  bunch  lengthening  and  energy 
spreading  as  observed  on  ACO  and  on  almost  every 
storage  ring.  The  most  striking  example  of  this  ano¬ 
malous  behaviour  was  the  observed  bunch  shorten¬ 
ing  and  electron  energy  spread  reduction  at  the  initial 
laser  start-up.  Although,  as  we  shall  see  later,  the 
bunch  lengthens  when  the  laser  saturates  as  predicted 
by  the  theory. 

The  rather  low  average  output  power  is  therefore 
the  result  of  the  weak  total  synchrotron  radiation 
power  at  166  MeV,  the  low  cavity  efficiency  and  the 
high  sensitivity  of  the  optical  klystron  to  the  energy 
spread 

Figure  7a  shows  the  observed  F.E.L.  time  structure 
over  a  200  ms  time  interval.  The  irregular  time  struc¬ 
ture  of  the  laser  is  the  result  of  the  large  difference  in 
magnitude  between  the  fast  laser  rise  time.  rm,  of  the 
order  of  50  to  100  ps,  and  the  slow  response  time  ts 
of  the  electron  beam,  which  is  200  ms  at  166  MeV.  It 
has  long  been  anticipated  that  at  the  instant  of  laser 
turn-on  in  a  storage  ring  laser,  the  power  would 
dramatically  overshoot  the  equilibrium  value  and 
perharps  oscillate  until  the  electron  beam  had  a 
chance  to  stabilize  at  a  level  which  saturates  the  gain. 
Our  results  show  that  such  an  equilibrium  may  be 
very  difficult  to  obtain  due  to  fluctuations  in  the  gain. 
Elleaume  [19]  has  recently  studied  the  stability  of  the 
laser,  and  found  that  although  the  equilibrium  situa¬ 
tion  is  stable,  it  is  only  weakly  damped  The  laser  ope- 
ration  is  resonant  at  a  period  Tt  of  the  order  of 
2  7r  v'r,  xjl  ( 1 5-25  ms  for  these  experiments  as  shown 
in  Fig.  7a).  and  even  a  very  small  fluctuation  of  the 
electron  energy  spread,  the  bunch  length,  or  the  elec¬ 
tron  beam  position  with  this  period  will  cause  the 
laser  to  pulse  on  and  off. 
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In  more  than  IS  successive  experimental  runs, 
involving  ring  currents  ranging  from  80  to  15  mA, 
continuous  operation  has  never  been  observed.  The 
laser  always  pulses  in  an  irregular,  pseudo-periodic 
fashion  with  a  period  of  10  to  25  ms.  Sometimes,  the 
peaks  are  more  or  less  modulated  at  the  ring  synchro¬ 
tron  frequency  (  ~  13  kHz). 

As  a  consequence  of  this  high  sensitivity  to  gain 
fluctuations,  the  laser  may  be  triggered  by  an  external 
low-frequency  signal.  Stable  pulsed  operation  has 
now  been  achieved  by  modulating  either  the  transverse 
beam  position  or  the  R.F.  frequency.  Figure  7b  shows 
a  single  trace  of  the  laser  output  when  the  transverse 
position  of  the  electron  beam  is  periodically  displaced. 
The  peak  power  has  been  increased  by  a  factor  of  4 
over  that  of  figure  7a;  the  horizontal  time  scale  is  the 
same.  This  result  was  obtained  by  turning  on  the  gain 
long  enough  for  the  pulse  to  develop,  heat  the  beam, 
and  decay,  and  turning  it  off  again  for  a  long  time  to 
allow  the  beam  to  cool  down.  A  vertical  displacement 
of  160  pm  is  enough  to  turn  off  the  laser.  The  beam 
motion  is  adiabatic  with  respect  to  the  betatron 
motion  in  the  ring,  so  the  transverse  dimensions  are 
not  affected  by  the  periodic  orbit  distortion.  Stable 
driven  pulses  have  also  been  obtained  by  modulating 
the  R.F.  frequency  back  and  forth  from  the  synchro¬ 
nism  condition,  but  this  procedure  presents  the  disad¬ 
vantage  of  inducing  synchrotron  oscillation  of  the 
electron  bunch. 

When  driven  by  beam  displacement  as  shown  in 
figure  7b.  pulse  widths  less  than  1  ms  are  obtained 
with  periods  from  1  to  60  ms.  No  change  in  the  average 
power  is  observed  over  this  range,  although  the  peak 
power  can  be  adjusted  by  more  than  a  factor  of  tea 
A  peak  power  100  times  larger  than  the  average  power 
has  been  obtained  so  far. 

This  pulsed  operation  of  the  laser  allows  another 
test  of  the  saturation  theory.  Figure  11  presents  the 
laser  pulse  structure  together  with  a  signal  propor- 
tionnal  to  the  bunch  length.  The  bunch  length  mea¬ 
surement  technics  is  described  in  [17].  A  negative 
variation  corresponds  to  a  bunch  lengthening  One 
clearly  sees  a  fast  bunch  lengthening  induced  by  the 
laser  pulse  followed  by  a  slow  relaxation  time.  The 
bunch  length  curve  appears  to  be  the  integral  of  the 
laser  pulse  shape.  These  results  were  precisely  pre¬ 
dicted  to  occur  with  energy  spread  [18]  and  therefore 


Fig.  11. — Time  recording  of  the  laser  intensity  (upper 
trace)  together  with  the  bunch  length  evolution  (lower 
trace).  A  negative  signal  corresponds  to  a  lengthening  of 
the  bunch.  The  total  time  interval  is  100  ms  for  the  left  picture 
and  5  ms  for  the  right  one. 

bunch  length  since  energy  spreading  automatically 
induces  a  bunch  lengthening  in  a  time  of  the  order 
of  the  synchrotron  period  (0.08  ms)  which  is  much 
shorter  than  the  pulse  duration. 

Therefore  the  laser  time  structure  can  be  understood 
by  a  simple  model  [  1 9]  where  the  fast  increase  in  energy 
spreading  is  responsible  of  the  off-equilibrium  beha¬ 
viour  of  the  system.  The  energy  spread  evolution 
has  been  verified  experimentally  through  the  time- 
resolved  recording  of  the  electron  bunch  length. 
Moreover  we  have  shown  that  this  off-equilibrium 
laser  behaviour  offers  an  opportunity  to  drive  the 
laser  intensity  by  an  external  trigger  and  to  obtain 
a  pulsed  laser  without  any  lost  of  average  power. 

The  experiments  will  continue  in  the  future  to 
analyse  the  short  time  structure  of  the  laser  with  a 
streak  camera.  Also  higher  gain/loss  ratios  are  expect¬ 
ed  to  be  obtainable  at  higher  energy  and  from  the  use 
of  better  mirrors,  a  second  R.F.  cavity  and  positrons 
instead  of  electrons.  This  should  extend  the  range  of 
tunability  of  the  laser  towards  the  U.V.  regioa 
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Abstract 

A  storage  ring  f ree-electron  laser  oscillator  has  been  operated  at  a  visible  wavelength 
X  =  6500  A 

Introduction 

The  purpose  of  the  Orsay  experiment  is  to  test  the  feasability  of  a  visible  Storage 
Ring  Free-Electron  Laser  (F.E.L.).  The  realization  of  the  permanent  magnet  undulator  NOEL 
was  achieved  in  January  1982  (D.  Detailed  results  concerning  the  spontaneous  emission, 
the  gain  measured  using  an  external  laser,  the  laser  induced  bunch  lengthening  have  been 
reported  at  the  1982  Bendor  meeting  (2)  .  Improvement  of  the  original  low  gain  (a  few 
parts  in  10-4  per  pass)  was  achieved  in  July  1982  by  modifying  the  undulator  into  an 
optical  klystron  (3,4)  (O.K.).  The  O.K.  optimization  and  early  results  on  the  spontaneous 

emission  have  also  been  reported  at  the  Bendor  meeting(2).  The  gain  has  been  measured  using 
an  external  argon  laser  (5)  and  laser  induced  bunch  lengthening  has  been  reported  (*>). 

All  the  measurements  turned  out  to  be  in  a  very  good  agreement  with  the  classical  theory (4>  . 
The  maximum  gain  (at  the  best  electron  beam  current  and  energy)  is  now  about  10-3  per  pass. 
Therefore,  operation  of  the  F.E.L.  above  threshold  requires  very  high  reflectivity  mirrors 
having  losses  less  than  5  10-4  (or  reflectivity  >  99.95  %)  .  Mirrors  having  losses  less  than 
10-4  have  been  tested,  unfortunately  they  have  shown  reflectivity  degradation  mainly  due  to 
the  U.V.  part  of  the  undulator  or  optical  klystron  synchrotron  radiation  (2<7)>  To  mj.nimize 
this  degradation,  experiments  were  carried  out  at  150,  160  and  166  MeV  electron  energy  (as 
opposed  to  the  original  240  MeV  optimization  energy  for  the  undulator)  in  Spring  1983.  Des¬ 
pite  an  initial  quick  degradation  from  10-4  up  to  3.5  10-4  average  losses  per  mirror,  sta¬ 
bility  of  the  cavity  Q  was  observed  at  this  last  value.  The  gain  at  150  MeV  turned  out  to  be 
less  than  6  10-4  preventing  from  crossing  the  oscillation  threshold.  However,  operation  of 
the  storage  ring  at  160  MeV  and  then  at  166  MeV  brought  an  improvement  of  the  electron 
beam  quality  (mainly  bunch  length  and  energy  spread)  and  allowed  the  oscillator  to  cross 
threshold . 


Description  of  the  experiment 

We  have  used  the  Orsay  storage  ring  ACO,  an  older  machine  first  operated  in  1965.  Its 
characteristics  for  this  experiment  are  given  in  the  following  table. 


Energy 

Circumference 

Bunch  to  bunch  distance 

Electron  beam  current  for  oscillation 

RMS  bunch  length  oi 

RMS  bunch  transverse  dimensions  o  ,a 
RMS  angular  spread  o'  ,  o'  x  ^ 

RMS  relative  energy  spread  y 
Electron  beam  life  time 


160-166  MeV 
22  m 
11  m 

16  to  100  mA 
0.5  to  1  ns 
0.3  to  0. 5  mm 
0.1  to  0.  2  mrad 
0.9  to  1.3  10"3 
60  to  90  minutes 


The  optical  klystron  main  characteristics  are  : 
Full  length  1.33  m 
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t'r.dulators 

Number  of  periods 
Period 

Transverse  pole  width 
Minimum  magnetic  gap 
Maximum  field 
K 

K  in  the  F.E.L.  oscillator  experiment 

Dispersive  section 
Length 

Maximum  field 
Nd(4) 

Nd  in  the  F.E.L.  oscillator  experiment 


2x7 
78  mm 
1 00  mm 
33  mm 
C  .  3 1  T 
0  to  2.3 
1.1  to  1.2 


234  mm 
1.58  T 
70  to  100 
95 


The  characteristics  of  the  optical  cavity  used  in  this  experiment  are  the  following 


Length 

Mirror  radius  of  curvature 
Rayleigh  range 
Wavelength  of  maximum  Q 
Average  mirror  reflectivity  at  63£8  A 
Round  trip  cavity  losses  at  6328  A 
Mirror  transmission 


5.5  m 
3  m 
1  m 

620  to  680  nm 
99.965  % 

7  10-4 
3  10"5 


Optical  measurements  were  made  using  the  experimental  set-up  shown  in  figure  1. 


PIN  HOLE 


0  75  p 


PHOTOMULTIPLIER 
R  955  r. 


1  IAa/wa/ 


Figure  1.  Experimental  set-up 

This  set-up  was  originally  used  to  investigate  the  spontaneous  emission*®^.  The 
transmitted  part  of  the  synchrotron  radiation  stored  in  the  cavity  is  refocused  by  an 
aluminium  mirror  Mj  of  1.75  m  focal  distance.  A  pinhole  of  75  p  diameter  pinhole  is  placed 
at  the  mirror  focal  distance.  The  light  passing  through  the  pinhole  is  sent  to  a  mono¬ 
chromator  by  an  aluminium  mirror  M2.  Typical  monochromator  bandwith  was  0.3  A.  The  light 
coming  out  from  the  monochromator  is  detected  in  a  photomultiplier,  the  signal  being  sent 
to  a  strip-chart  recorder.  The  pinhole  is  mounted  on  a  horizontal  and  vertical  translation 
plate  for  transverse  profiles  records. 

Results 

Laser  oscillation  was  obtained  after  a  careful  alignment  of  the  electron  beam  on 
the  cavity  axis  (within  0.1  mm  all  along  the  1.3  m  O.K.  length)  and  maximization  of  emission 
as  a  function  of  the  storage  ring  radio  frequency  cavity  voltage  and  the  optical  cavity 
length.  Laser  operation  requires  very  precise  synchronism  between  light  pulse  reflections 
and  electron-bunch  revolution  frequency.  To  avoid  backlash  and  mirror  misalignment,  fine 
tuning  of  the  cavity  length  was  performed  by  slightly  changing  the  R.F.  frequency  instead 
of  translating  the  mirrors. 

Fig.  2  shows  two  "detuning  curves"  of  laser  power  (normalized  to  the  maximum)  versus 
R.F.  frequency  variation  and  equivalent  mirror  displacement. 


imwwo  we  no  us 

RF  frequency  (Hi) 


Figure  2.  Normalized  laser  power  as  a  function  of  RF  frequency  and  equivalent  mirror  dis¬ 
placement.  Curve  (a)  is  recorded  close  to  the  maximum  gain/loss  ratio  and  curve  (b)  close 
to  the  laser  threshold.  The  shift  between  the  two  curves  is  probably  due  to  some  slight 
cavity  length  drift. 

Curve  (a)  has  a  3.7  ym  F.W.H.M.  of  equivalent  mirror  displacement  ;  curve  (b)  as 
recorded  much  closer  to  laser  threshold,  has  only  al.6  ym  F.W.H.M.  The  shift  in  displa¬ 
cement  between  curves  (a)  and  (b)  is  probably  due  to  some  slow  cavity  length  drift  (e.g. 
a  temperature  drift  of  0.02°C/30  minutes  would  be  sufficient).  In  this  experiment,  the  gain 
was  not  proportional  to  the  electron  current,  mainly  due  to  the  anomalous  bunch  lengthening 
and  energy  spreading  at  high  current  ;  these  effects  make  the  gain  versus  ring  current  reach 
a  maximum  and  then  decrease.  One  consequence  of  this  is  that  the  ratio  of  gain  to  cavity 
losses  always  remained  just  above  1  for  laser  operation.  Wider  detuning  curves  are  expected 
for  higher  gain/losses  ratios. 

Fig.  3  presents  the  horizontal  and  vertical  transverse  profiles  of  the  laser  output. 

They  are  in  a  very  good  agreement  with  the  expected  cavities  TEMoo  profile.  The  slight  discre¬ 
pancies  might  arise  from  some  residual  instability  of  the  laser  too  close  to  threshold  or 
non  uniform  mirror  reflectivity. 


Figure  3.  Experimental  laser  horizontal  and  vertical  transverse  profiles  and  the  calculated 
cavity  TEM00  profile. 

Fig.  4  shows  the  laser  time  structure  in  a  200  ms  total  time  scale.  Although  not  seen  in 

:’ig- - ,  the  laser  reproduces  the  pulsed  27  7Hz  structure  of  the  electron  laser.  As  it  can  be  seen  on  fig.  A, 
the  laser  adepts  a  quasirar.icr.  pulsed  structure  over  a  long  time  scale  with  a  typical  rise  time  around  200  us 
(corresponding  to  a  gain  minus  losses  of  2.  1C-4  per  pass)  and  a  quasirar.dcm  period  of  about  20  to  AO  ms. 
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20  ms 


Figure  U.  Laser  time  structure  over  a  200  ms  interval. 

Each  subpeak  is  also  usually  more  or  less  modulated  at  13  KHz  very  dnse  to  the  theoretical  R.F.  synchrotron 
frequency.  It  should  be  noted  that  this  time  structrure  depends  on  the  electron  energy  and  also  how  far  the 

laser  is  from  threshold.  .  ,  ,  J  ^  .  .  .  ,  , 

Fig.  5  presents  two  spectra  :  a)  is  recorded  without  amplification  (optical  cavity  detuned),  b)  is  recorded 
at  laser  operation  (cavity  tuned).  Usually,  the  laser  oscillates  simultaneously  at  three  wavelength  with  a 
wrontly  dominant  central  line.  Each  wavelength  is  located  at  a  maximum  of  the  gain  versus  wavelength  curved'. 
Laser  tunability  was  obtained  from  640  nm  to  655  nm  by  changing  the  magnetic  gap.  In  fact,  the  range  of 
tunability  is  limited  for  the  moment  by  the  mirror  reflectivity. 


Figure  5,  Spectra  ot  the  cavity  output  radiation  under  two  conditions  ;  (a)  cavity  detuned  (no  amplification), 
(b)  cavity  tuned  (laser  on). 

Fig.  5  shows  an  enlargement  of  the  main  laser  line  spectrum  recorded  by  using  a  one  dimensional  C.C.D. 
detector  instead  of  the  usual  monochromator  exit  slit.  The  aperture  tiige  is  3ms.  Each  narrow  square  peak  in  this 
curve  is  recorded  on  only  one  C.C.D.  element  and  corresponds  to  a  0.3  A  spectral  width  of  the  same  order  as 
monochromator  resolution.  We  conclude  thatthereisa  residual  inhomogeneous  contribution  to  the  laser  line- 
-width  probably  connected  with  tjje  long  time  scale  laser  pulsed  structure  (see  fig.  4 ) ,  so  that  typical  laser 
lines  are  gaussian  with  about  3  A  F.W.H.F!.  if  averaged  over  a  time  scale  of  =  1  sec. 

The  central  wavelength  of  any  line  is  always  equal  to  the  wavelength  of  maximum  emission  of  spontaneous 
emission  with  the  cavity  completely  detuned  (no  amplification)  plus  0.15  of  the  wavelength  ;ihterfringe 
distance  (see  fig.  5  )  instead  of  0.25  as  predicted  from  Madey's  theorems).  This  discrepancy  is  probably  due 
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to  the  transverse  multimode  content  of  the  snontaneous  emission  stored  in  the  cavitv.  wheareas  the  measured 
transverse  rrofiles  of  the  laser  cutout  show  that  laser  operation  is  only  achieved  on  the  TE.”  oo  mode. 


A  typical  75  uW  average  output  power  has  been  recorded  at  50  mA  current  of  166  JfeV  electrons.  This 
corresponds  to  a  typical  6C  mW  output  peak  power  over  the  1  ns  electron  bunch  length  (including  effects  of 
both  long  and  short,  time  scale  pulsed  structure  of  the  laser),  and  a  2  kW  intracavity  peak  power.  This  gives 
a  2.4  x  10-2  efficiency  with  respect  to  the  3.1  watt  ot  total  synchrotron  radiation  power  generated  in  the 
whole  storage  ring  at  this  energy  and  current.  This  efficiency  is  about  0.4  the  Renieri  typical  efficiency 
(10)  . 


c 

n  -  —  x  rwir  ''•here  ^ 


mirror  transmission _ 

round  trip  cavity  losses 


and  at„E  =  i'?  10"3  is  the  R.M.S.  relative  energy  spread.  In  this  experiment  hmir  -  0.043. 

CL 

This  low  mirror  efficiency  is  mainly  due  to  the  high  absorption  in  the  mirror  dielectric  layers  compared  to 
its  transmission.  The  rather  low  Renieri  average  output  power  is  therefore  the  result  of  high  sensitivity  of 
the  O.K.  to  energy  spread,  poor  mirror  efficiency,  and  weak  total  synchrotron  radiation  power  at  166  JfeV 
(power  is  proportional  to  the  fourth  power  of  the  electron  energy).  As  shorn  in  Fig.  6,  the  laser  output 
power  decreases  with  the  electron  beam  current  and,  except  in  regions  too  close  to  threshold  is  usually  pro¬ 
portional  to  the  current  as  predicted  by  the  same  model. 


Figure  6.  Laser  power  together  with  electron  beam  current  versus,  time. 

However,  important  discrepancies  appear  in  the  structure  for  long  time  scales  (see  Fig.  3)  which  is 
predicted  to  be  constant.  Moreover,  bunch  shortening  instead  of  bunch  lengthening  has  been  seen  at  laser 
turn  on  :  such  an  effect  lias  already  been  seen  at  much  lower  current  by  stimulating  the  F.E.L.  with  an  argon 
laser 


Conclusions 


In  summary,  this  work  demonstrates  for  the  first  time  the  feasability  of  the  storage  ring  Free-Electror. 
Laser.  Future  experiments  will  continue  to  analyse  the  saturation  mechanism  and  the  long  and  short  tin?  scale 
structure.  Higher  gair./losses  ratios  are  expected  from  the  use  of  better  mirrors,  a  second  RF  cavity,  and 
positrons  instead  of  electrons. 
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